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КВАНТОВА ФІЗИКА ЖИВОГО TA ЇЇ 
ТЕОРЕТИЧНІ ОСНОВИ 


ПАРАДИГМА ПОВЗДОВЖНІХ ЕЛЕКТРОМАГНІТНИХ 
ХВИЛЬ В КОНЦЕПТУАЛЬНИХ ОСНОВАХ ФІЗИКИ 
ЖИВОГО: КРИТЕРІЇ ТА УМОВИ ФОРМОУТВОРЕННЯ, 
ІНФОРМАЦІЙНО-ДІАГНОСТИЧНА ЗНАЧИМІСТЬ 
Ю.В. ЧОВНЮК, Т.М. ОВСЯННІКОВА, Б.Ф.РУДЬКО 


Науково-дослідний Центр квантової медицини "Відгук" МОЗ України 


PARADIGM OF LONGITUDINAL ELECTROMAGNETIC 
WAVES WITHIN THE FRAMEWORK OF PHYSICS OF 
THE ALIVE CONCEPT: CRITERIA AND CONDITIONS 
FORMATION, INFORMATION-DIAGNOSTICAL VALUE 


Yu.V. CHOVNYUK, T.N. OVSYANNIKOVA, B.F. RYUD’KO 
Scientific-research Centre of Quantum Medicine “Vidhuk” 


Abstract. The objective of the given paper consists in substantiation of the paradigm 
of longitudinal electromagnetic waves within the framework of physics of the alive 
concept, evaluation and determination of the criteria of such waves formation in 
bioobjects, their “physical nature”, definition of informative and diagnostical value of 
the above-mentioned electromagnetic waves for the analysis of the main characteris- 
tics of the living matter. 

Ключові слова: повздовжні електромагнітні хвилі, фізика живого, діагностика 


У [1] обговорюється гіпотеза про існування повздовжніх електромагніт- 
них хвиль у будь-яких середовищах (у т.ч. у живій матерії). Проблема 
повздовжніх електромагнітних хвиль спочатку вперше була сформульова- 
на в квантовій електродинаміці. Теорія В. Фока і Б. Подільського, відповід- 
но до якої в лагранжиан електромагнітного поля вводився "член, що фіксує 
калібрування", привела до розширення квантової електродинаміки, тому 
що дозволила погодити теорію електромагнітного поля з перестановочни- 
ми співвідношеннями. Крім того, у такій “модифікованій” квантовій елект- 
родинаміці з'явився новий екзотичний (із погляду фізики) об'єкт -- 
повздовжні електромагнітні хвилі, що безумовно, явно суперечили макс- 
велловскому уявленню про винятковість поперечності електромагнітних 
хвиль. Розроблена в квантовій електродинаміці модель фізичного вакууму, 
що є плазмою віртуальних електронів і протонів, явно суперечила цьому 
уявленню. Добре відомо, що в плазмі можуть поширюватися як поперечні, 
так і повздовжні (ленгмюрівскі) хвилі, отже, вакуум квантової електродина- 
міки, як середовище поширення, узагалі говорячи, допускає збудження в 
ньому повздовжніх хвиль. У рамках теорії поперечних електромагнітних 
хвиль повздовжні хвилі є “нефізичними”, хоча вони і входять в обов'язково- 
му порядку в математичний апарат квантової електродинаміки. Автор [1] 
приводить систему доказів, що спростовують нефізичність повздовжніх 
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хвиль на основі не зв'язаного з теорією Максвелла аналізі релятивістської i 
каліброваної інваріантості, самопогоджуванності теорії, а також встановив 
матеріальне джерело таких хвиль (фізичний вакуум). Використовуючи 
нетрадиційний підхід (заснований не на лагранжевому формалізмі), 
Н.П. Хворостенко виявив взаємозв'язок між рівняннями Максвелла і Діра- 
ка 1 виключила в спінорному уявленні електромагнітних полів умову їхній 
поперечності. Головна перевага таких хвиль - - це відсутність їхньої взаємо- 
дії з електричними зарядами. 

Теоретичне й експериментальне дослідження повздовжніх електромаг- 
нітних хвиль було продовжено в роботах тульської наукової школи [2,3]. 
Відповідно до зазначених досліджень рух тіла по інерції (наприклад, атома 
речовини, ядра атома іт.д.) можна розглядати як результат функціонування 
системи самопідтримки руху, у якій інертність ядер і електронів, що склада- 
ють конкретне тіло, обумовлюється бездиссипативним обміном енергією з 
фізичним вакуумом (ефіром). Даний процес носить складний коливальний 
характер і супроводжується створенням в окрузі частки плинів, що рухаєть- 
ся, ефіру. При цьому потік вакууму, що рухається, викликаний переміщен- 
ням ядра, підтримує існування руху тіла. 

Припинення переміщення об'єкта (за рахунок зовнішнього впливу) не- 
минуче приведе до трансформації енергії ефірного плину в приядерній об- 
ласті в енергію деякого хвилястого процесу у фізичному вакуумі. Саме да- 
ний хвилястий процес, на думку учених тульської школи, відповідає за ме- 
ханізм передачі енергії в ланцюжку атомів, інший упорядкованій атомній 
(молекулярній) структурі. При цьому існує два різновиди енергообміну між 
елементарними об'ємами (атомами): “повздовжній” і “поперечний”. Kox- 
ний із зазначених варіантів енергообміну має свою власну швидкість поши- 
рення. Так, у випадку "повздовжнього" енергообміну швидкість його пере- 
дачі U, близька до швидкості повздовжньої магнітної хвилі, що є торсіон- 
ним полем (полем крутіння). По оцінках зробленим у [2,3], U, є величиною, 
що значно перевищує швидкість світла у вакуумі. При "поперечному" енер- 
гообміні швидкість його передача U, близька до швидкості поширення 
повздовжньої електричної хвилі. По оцінках тих же авторів U, складає 
(15. 22) :10* м/сек. Аналіз отриманих експериментальних і теоретичних 
результатів представників тульської наукової школи дозволив iM висунути 
гіпотезу про те, що жива речовина здатна до генерації повздовжніх елект- 
ромагнітних хвиль [2]. Тому, по їхньому твердженню, біооб'єкти мають по- 
тенційну можливість за допомогою повздовжніх випромінювань, а у випад- 
ку створення синхронізованих кооперативних повздовжніх сигналів — 
спроможність управляти переміщенням організму як цілого. 

Ціллю цієї роботи є обгрунтування парадигми повздовжніх електромаг- 
нітних хвиль у коцепції фізики живого, проведення оцінок і встановлення 
критеріїв формоутворення подібних хвиль у біооб'єктах, їх “фізичності”, 
визначення інформаційної і діагностичної значимості таких електромагніт- 
них хвиль для аналізу основних характеристик живої матерії. Для реалізації 
зазначених цілей роботи використані підходи [4,5], справедливі в рамках 
максвелловскої теорії електромагнітного поля, засновані на аналізі діелект- 
ричної проникності біооб'єкта, а точніше її просторової дисперсії в області 
слабкого поглинання, де узагалі говорячи, може проходити через нуль. 
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У рамках моделі провідного середовища |5) поширення в ній хвиль іс- 
тотно відрізняється від поширення в середовищі, що не проводить. Якщо 
Ha ряду з діелектричною і магнітною проникностями є і и средовище 
характеризується так само провідністю о, то до рівнянь Максвелла варто 
додати закон Ома: 


ј=о:Ё (1) 


З урахуванням (1) можна представити рівняння Максвелла у виді (вико- 


ристана гауссова система одиниць): 
4 (в) = 0; rotE +". 28 
° (2) 


У діелектрику, що He проводить, змінними в час! поля є поперечні [5], 
тобто вектори Е i Н перпендикулярні напрямку просторової зміни поля. У 
граничному випадку нульової частоти, як відомо з електро- і магнітостати- 
ки, статичні поля в діелектрику повздовжні в тому розумінні, що вони утво- 
рюються зі скалярних потенціалів і тому паралельні напрямку просторової 
зміни останніх (це особливо стосується тих типів електромагнітних хвиль, 
що поширюються в живій матерії, що мають низькі частоти (порядку декіл- 
ькох герц)). 

При відмінній від нуля провідності, котра безумовно має місце в біо- 
об'єктах (у т.ч. біоклітин тканин організму людини й ін.), картина карди- 
нально змінюється . Розглянемо для простоти поля залежні тільки від однієї 
просторової координати 5. Разкладемо E i Н на повздовжню і поперечну 
складові: 


E(t) - E,(&7) +E, (&1), 


Е : - (3) 
H(&t)=H,(&t)+H, (51). 

З властивості оператора ротора випливає, що поперечні складові E i H 
задовольняють двом роторним рівнянням системи (2), що приводить до по- 
перечних хвиль, у той час як повздовжні складові задовольняють 
рівнянням: 


дНи дНи я 
дЕ ' 8t і 
дЕ, д 4no (4) 
=0, Е «бе в, = 0. 
[oa ot = 


З першої пари рівнянь випливає, що єдиним можливим говздовжнім 
магнітним полемо є статичне однорідне поле. У цьому відногтс лічі провідне 
середовище поводиться так як і те, що не проводить. 
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Однак друга пара рівнянь у (4) показує, що повздовжнє електричне по- 
ле однорідне в просторі, але змінюється в часі: 


E,(E,t)=E,-e * (5) 


Отже, у відсутності зовнішніх струмів статичне повздовжнє поле не 
може існувати в провідному середовищі. Для гарних провідників, таких, як 
мідь, с має розмір порядку 10" сек", так, що коливання загасають за дуже 
короткий час. У той же самий час с біооб'єктів (біоклітин) порядку 
(10"...10") сек", тому загасання повздовжніх електромагнітних хвиль 
відбувається з характерною частотою: 

4nc 
ъ= (6) 
складові десятки гігагерц, близької один по одному до f, =0,, T = 60 ГГц, 
зазвичай використовуваної в методі мікрохвильової резонансної терапії 
(МРТ) професора С.П. Сітька. 

Розглянемо тепер біосередовище в рамках, що слабопоглинає 14), при- 
чому нехай область слабкого поглинання простирається в деякому широко- 
му інтервалі частот від о, до о, (причому о, >> o), i припустимо, що o такі, 
що O, << @<< @,. Останнє припущення в біооб'єктах цілком виправдане, якщо 
вважати, що (9/27) = 60 ГГц, (о, /2n) = (1. 10)Гц, (о, /2n) = (10°...10° ГГц, 
Відповідно до співвідношень Крамерса-Кроніга (1927) [4] для комплексного 
уявлення діелектричної проникності середовища e(o) (e(o) == (©) +ie "(o), 
де e'(@)Ree(o), é"(@) = Ime(o),i = +/-1), можна одержати: 


e'() -l= В г ах, (7) 


х-ф 
А 1) *ze'(x)-1 
s*(a) -1 = 2 s dx, (8) 


де знак інтеграла означає, що інтеграл від полюсного вираження розуміється 
взмісті його головного значення. У зв'язку зтим, що функціяє (сунепарна [4], 
формулу (7) можна привести до виду: 


2 х= "(х)ах 
ст (9) 


о? 


Область інтегрування в (9) легко розбити на дві частини: x « o, іх» Q,. 
Ця розбивка дозволяє в першій частині інтеграла (області інтегрування) 
знехтувати в знаменнику підінтегрального виразу х у порівнянні 3 ©, a 
другий — о y порівнянні 3 x: 
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є"(о) =1+ й ‘fe"(x) ax - ES ЕО: (10) 


x 
x то 
тобто функція e'(@) у розглянутій області має вигляд a-b/o, ne a i b — 


позитивні постійні, обумовлені з наступних співвідношень: 
A 7r gay г 
а= 2 fe (x) = 
о x; 


2 gi (11) 
b= = i "(х)ах 


Другу з зазначених констант (b) можна виразити через силу осциллято- 
рів N,, відповідальних за поглинання в області від 0 до o, [4], і тоді: 


є"(о) заз Ме" (12) 


то? 


де e, m — заряд i Maca осциллятора (електрона). 

Проаналізуємо цей вираз (12). Видно, що в досить широкій області 
слабкого поглинання діелектрична проникність e(o), а точніше її дійсна 
складова є"(о), узагалі говорячи, проходить через нуль. Потрібно в цьому 
зв'язку нагадати, що прозорі в буквальному значенні слова є біосередови- 
ща, у котрих e(o) не тільки дійсне (c(@) =='(о)), але і позитивне; при He- 
гативному значенні «(о хвиля загасае в глиб біосередовища, хоча в ній і не 
відбувається диссипації енергії. 

Якщо частота o, при котрій &=0 (тобто o, = Ма), "задіяна" у спектрі 


електромагнітних хвиль, що взаємодіють із біооб'єктом, то індукція D то- 
тожно обертається в нуль у такому біосередовищі і рівняння Максвелла 
допускають існування змінного електричного поля, що задовольняє усього 
лише одному рівнянню: 


rotE =0 (13) 


при рівному нулю магнітному полі. Іншими словами, у цій ситуації біосередо- 
вище "підтримує" існування повздовжніх електричних хвиль. Для визначення 
швидкості їхнього поширення необхідно враховувати дисперсію діелектрич- 
ної проникності биосередовища не тільки по частоті, але і по хвильовому 
векторук. 

Розглянутий підхід дозволяє обгрунтувати з фізичної точки зору резо- 
нансний ефект методу МРТ професора Сітько С.П. відповідно до зазначе- 
ної методики, вплив на біооб'єкт (у БАТ шкіри людини) повинен робитися 
на деякій характерній частоті (f, 21-60 ГГц). Це означає, що в широкій 
області прозорості біооб'єкта є вузька область ("лінія") поглинання навко- 
ло деякої частоти oy. Розглянемо детальніше округу цієї частоти, що задо- 
вольняє умові: 
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у << |o - o| << oy (14) 


де y — ширина лінії. У цій частотній області в підінтегральному виразі в (7) 
можна замінити х на 0), скрізь, крім швидко змінної функції e"(x). Toxi 


одержимо: 


(Que re (15 


c 


де інтегрування робиться по лінії поглинання (C):(@,-Y, @,+Y). Добротність 
електромагнітних коливань, що ефективно поглинаються в біосередовищі, 
може бути оцінена співвідношенням: 


E (16) 


Якщо 0-50, TO €'(0)—0 за умови, що £"(x)0 в област! Ee [@,—y,0,+y], тоб- 
TO €”(@)~0(@,—@), (нуль більш високого порядку, ніж (0,-0)). У цій ситуації 
(€’(@,), €”(@,) )-з0 і у-з0, тобто добротність таких електромагнітних коливань 
істотно зростає, що може свідчити на користь ефективності впливу 
повздовжніх електромагнітних хвиль на біосередовище (Q~10°). 

Послідовне формулювання умов, котрі передумовлюють зв'язок між 
частотою і хвилевим вектором (т.3. закон дисперсії) цих повздовжніх елек- 
тромагнітних хвиль (ПЕМХ), потребує врахування просторової дисперсії 


діелектричної проникності біосередовища є. Ця умова така: 
є (о, к) «0, (17) 
де 1ндекс | означае повздовжню (HEMX) хвилю. При малих значення хвиле- 
вого вектора К розв'язок рівняння (17) мае вигляд: 
-— 
о(к) = 01, жоК (18) 
где о. — константа, а о; — значення частоти, що обертає в нуль проникність 


€(@)=€,(0,0). При цьому швидкість розповсюдження ПЕМХ визначається Ta- 
ким чином: 


|= 
ug 
І 


U = =a-k 
ok (19) 


i пропорційна першому стушню хвилевого вектора К. Швидкість таких 
ПЕМХ може бути достатньо малою (суттєво менша швидкості світла в цьому 
біологічному середовищі). 

В рамках квазікристалічної моделі біотіла можна розглядати просторо- 
ву дисперсію біосередовища як оптично не активну. Тоді діелектрична про- 
никненість такого середовища стає тензором другого рангу ғ, (с, К), а по 


— 


ступіням k потрібно провести розклад у ряд і врахувати, як мінімум, 
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квадратичні члени. Для подальшого використання зручніше писати цей 
розклад для зворотнього тензора т, = ez , (i,k) = (1,2,3). Запишемо його y 
вигляді: 


Nik = n (w) + Pikim (W)K Kn (2, т) E (123) (20) 


Тензор В, (четвертого рангу) можна вважати симетричним по другій 
парі індексів, поскільки він домножується на симетричний добуток КК . 
Враховуючи симметрію тензора €,,(@,k) в магнітному полі В: 


вл (в; B) = eu (8-k;-2), (21) 


присутньому в розглянутому біосередовищі, можливо встановити симметрию 
тензора Ві при В «0: 


Bu, = Ва» = В emi (22) 


Він зовсім не повинний бути симетричним по відношенню до переста- 
новки обох пар індексів. У відсутності поглинання тензор Ві» дійсний. 

В ізотропному біосередовищі тензор В, повинний відтворюватися 
тільки через одиничний тензор, т.ч. він має вигляд: 


Ваһ = Відждіт «P (6,6 +845 im (23) 


де 5, i йому подібні — символи одиничного тензора, причому тензор Pium мае 
тільки дві незалежні компоненти. В ізотропному біотілі (в рамках кристаліч- 


ної моделі) також і n = 8, i, таким чином, тензор (20) приймае вигляд: 
= (30) «ge? kk 
"а =| n +B, 04 + BK; k,, (24) 


B відповідності з загальним виглядом діелектричного тензора B 1зотропному 
середовищі з просторовою дисперсією: 


== КК ‚К 
e, (o. K) EODD = a)r (25) 


B (25) є, — означає компоненти діелектричної проникності середовища 
для поперечних і повздовжніх електроматнітних хвиль, відповідно. Вони 
залежать тільки від абсолютної величини хвильового векторах i Bit о. Якщо 
напруженність електричного поля Е направлена по хвилевому вектору, то 
індукція D ==, Е; якщо ж Elk, ro  - є, Є. Відповідно величини є, іє, прий- 
HATO називати повздовжньою i поперечною діелектричними проникностя- 
ми (біосередовища). При К — 0 вираз (26) повинний наближатися до зна- 
чення є(с2)9,, що не залежить від напрямку k; при цьому зрозуміло, що: 


є (о, 0) = є (о, 0) = e(o) (26) 
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Розповсюдження електромагнітних хвиль в біосередовищі визначаєть- 
ся рівняннями вигляду: 
-à -1 EN, EN Кб: 35 
(п 5 up a n =, Пэн, В, =e Des Па (27) 
Але при підстановці (23) y ui рівняння анізотропний член з 8, випадає за 
причини ортогональності векторів Dik в плоскій хвилі (поперечній елек- 
тромагнітній), т.ч. біосередовище залишається оптично ізотропним. В ку- 
бічних квазікристалах (модель біосередовища) тензор B, не зводиться до 
одиничного тензору; в залежності від кристалічного класу він має для цих 
квазікристалів 3 або 4 незалежні компоненти. Без урахування просторової 
дисперсії кубічні квазікристали оптично ізотропні; врахування квадратич- 
ної по К дисперсії приводить до появи у них нової властивості — оптичної 
анізотропії, що було встановлено Лорентцом іще в 1878 році. 
. . (0) _ / 0) . Xs : 
В кубічному кристалі n; ^ 28, /є" ^ і розклад (20) приймає вигляд: 


1 
Nik = NO! Sic + Dank, (28) 


0 . ds 
Вводячи; = c a pics X, декартової системи координат х,, X,, X, направ- 
ляючи вздовж хвилевого вектора К, можна отримати: 


ел т ^g Pas D, =0 (29) 


Ці рівняння схожі з рівняннями для електромагнітних хвиль, що розпов- 
сюджуються в некубічних квазікристалах без урахування просторової дис- 
персії. Їх визначник представляє собою квадратне по п” рівняння, що виз- 
начає показники заломлення двох хвиль з однаковим напрямком i різними 
поляризаціями. Таким чином, просторова дисперсія в кубічному ква- 
зікристалі, що моделює біосередовище, знімає виродження по поляризації 
-- швидкості двох хвиль стають різними і залежними від напрямку 
(причому не має значення які з компонент t(o,k) дають вклад ви’ — 


повздовжня і поперечна є, є, — незалежно, чи дві одночасно). 
Аналіз дисперсійних співвідношень для ПЕМХ (17), (18) для різних 
знаків о показує, що при o>0 групова швидкість цих хвиль до /9к > 0, тобто 


енергія в ПЕМХ переноситься "вперед" (вздовж напрямку k), в випадку 
0<0, группова швидкість від'ємна i виникають "зворотні" ПЕМХ, в яких 
енергія в хвилі передається в напрямку, протилежному K. 

Відмітимо, що співвідношення (18) справедливі для випадку, коли: 


Qa 


1 
ak? ~ op k~ |.) =k, (30) 


Якщо К << k,, то розклад (18) стає невірним і потрібно враховувати 
наступні (по k) члени розкладу в (18), о(К), при яких виконується (17). 
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Розглянимо середовище в ізотропному наближенні, що немає природ- 
ної оптичної активності (тобто гіротропне). Оскільки ізотропне середови- 
ще залишається оптично ізотропним і при урахуванні просторової диспер- 
сії, це означає, що закон дисперсії ПЕМХ в такому середовищі дається 
рівнянням: 


n? ==, (в, К) (31) 


Причому під ¢(,k) слід розуміти повздовжню складову діелектричної 
проникності (тензора є, (09) (25)) є, Оскільки &(®,К) задовольняє тим же 
співвідношенням Крамерса-Кроніга, що і функція e(o) без просторової 
дисперсії, можна стверджувати, що поблизу лінії поглинання функція 
€,(@,k) має вигляд: 


e(o) = є (©) = m FOS (32) 


де Є — контур лінії поглинання, причому у << o - o |«« o», a ү ширина цієї 


лінії, 0, — деяка частота, поблизу якої відбувається резонансне поглинання 
електромагнітної енергії (наприклад, НЗВЧ-діапазону електромагнітних 
хвиль, що застосовуються в методі МРТ професора С.П. Сітько), 
е; (о) ==! (©) чіє"(о),і = МЕТ 

Представимо співвідношення (32) в іншому вигляді, враховуючи, що 
сталі залежать від К (при урахуванні ефекту просторової дисперсії); запи- 
шемо це співвідношення у вигляді: 


s (o, К) sal о A(k) > 0 (33) 


[© (к) - о] 


Поблизу нуля (17) є, (o,k) величина А відносно мала (це можуть бути i 
якісь наближені правила відбору, зменшуючі матричні елементи, що визна- 
чають значення А). Тоді потрібно врахувати, поряд з полюсним членом, 
також i постійний (незалежний від ©) член B &; його можливо позначити a(k) 
і вважати, що a(k)>0, тобто на далекій відстані від лінії поглинання 
біосередовище оптично прозоре (В ПЕМХ). Справедливість наведеного 
вище припущення, тобто одночасного урахування постійного і полюсного 
членів в €,(@,k), потребує, щоб, вони були близькими одне до одного в 
області частот о-о,| << Фу, де тільки і може бути застосований полюсний 
вираз; іншими словами, повинно бути: 


А << ао, (34) 


Оскільки К вважаємо малим, можна в рамках (18), (30) розкласти функ- 
цію по його (хвилевого вектора) ступеням. При цьому досить замінити a(k) i 
A(k) сталими значеннями а = a(0) > 0, А = 4(0) > 0, залишивши поправочний 
член тільки в розкладі функції о, (К) (о, =@,(0)) в малій різниці 0,-0: 


a(k) афу +vk*, v= (1/2)a (35) 
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Тоді проникність: 


A 


e, (o, k) =а+— —— 
2 
00 +vk* -0 


(36) 


Цей вираз як функція частоти проходить через нуль в області o, (К) < o 
тобто: 
A 


e,(0,k) -Q 0-0, + vk? je (37) 


При k=0 це виникає о=0,=00+А/а. При k=0 діелектричні проникності 
є, і є, співпадають, саме цьому ©, Є гранична (при к-»0) частота 
повздовжньої електромагнітної хвилі. Корень функції €(@,k) прик є 0 (37) 
має фізичний зміст для "прямої" і "зворотньої" ПЕМХ). Розглянемо (31) з 
урахуванням (36). Тоді дисперсійне рівняння (31) приймає для ПЕМХ 
біосередовища такий вигляд: 


2 2 
(пі ви -очо,) = A, (38) 
vo^ lao! ve, _ laa 
де введено позначення f = — =——— а —0 =L ця величина може 
c 2 с? c 2c 


мати два знаки (В>0 для “прямих” ПЕМХ (ППЕМХ), В<0 для “зворотніх” 
ПЕМХ (ЗПЕМХ)). Рішення цього рівняння можна розглядати як результат 
перехрещення двох гілок спектра - - звичайної світової (КВЧ-діапазону) елек- 
тромагнітної хвилі з и; і ПЕМХ зи? = (o — o, )/B, зв'язаної з полюсом діелек- 
тричної проникності e(o,k). Потрібно підкреслити, що в мікроскопічній теорії 
полюс діелектричної проникності відображає існування бозевських елемен- 
тарних збуджень в біодіелектрику — ексітонів; знак В — сталої при цьому, 
співпадає зі знаком ефективної маси ексітона. Відповідну гілку спектра хвиль 
(ПЕМХ) потрібно називати ексітонною. Область спектру поблизу віртуаль- 
ного перехрещення обох гілок називають поляритонною. 

Розповсюджуваним в біосередовищі хвилям відповідають тільки додат- 
ні корені nj. Потрібно відмітити, що суто уявні значення п (п) « 0) відпо- 
відають хвилям, для яких біосередовище не прозоре (хоч в ній і немає 
поглинання), тобто можна сказати, що хвиля в цьому випадку відбивається 
від біосередовища. При В» 0 крива проникає в область о> о,, тим самим 
утворюючи тут допоміжну хвилю, котрої не було б без урахування просто- 
рової дисперсії; при o > o, в біосередовищі можливе розповсюдження двох 
різних ППЕМХ. Точка, в якій dn; (о)/ао =+% і два корені зливаються, 
маємо два додатні корені і в біосередовищі можуть розповсюджуватись дві 
різні ЗПЕМХ: 


В Om — Dy = -pl a+2 В (39) 


Те ж саме відноситься до області між точкою m' 3 координатами: 
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= ок -o = ff 2-2/5 (40) 


2 = 
Ni,» =a + 


і точкою O=O: 

Якщо величина А не достатньо мала, збереження постійного члена в 
(36), строго кажучи не послідовне. Опустивши цей член (тобто прирівняв- 
ши в цих формуах a=0), приходимо до картини, що горизонтальною асимп- 
тотою віх кривих стає вісь абсцис (замість лінії п? = а). Область O20, випа- 


дає з розгляду. 

Зі сказанного можна зробити такі висновки: 

1.При розгляді прозорих біосередовищ, що підтримують повздовжні 
електромагнітні хвилі з урахуванням просторової дисперсії діелектричної 
проникності (повздовжньої — &), в залежності від знаку останньої (sign) 
виникає два різновиди ПЕМХ: К 

а) прямі (EMX), що переносять енергію вздовж хвилевого векторах, 
у яких групова швидкість V, = 0a/dk < 0; 

б) зворотні (ЗПЕМХ), що переносять енергію в напрямку, зворотньому 
хвилевому вектору К, у яких групова швидкість V, = ook > 0. 


2. ППЕМХ, ЗПЕМХ надзвичайно перспективні для використання і 
пояснення ряду ефектів резонансного поглинання енергії біооб'єктами, 
при впливі на останні електромагнітним полем в діапазоні КВЧ (несучої 
частота f,~60 ГГц), що застосовується в методі MPT проф. С.П. Сітько. 
Головна особливість ПЕМХ полягає в тому, що вони можуть невзаємодіяти 
з зарядами біотіла, а значить проникати в глибину біотіла на відстані, що 
набагато перевищують товщину скін-шару. Окрім того, швидкість ПЕМХ 
може суттєво відрізнятися від швидкості світла с в біосередовищі, що 
опромінюється НЗВЧ-полем, інколи може бути набагато менше останньої. 
Знак і величина групової швидкості ПЕМХ залежать від дисперсії 
повздовжньої діелектричної проникності не тільки частотної, але й 
просторової. Наявність "повільних" ПЕМХ, що підтримуються біосередо- 
вищем, дозволяє їм ефективно взаємодіяти з біосередовищем і підвищувати 
до оптимальних основні показники і характеристики НЗВЧ-електромагніт- 
ного поля, впливаючого на організм людини (т.е. лікувальний ефект). 

Розглянемо ситуацію поблизу лінії поглинання (частота ї,= 0/27~60 ГГц) 
в ізотропному біосередовищі. Будемо вважати, що ізотропність біосередо- 
вища обумовлена двома фізичними механізмами: 1) наявністю постійного 
магнітного поля Н (планети Земля), внаслідок чого тензор діелектричної 
проникності біосередовища перестане бути симетричним і залежить від 
цього поля: 


ex З вк (в б Ë) (41) 


2) наявністю слабкої просторової дисперсії біосередовища, внаслідок 
чого є, можна розкласти по степеням К. Як правило для конденсованих 
діелектричних середовищ це є розклад по ступіням a/A, де а — атомні 
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розміри, A — довжина хвилі КВЧ поля. 3 точністю до членів першого 
порядку малості такий розклад має вигляд: 


Ek (о, k) = = (а) th gk; із NER (42) 


ik 
(при 90 компоненти цього тензора, не пов'язані з будь-яким розкладом по о 
(просторова дисперсія є, враховується поряд з частотною), прямують до 
постійних значень). Відповідна проникність (42) зв'язок між вектором елект- 
ричної індукції D i E в монохроматичному (~e™ ) полі в координатному 
зображенні має вигляд: 


0 "m" 
де eO =e i (6,0), а Yu — деякий тензор третього рангу, залежний від частоти 


ôE, 


1 


D, = 2078. + Yin * (43) 


i 


При наявності просторової дисперсії діелектрична проникність 
біосередовища є тензором, (а не скаляром) навіть в ізотропному 
біосередовищі (яку і надалі будемо розглядати): виділений напрямок ство- 
рюється хвилевим вектором. Оскільки припускається наявність у біосере- 
довищі природної оптичної активності або природної гіротропії, то це озна- 
чає, що в біосередовищі відсутній центр симетрії. Добре відомо, що в живій 
матерії присутні оптично активні рідини, що складаються із речовин з дво- 
ма стереоізометричними формами (цукор), причому обидва ізомери повин- 
ні знаходитися в рідині в різних кількостях; така рідина не має центру 
симетрії. 

В ізотропному біосередовищі співвідношення (43) можливо представи- 
ти у вигляді: 


ENG 
Eik = Е 


+i Е Yip (44) 


де хвилевий вектор плоскої хвилі представлений у вигляді k = o /с п. Вводячи 


замість антисиметричного тензору другого рангу у, п, дуальний йому аксиаль- 
ний вектор гірації g: | 


O 
e Yu; = Cis 1 (45) 


де Cu — одиничний антисиметричний тензор третього рангу, можна є, запи- 
сати у вигляді: 
a Oh s 
tg = Е +181 (46) 
Запис (46) формально співпадає з таким для опису магнітооптичних 
ефектів (тобто залежності є, від Н). Різниця полягає в тому, що в останній 
залежності вектор g залежав тільки від властивостей біосередовища 


(i прикладеного магнітного поля Н), між тим як в (46) вектор гірації зале- 
жить 1 від хвилевого вектора поля. В відповідності з (45) компоненти 
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вектора Є в оптично активному біосередовищі являються лінійними функ- 
ціями компонент M: 


8; = Syll (47) 


Підставляючи (47) в (45), знайдемо: 
о ==. 
Р Ум = CikmE mi > (48) 
звідки зважаючи на довільності и: 


Ф 
УЗ Ум = бт ть (49) 


чим встановлюеться зв’язок мж компонентами 1стиного тензора тре- 
тього рангу Yy i псевдотензора другого рангу gy 
В ізотропному біотілі псевдотензор у, зводиться до псевдоскаляру: 


якою Yu = NL (50) 


де д, — символ Кронекера. 

Відмітимо, що псевдоскаляр є величина, що змінює знак при інверсії ко- 
ординат. Дві стереоізомерні рідини, що входять до складу біосередовища, 
формально переходять одне в одне в результаті операції інверсії; тому їх 
значення величини f мають протилежні знаки. 

Загальний вигляд тензора діелектричної проникності біосередовища 
при наявності в ній просторово-частотної дисперсії, природної оптичної 
активності і магнітооптичних ефектів (за рахунок магнітного поля Н) мож- 
на записати так: 


ca (6, й) = (ов "ya, ak 


k kk 
ik i^*k ; 
; * (о, k, H) p + ieu, (51) 

Магнітооптичний ефект B ізотропних біотілах представляє інтерес з 
приводу свого своєрідного характеру і порівняно великої величини. Залеж- 
ність Mix Di E при урахуванні цього ефекту подається формулами (не беру- 
чи до уваги величини другого порядку малості): 


Е--5+Бхб] D -sE «i[E 3], (52) 


де = — діелектрична проникненість ізотропного біосередовища у відсутності 
магнітного поля Н, С — вектор оптичної активності, котрий в тому ж набли- 
женні пов'язаний з $ співвідношенням: 
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] x 
Со ння (53) 
= 
Залежність g (або С) від зовнішнього магнітного поля зводиться в ізо- 


тропному біосередовищі (з урахуванням того, що Н порівняно слабке поле) 
до слідуючого: 


à = Я, | (54) 


де f — скалярна постійна, котра може мати будь-який знак (плюсчи мінус). В 


той же час, у оптично активного ізотропного біотіла вектор гірації g=fn,i 
зв'язок між індукцією i напруженністю електричного поля хвилі подається 
формулою: 


Dac Е + if[E x я]. (55) 


Відмітимо, що в (51) є, іє, залежать тільки від абсолютної величини хви- 
левого вектора (i від œ, Н ): Якщо напруженність E направлена по хвилево- 
мувектору (повздовжні електромагнітні хвилі), то індукція D == ‚Е; якщо ж 


ЕЕ, To D -єб. Відповідно, величини є, i €, називають повздовжньою і попе- 
речною проникностями. Прик -» 0 ) є, (б, К іє, (о, Kk) прямує до значення, не 
залежного від напрямку k, але від Н залежність може бути: 


є (о, я) = s, (з, Я) = (00, Я). (56) 


Введемо п, = Ve — це коефіцієнт заломлення б1осередовища у відсут- 


ності поля Ні прик -» 0. 

Зміна показника заломлення біосередовища при урахуванні природної 
активності речовини і магнітооптичних ефектів є малою величиною. Тому 
при ii визначенні можна покласти в малому члені [E х 8] в (55) ів другому 


рівнянні (52) пен, = є"), Тоді задача про розрахунок різниці п-п, фор- 


мально співпадає як для випадку розрахунку зміни п під впливом магнітно- 
го поля так і для розрахунку змін п під впливом ефектів природної оптичної 
активності. Безумовно, суть вектора g в кожній з цих задач різна, так як в 
першому випадку 8 залежить тільки від властивостей біосередовища і 
прикладеного магнітного поля Н, ав другому — Є залежить і від хвилевого 


вектора поля К. Тому можна записати: 


піз = Hy E & = п + (fH + f)n,, (57) 


де покладено, що H = Hn. Цим двом значенням n? відповідають слідуючі 
відношення обох компонент вектора Е (або D): 


E, - HE, (58) 
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тобто кругова ліва і права поляризації хвиль. Рівність абсолютних значень Е, і 
E, (D, i Dj) при зсуві фаз між ними 7/2 означає кругову поляризацію хвилі з 
напрямком обертання вектора Е (Р) відповідно по і проти часової стрілки, 
якщо дивитись вздовж напрямку хвилевого вектора К (або, як принято гово- 
рити, відповідно ліво- і правополяризовані хвилі). Природньо, що при ураху- 
ванні анізотропії біосередовища хвилі стануть еліптично поляризовані. Вико- 
ристовуючи подання вказанних ефектів через вектор С можна записати: 


Ts 
ТИТЕ. зн з) (59) 


де використане Te ж наближення (другого порядку малост! по п (показнику 
заломлення)) iG = F? (для оптичної активності біосередовища), а H = Hn. 
Поблизу лінії поглинання в гіротропному біосередовищі діелектричну 


проникненість є? (без урахування просторової дисперсії і гіротропності 
середовища) представимо у вигляді полюсного виразу: 


Бетте (60) 


Тут ми не будемо вимагати спеціальної малост! коефіцієнта А і відповід- 
но цьому He пишемо постійного члена. Для зв'язку між E i D сд користува- 
тися формулою, що враховує зв'язок типу (52) i (55). В ізотропному біосере- 
довищі можна записати: 


B=, B+ F- 27H [5х] (61) 


де вважається, що: С «(в 2t a Hain, Н = Hn. Такий запис вектора 
Е 


оптичної активності біосередовища дозволяє врахувати одночасно не тільки 
природню оптичну активність біосередовища, але й магнітооптичні ефекти в 
ньому. Відмітимо, що поблизу лінії поглинання компоненти тензора ez (3B0- 
ротнього тензору =,), проходять через нуль і немає підстав для порушення 
збігу його розкладу по хвилевому вектору. 

Дисперсійне рівняння має вигляд: 


1 ~ 
Ра ge (62) 


@ 


і 


ne nl =e), Підставляючи сюди (0) is (60), отримаємо рівняння: 


ЕСЕ UT 
peje 
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Поблизу лінії поглинання (00,) можна вважати: 


f = 2; 
Е > [ctm (64) 


Один із цих коренів рівняння (63) n? існує як при O< Op, TAK i IPH O> Op 
де без урахування просторової дисперсії істотних значень п. Два інших 
існують тільки при о< о, — лівіше точки m, в якій: 


2 2 
і -34 5). п? - (5) (65) 


Співвідношення (60) справедливі для повздовжніх і поперечних елек- 
тромагнітних хвиль при урахуванні просторової дисперсії біосередовища, ї 
гіротропності, для діелектричної проникності є, €, але відносяться тільки 
до частот, досить далеких від центру лінії поглинання: o-o] >> ү, He y — 


ширина лінії. При о-о, < y необхідно враховувати поглинання, тобто 
уявну частину діелектричної проникності, що суттєво може поміняти картину. 


Висновки 


1. Дослідження ізотропного біосередовища, що має властивості гіро- 
тропності за рахунок наявності магнітооптичних ефектів і природної 
оптичної активності, обіцяє великі інформаційно-діагностичні 
переваги, особливо при дослідженнях поблизу лінії поглинання (як в 
методі МРТ проф. Сітько). 

2. Наявність зазначених ефектів, що визивають гіротропію біосередо- 
вища, як в поперечних так і в повздовжніх електромагнітних хвилях, 
призводить до генерації в біосередовищі хвиль, що розповсюджу- 
ються з суттєво різними швидкостями. При цьому швидкості розпов- 
сюдження можуть бути набагато менші швидкості с (в вакуумі) і фа- 
зової швидкості с/п, в біосередовищі. Ця обставина дозволяє стверд- 
жувати, що такі "повільні" електромагнітні хвилі ефективно 
взаємодіють з тканинами біосередовища ( Е-з0, поблизу точки т 
n; , >>1). 

3. "Повільні" електромагнітні хвилі — правополяризовані 1 їх більше, 
так як існують два типа “повільних” хвиль з п2 i пі, в той час як 
“швидких” єлектромагнітних хвиль всього одна з MŽ, і вона лівополя- 
ризована. Таким чином, ця обставина вказує на те, що з живою мате- 
рією (біооб'єктом) ефективніше взаємодіють правополяризовані 
хвилі. Окрім того, природа потребує різновидності, яке більше при 
правій поляризації хвиль, падаючих на біотіло і поглинених ним. 
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ПАРАДИГМА ПОВЗДОВЖНІХ ЕЛЕКТРОМАГНІТНИХ ХВИЛЬ В КОНЦЕПТУАЛЬНИХ 
ОСНОВАХ ФІЗИКИ ЖИВОГО: КРИТЕРІЇ ТА УМОВИ ФОРМОУТВОРЕННЯ, 
ІНФОРМАЦІЙНО-ДІАГНОСТИЧНА ЗНАЧИМІСТЬ 


Ю.В. ЧОВНЮК, Т.М. ОВСЯННІКОВА, Б.Ф. РУДЬКО 


Метою цієї роботи є обгрунтування парадигми повздовжніх електромагнітних хвиль у коцепції 
фізики живого, проведення оцінок і встановлення критеріїв формоутворення подібних хвиль у біооб'єк- 
тах, ix “фізичності”, визначення інформаційної і діагностичної значимості таких електромагнітних 
хвиль для аналізу основних характеристик живої матерії. 


ПАРАДИГМА ПРОДОЛЬНЫХ ЕЛЕКТРОМАГНИТНЫХ ВОЛН В КОНЦЕПТУАЛЬНЫХ 
ОСНОВАХ ФИЗИКИ ЖИВОГО: КРИТЕРИИ И УСЛОВИЯ ФОРМООБРАЗОВАНИЯ, 
ИНФОРМАЦИОННО-ДИАГНОСТИЧЕСКАЯ ЗНАЧИМОСТЬ 


Ю.В. ЧОВНЮК, Т.Н. ОВСЯННИКОВА, Б.Ф. РУДЬКО 


Целью данной работы является обоснование парадигмы продольных электромагнитных волн в 
концепции физики живого, проведение оценок и установление критериев формообразования подоб- 
ных волн в биообъектах, их "физичности", установление информационной и диагностической значимо- 
сти таких электромагнитных волн для анализа основных характернистик живой материи. 
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ОНОВЛЕНИЙ ПІДХІД ДО ТЕОРІЇ ФІЗИЧНИХ 
ПРОЦЕСІВ В СЕРЕДОВИЩАХ З ПАМ'ЯТТЮ 
Р.В. ЛУЦИК, Ю.Л. МЕНТКОВСЬКИЙ 


Київський державний університет технологій та дизайну 


AN IMPROVED APPROACH ТО THE PHYSICAL 
PROCESSES THEORY IN MEDIA WITH MEMORY 
R.V. LOOTSYC, Yu.L. MENTKOVSKY 


Abstract. Using characteristic examples with evident ways of generalisations a com- 
parative analysis of the solutions of the evolutional boundary problems (EBP) and the 
hereditary determinative relations (HDR) of the material media physics is performed, 
on the base of the unical (for EBP and HDR) causality principle the 
Bolzman-Picar-Volterra theoretical approach to the systems with memory (heredi- 
tary systems) is improved and generalised. The actuality of such systems has been es- 
sentially increased during a last few years [2,3]. A possible way of cardinal 
simplification of the hereditary integral-differential equations is demonstrated, the ac- 
curacy of which may be principally controlled and recent intuitively discovered mod- 
els are verified. This paper was divided into two parts. 


ЧАСТИНА 2 


53. Формалізація фундаментального фізичного принципу 
причинності 

Приклади (4), (4a) та (5), (5a), а також інші наводять на думку, що B од- 
нопараметричних задачах із зосередженими параметрами закон розгортан- 
ня процесу у часі під дією внутрішніх і зовнішніх чинників можна шукати 
серед функціоналів типу: 


y(t) = fo(t3(:))&: ат) 
де 


(хі є) (18) 


— це певна узагальнена функція змінних /, т, причому, /(т) репрезентує вза- 
ємодію системи з оточенням (типу модифікованого джерела I (9 (16) у випад- 


ку ЕКЗ виду (3), (5), або "сигналу" x(t) увипадку CBC (4)). 

Формалізувати фундаментальний фізичний принцип причинності мож- 
на, наприклад, наступним чином. 

-- Оскільки наслідок не може випереджати причину, то варіювання 
факторів, що визначають, наприклад, взаємодію системи з оточенням, 0/(т) 
не повинно впливати Ha y(f) (17), коли 


tot (19) 
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Отже, повинно бути 


&y(t )- Т0) (так «0, якщо > f. (20) 


—® 


В лінійних по взаємодії з оточенням теоріях інтегранда Q(t,«| (т) no- 
винна мати вигляд 


Q(t x| i()) E К(г, т) С); (21) 
Значить, в таких теоріях 


чо = R(^*), | (22) 


а варіація в (20) набуває вигляду 


dy(t) = ©) (г) т =0, аб 
ЯКЩО T >É. 


Оскільки ж варіація 5} (т) — є довільною, TO з (23) випливає, що 
R(tx) = 0, Kom « >t. (24) 


Таким чином, в класі узагальнених функцій ядро R(t,x) (24) повинно ма- 
ти наступну загальну структуру: 


R(t, t) = a(t, t)8(t - т) -+ e(t ~ 1) b(t, ws (25) 

Цей буде відбивати в загальному вигляді фундаментальний і унвер- 
сальний фізичний принцип причинності для систем з одним зосередженим 
параметром. 

Доданок з 6-функщею Дірака відповідає миттєвій реакції системи на 
зовнішню дію (що в нерелятивістській фізиці не суперечить принципу при- 

чинності); доданок з 0- "функцією Хевісайда ' ‘відрізає” майбутнє, котре He 
може впливати на поточне (“саморух” тут включено до b(t, т)). 

Нижче показано, що аналогічні, хоч і більш загальні побудови, можна 
виконати і для систем з розподіленими параметрами. Звичайно, описаний 
підхід можна узагальнити також на системи з кількома параметрами стану 
(або характеристиками процесу) [4]. 

Конкретизація функцій a(t,t) та b(t,x) (25) вже не може бути універсаль- 
ною. Вона класифікує співвідношення, вбираючи в себе із запасомі ЕКЗ ти- 
пу (3), i CBC типу (4). Дійсно, покладаючи 


J(t) = (2); а(х) =0; т) =К(ь*); (26) 


(див. (16), (5a) та (14а)), отримуємо розв'язок (5a), (5), (6) ЕКЗ (3). Обравши ж 
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a(t, t) =]; b(t,x) = К (т) : (t) - x(t); (27) 
(див. (4a) та (146)) одержуемо CBC (4). 

Таким чином, показано, що розв'язки (5) ЕКЗ (3) i CBC (4) дійсно є 
окремими випадками єдиного фундаментального принципу причинності 
(25), (23), (21), (17) ірівноправно визначають закон еволюції параметрів y(f) 
під дією внутрішніх і зовнішніх чинників (x(7) у випадку (4) та Г (t) (16) y ви- 
падку (3), (5)). Як ЕКЗ (3), такі СВС (4) належать (неявно і відповідно явно) 
до еридитарних (спадкових) теорій. Тому слід було б, мабуть, скорегувати 
термінологію і називати, скажімо, ЕКЗ (3) ra CBC (4) каузальними (причин- 
ними) крайовими задачами та співвідношеннями відповідно — ККЗ (3) та 
КС (4). Або щось подібне. Повторюємо ж бо, що адекватна термінологія по- 
легшує "роботу" інтуїції. 


54. Стосовно ЕКЗ (1) та СВС (2) для систем із розподіленими 
параметрами 


Ці проблеми обговоримо також на конкретних і виразних прикладах, які 
допускають очевидні (для фізика-теоретика) узагальнення Аби з самого 
початку не ускладнювати теорію врахуванням граничних умов, розглядати- 
мемо еволюційну задачу у безмежному просторі. За характерний і виразний 
приклад оберемо проблему теплопровідності в тілах, в котрих можна знех- 
тувати супутніми реологічними процесами; тобто в яких немає помітних по- 
токів речовини (наприклад, дуже слабке розширення при нагріванні, немає 
достатньо рухомих складових типу рідин тощо). Для таких тіл достатньо 
розглянути одновимірну ідеалізовану задачу: задачу про поширення тепло- 
ти вздовж однорідної нескінченно довгої тонкої дротини (певні вихідні мір- 
кування щодо систем обмежених розмірів буде висловлено наприкінці), 
теплоізольованої з боків (узагальнення на дво- i тривимірні тіла не привно- 
сить суттєвих принципових ускладнень в теорію, лише дещо захаращуючи 
формули). Отже, розглянемо задачу Коші для одновимірного рівняння теп- 
лопровідності, тобто рівняння, доповненого початковою умовою: 


ar , SUE | 
Тат о) (28) 
T(x0) = T, (x). (29) 


Tyr T(x,f) — температура дротини, як функція координати (x) i часу (t); 
a = x/cp — коефіцієнт температуропроводност! (A --- теплопровідність дроти- 
ни; c — il питома теплоємність; р — густина речовини -- усі ці величини BBA- 
жаємо тут сталими); с = 4, /ср — "локальне темпего" "рне джерело” (а, — пи- 
тома потужність внутрішнього (прихованого) тепловиділення за рахунок, 
скажімо, фазових перетворень в дротині). 

Відсилаючи за браком місця по доказ до літературних джерел (див., 
напр., [8. 9]), запишемо розв'язок ЕКЗ (28), (29) у вигляді 


T(x, t) = [n (х) (20; x, x") dx' + [а je: т)К(о т; x, x’). (30) 
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Тут 


б) 


= . чо 1 4a? (r1) 
K(541; x E e 31 
( у oe) |) И 


Перший доданок є розв'язком однорідного рівняння теплопровідності 
(28) (колис = 0) з початковою умовою (29); другий доданок є розв'язком He- 
однорідного рівняння (28) (с = 0) з нульовою початковою умовою (Т, (x) = 0). 
Тут — повна аналогія з розв'язком ЕКЗ (3), оскільки прийнято, що f, = 0. 
Дал! можна буквально повторити усе, що було сказано стосовно EK3 (3). 
А саме, запровадимо функцію Гріна 


G(t, я; X, x") = e(t - t)K(t, т; X, х) (32) 
та нове джерело 
a(x) = 0(c)o(x,t) + B(A). (33) 


Нагадаємо, що тут узято & = 0. 
Тоді (30) можна подати у вигляді 


T(x.f) = [de а(х xT (x80) +8. (т), 
бо (x t) = ө(2)с(х, г) 5 
цілком аналогічному (5a). Значить, i принцип причинності для (34) виглядати- 


ме аналогічно попередньому (17)-(26). 
До речі, можна показати, що функція Г ріна G(t, UE x) задовольняє рівняння 


(34) 


дСІ, vx, x’) Е 
ot 


а 200655) = 8(t-1)8(x-2’), (35) 


аналогічне (15). 

Перш, ніж переходити до порівняльного аналізу розв'язку (30) ЕКЗ 
(28), (29) з відповідними СВС, прокоментуємо фізичний зміст зображення 
(30), оскільки на нашу думку у традиційному трактуванні формул типу (30) 
не все гаразд. Скажімо, в [8] вказано на певну неточність вихідних посилань 
теорії поширення теплоти, оскільки кожна точка-джерело x' впливає на всі 
(тобто як завгодно далекі від неї) точки стержня х миттєво. В такому катего- 
ричному твердженні є певна непослідовність, хоч зовні воно правдоподіб- 
не. Справді, при виведенні рівняння теплопровідності суттєво враховуєть- 
ся закон Фур'є |З, 9]: 


із WT (7,1), 


B якому j — щільність теплового потоку, що визначає, зокрема, швидкість по- 
ширення теплоти, а вона € скінченою: / | < œ. Тому формула (30) є досить по- 


слідовною; вона враховує скоріше той внесок від точок-джерел х’, який на да- 
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ний момент часу / акумульовано вже в точці х. (Поширення теплоти було ж бо 
1 до початкового моменту (t=0), коли УТ = 0.) 

Звичайно, рівняння (28) не може бути універсальним. Воно стосується 
квазістаціонарних, тобто, не дуже швидкоплинних процесів поширення 
теплоти. Це випливає з аналізу його виведення (див., напр., [9, 31). Але для 
таких воно є, нанашу думку, послідовним. Отже, розв'язок ЕКЗ має спадко- 
вий (еридитарний) і в одночас нелокальний характер. 

За аналогією з теорією систем із зосередженими параметрами можна 
очікувати, що і в загальному випадку між ЕКЗ і СВС існує паралель: СВС 
для систем із розподіленими параметрами мають бути, взагалі кажучи, ери- 
дитарними і нелокальними. 

На прикладі безмежного в'язкопружного прямого стержня побудуємо 
тепер відповідне CBC для порівняння з розв'язком (34) ЕКЗ (28), (29). 

В теорії невеликих однорідних в'язкопружних стержнів, як відомо, з ус- 
піхом використовують СВС типу Больцмана-Вольтерри: 


e(t) = 20 + [Ks t) °С) dt (36) 


Tyre(t) = em — відносний розтяг (l(t), 1, — відповідно поточна ino- 
0 
чаткова довжини стержня), с(1) = F(r)/S — напруженість y стержні (F(t) — 
сила, що діє на стержень вздовж його осі, 5 — площа поперечного перерізу 
стержня, яка за малих розтягів є практично сталою), K(51) — ядро 
повзучості. 

Як бачимо, (36) має форму (4), коли покласти y(t) 2 £(t);x() = с(д/Е.От- 
xe, (36) e типовим представником CBC (4). Для конкретності подальші по- 
будови прив’яжемо саме до (36). 

Розглянемо ідеалізовану задачу малих локальних розтягів безмежно до- 
вгого в'язкопружного стержня. Нехай в результаті створених локальних на- 
пружень o(x, 1) відбулися локальні розтяги; причому, довільна точка Р стер- 


жня з координатою х перемістилася в точку x' на малу відстань u(x): 
Р{х} з pdy +u(x, г)). (37) 


Нехай близька до Р точка О 3 координатою x, перемістилася Ha малу від- 
стань u(x, ): 


O{x,} > оҷ, +и(х,, г)). (38) 
Оскільки вважаємо, що 
Ax =x, -x (39) 


— це мала величина, TO, очевидно 
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ди 
Au = u(x,t) - u(x, г) s — Ax. (40) 
Ox 
Розтягом e(x,f) в даному випадку буде 


ae (41) 


Постулюемо, що для в’язкопружного стержня можна записати ериди- 
тарне співвідношення типу (36), але через Te, що розглядається динамічний 
процес, це співвідношення має бути ще й нелокальним (динаміка напру- 
женьзалежить від руху всіх ділянок стержня і навпаки). Тому слід покласти 


e(x; t) = ця + [dx jx k (tn: X, х) Be) (42) 


EG) і = E(x’) 


Твердження (42) є, звичайно, вірогідною гіпотезою, яка потребує прак- 
тичної перевірки на якомусь виразному експерименті. Тому тут все ще 
попереду. 

Спорідненість виразів (30) і (42) очевидна. Звідси неважко за аналогією 
з (4a) записати функціонал (42) у вигляді 


(x,t) = je асо ~1)8(x ох) +I (1,75 xx’) -9(2)) d (43) 


Позначення за аналогією з (4a) зрозумілі. Функціонал (43) є спорідненим з 
(34). Далі можна конструктивно сформулювати (як і у випадку систем із зосе- 
редженими параметрами) єдиний для (34) i (43) принцип причинності, з якого 
ці співвідношення випливатимуть як окремі випадки. Це робиться просто, то- 
му виклад опускаємо. 

Ясна річ, переходячи від ідеалізованих нескінченних об'єктів до реаль- 
них (скінчених), треба буде, як правило, враховувати не лише початкову, 
але й граничні умови. Це з нашої точки зору стосується не тільки ЕКЗ, але й 
СВС, оскільки, як бачимо, вони споріднені. Таким чином, ідеї, що сходять 
ще до робіт Е. Пікара та В. Вольтерри, слід доопрацювати та узагальнити в 
дусі викладеного вище. Тим більше, що СВС за певних умов повинні пере- 
ходити у розв'язки певних ЕКЗ. Не слід також протиставляти ЕКЗ типу (1), 
(3) та (28), (29) i CBC tany (2), (4), (43) як нібито не еридитарні та відповідно 
еридитарні. Всі вони еридитарні і, взагалі кажучи, не локальні. Тому загаль- 
ною математичною формою фундаментальних фізичних законів є, очевид- 
но, сукупність інтегродиференціальних рівнянь і співвідношень, що задо- 
вольняють принципу причинності. 

Являє неабиякий науковий інтерес дослідити умови переходу від ериди- 
тарних співвідношень до миттєвих і формально-миттєвих та від нелокаль- 
них закономірностей до локальних, як кардинально простіших і зручних 
для практичних застосувань. На прикладі систем із зосередженими пара- 
метрами коротко обговоримо і дану проблему. 


27 


Р.В. ЛУЦИК, Ю.Л. МЕНТКОВСЬКИЙ 


55. Відносно умов наближеного переходу від інтегродиференціальних 
до чисто диференціальних рівнянь 


В роботі [9] було показано, що вимушені коливання точкового тіла маси 
т, зв'язаного з нерухомою основою в'язкопружним стержнем описуються 
наступним інтегродиференціальним рівнянням 


¢ = em EU - ILC «9| TX (44) 


Тут x(t) — координата матеріальної точки (X — ii прискорення); E — по- 
чатковий модуль Юнга; 5 — площа поперечного перерізу стержня; і — ïo- 


го рівноважна довжина; R(t,t) — так зване ядро релаксації (напружень); 
f(t) — примусова сила. 

(В формі Гамільтона мали б систему з двох рівнянь першого порядку чи- 
сто еволюційного типу. Рівняння Ньютона (44) еквівалентне їй. Це приклад 
згаданого в 51 еквівалента еволюційним рівнянням.) 

В принципі, можна побудувати розв'язок інтегродиференціального рів- 
няння (44). Але цей розв'язок має занадто складну аналітичну структуру, 
яку важко якісно проаналізувати. Аналітичні труднощі стрімко зростають, 
якщо мова йде про цілу систему тіл з в'язкопружними зв'язками. Між тим, 
така проблема є технічно вельми актуальною в корабле- та літакобудуванні, 
та й скрізь, де важливою є проблема вібрацій. Для практичних цілей слід на- 
ближено спростити рівняння (44). Емпірично спрощені варіанти, відомі ін- 
женерним наукам. Та бажано мати грунтовний теоретичний підхід до проб- 
леми, який давав би змогу оцінити фізичну суть апроксимацій та визначити 
їхню точність. 

У роботах [5, 9) було запропоновано спеціальні функціональні розклади 
функціоналів типу Фредгольма і Вольтерри (до якихналежить і функціонал 
у фігурних дужках (44)). Було показано, що вираз у фігурних дужках (44) 
можна подати у вигляді наступного збіжного ряду: 


EJL PRO -2C)'o..()-x9(). (45) 
Tyr 
Q(t) =K, (£t); кут) = [К„.(6 Баку; 
Кот) = (т); x(n) - ZE. 


Обмежуючись трьома доданками збіжного розкладу (45), отримаємо замість 
(44) чисто диференціальне рівняння: 


KOO) = KORO (04) +70), «m 


в якому запроваджено скориговану змінну (ефективну) Macy 


(46) 
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0) = m+? о (48) 
змінний коефіцієнт пружності | 
Ко) = n (1+0, (#)); (49) 
змінний коефіцієнт в'язкого тертя 
ау т 0, (г). (50) 


Для максимального спрощення попереднього якісного аналізу рівняння (46) 
та побудови його розв'язків вдамося до популярної моделі, яка забезпечує 
прийнятну для інженерних цілей відповідність експерименту. 

А саме, покладемо 


К(1,т) = пе Коб), (51) 
Параметри п та В разом із модулем Юнга Е паспортизують стержень 
стосовно досліджуваного процесу вимушених коливань під дією примусо- 
вої сили f(t) (47) (п — коефіцієнт повзучості; т, = ИВ — час в'язкопружної 
релаксації). За такої (51) моделі неважко отримати. 


0100) -3(e*-) e)-3:(-e7)-e" | 


| 
B | 


пм 1 TU CE 2 | (32) 
RO 1 (1-е) a a 
Усі функції О, (г) мають релаксаційний характер i для 


Í 22 То = (53) 
B 


перетворюються практично Ha константи 


Q =-1; 0, =; Qa- 
В p? p? (54) 


Значить, згідно з (48)-(50) перетворюються Ha константи також кое- 
фіцієнти рівняння (47): 


(55) 
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Тобто при умові (53) рівняння (47) можна замінити на практично тотожне йо- 
му рівняння: 


X +28% + @7х = (9, (56) 
Но 


в якому 


TEN LUE | (57) 


Рівняння (56) є добре вивченим рівнянням вимушених коливань, набагато 
простішим і прозорішим за вихідне рівняння (44); причому, в ході спрощення 
одержано очевидну можливість оцінки точності усіх розрахунків (на якій в си- 
лу її очевидності не спинятимемося). 

Відомо, що наявність "в'язкого тертя", представленого в (56) доданком 
25х, веде до дисипації механічної енергії, яка зумовлена тут в'язкістю стерж- 
ня (див. у, за умови п 0 (55)). Отже, роль пружності (k, #0) і в'язкості 
(у, = 0) одночасно перебирав на себе в'язкопружний стержень, що виразно 
і просто випливає з рівняння (56). 

Як добре і широко відомо, за певного співвідношення між параметрами 


5 та о, рівняння (56), його розв'язок після певного часу релаксації (т, = 1/8) 
виходить на усталений режим, який можна представити виразом 


x(t) я dl. [re sin a(t —t)dr. (58) 
Ноо 5 
Тут 
o = 04-82. (59) 


(В такому вигляді вираз (58) зберігає ще мало істотні “уламки” перехідного 
процесу, суттєві лише для проміжку 0 < t < 1/8.) 

Зазначимо принагідне, що праву частину (58) теж можна наближено 06- 
числити за допомогою функціональних розкладів типу (45); що може слу- 
жити альтернативою традиційного зображення примусової сили f(t) у BH- 
гляді розкладів Фур'є. 

З наведеного вище випливає, що нелокальні теорії з пам'яттю є широко 
правомірними i що не слід занадто їх лякатися і уникати. Слід лише присту- 
пити до їхнього глибокого і систематичного вивчення та до розробки набли- 
жених методів спрощення. Одна справа емпірично спрощені рівняння, які 
так чи інакше доводиться підбирати; інша -- систематичні методи спро- 
щення, здатні чітко виявити фізичний сенс і точність наближень. Мораль, 
гадаємо, зрозуміла. 

Цими, як нам здається, прозорими міркуваннями тут і обмежимося. 
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На виразних прикладах, що допускають широкі узагальнення, виконано порівняльний аналіз 
розв'язків еволюційних крайових задач (ЕКЗ) та спадкових визначальних співвідношень (СВС) фізики 
матеріальних середовищ. На основі об'єднувального для ЕКЗ i CBC принципу причинності уточнюєть- 
ся та узагальнюється висхідний ще до Л. Больцмана, Е. Пікара та В. Вольтерри підхід до теоретичного 
описування систем з пам'яттю (ередитарні системи), актуальність яких в останні десятиліття істотно 
зросла [2,3]. Показано можливі шляхи кардинального спрощення ередитарних інтегро-диференціаль- 
них рівнянь, що допускають оцінку точності отриманих наближень та обгрунтовують інтуїтивно “нама- 
цані" прикладниками моделі. За технічними причинами роботу розділено на дві частини. 


ОБНОВЛЕННЫЙ ПОДХОД К ТЕОРИИ ФИЗИЧЕСКИХ ПРОЦЕССОВ В МАТЕРИАЛЬНЫХ 
СРЕДАХ С ПАМЯТЬЮ 


РВ. ЛУЦЫК, ЮЛ. МЕНТКОВСКИЙ 


Навыразительных примерах, допускающих широкие обобщения, проведен сравнительный ана- 
лиз решений эволюционных краевых задач (ЭКЗ) и наследственных определяющих соотношений 
(НОС) физики материальных сред. На основе объединительного для ЭКЗ и НОС принципа причинно- 
сти уточняется, восходящий к Л. Больцману, Э. Пикару и B. Вольтерре, подход к теоретическому описа- 
нию систем с памятью (эредитарные системы), актуальность которых в последние десятилетия 
существенно возросла [2,3]. Показаны возможные пути кардинального упрощения эредитарных интег- 
ро-дифференциальных уравнений, допускающие оценку точности получаемых приближений и обосно- 
вывающие интуитивно “нащупанные” прикладниками модели. По техническим причинам работа 
разбита на две части. 
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Abstract. The modern theory of cell-to-cell transmission is commonly based on ideas 
about chemical intermediaries (transmitter agents) securing transmission between 
two neurons in the synaptic cleft or between the motor neuron and muscle fiber in the 
neuro-muscular junction. About ten millions of transmitter molecules are releasing si- 
multaneously under the influence of one nerve impulse and then the chemical reac- 
tions are taking place. In order to describe the major events in synapse the simplified 
scheme is used. Kinetic equations corresponding to the scheme of biochemical reac- 
tion in synaptic cleft have non-linear form. Using some special approximations it is 
possible to investigate the main characteristics of such a system. Considering the syn- 
aptic transmission as the process which is isomorphic to the critical phenomena in bi- 
nary liquid mixtures near the critical mixing point it is possible to calculate the 
correlation function of the density (concentration) fluctuations and determine the 
size of the activation zone and other critical parameters. 

Keywords: critical phenomena, correlation function, correlation length, critical pa- 
rameters, synaptic cleft, synaptic transmission. 


Introduction 


Synaptic transmission is one of the most important problems of modern bio- 
physics. Such a great interest is caused by the role of synaptic transmission on dif- 
ferent levels of biological objects’ functionality. So it is natural to consider the or- 
dinary cholinergic synapse. In agreement with modern imaginations about the 
nerve impulse transmission the main consequence of events in such a cleft may be 
described as follows (see, for example, [1]): Acetylcholine (ACh) is synthesized 
by cholineacetyltransferase and is stored in spheroid vesicles, then releases and 
reacts with specific acetylcholine receptors (R). The formation of the transmit- 
ter-receptor complex (AChR) produces conformal changes in the postsynaptic 
membrane and therefore the changes in the membrane potential. F inally ACh is 
either inactivated by acetylcholinesterase (AChE) or is removed by diffusion. 

The scheme of biochemical reaction may be presented in such a way: 


5 у 
ACh + К AChR —-—» PR —»—» В (1) 
БАЕ 
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Here PR is the product of ACh destruction by the AChE ands, аге the coeffi- 
cients of reactions velocities. 

The process of ACh release is a cooperative one: about 107 ACh molecules 
are releasing simultaneously under the influence of one nerve impulse [1]. 

In order to analyse synaptic transmission let us simplify the base scheme of 
biochemical reaction as follows: 
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ACh + R—&— AChR —5—, R (2) 


Kinetic equations describing the temporal evolution of concentrations X 
and Y of R and ACh have the following form: 


dx | 3 
агт е9) m 
d 

P. = f(t) — ky xy (4) 


Here f (t) is the source function that describes the intensity of ACh releasing. 
In general case the analyzing of the system (3)-(4) is difficult because of its 
non-linearity. So let us use some special approximations. 


Stationary Source Function Approximation 


If f = const then one has from equations (3)-(4): 


ах _ (5) 
m k,(1 - x) - xy 
dy 
Z -f-k (6) 
di f-kxy 
so the stationary points are: 
x, =1-f/k, (7) 
and 
-ka fk 
зав ©) 


Using new parameters 4 =k, /k, and o = f/k, one can transform the formulae 
(7)-(8) to the following form: 


x, =1-Ф (9) 
апа 
у, -Xo/(1 - o). | (10) 


Figure 1 shows the dependence of stationary points on parameters À and o. It 
is important to stress that x, does not depend on. It could be seen from (10) that 
stationary value of the ACh concentration (i.e. y, ) increases to infinity when pa- 
rameter ф approaches to 1. Physical sense has the situation when y, changes in 
range from 0 to y, — some critical value. So from inequality 


У; S Y orit (11) 


one has 
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Fig. 1. Dependence of the stationary points on parameters A and q. 


у crit 
< И, ар. 
J Ed ua UE Ey E (12) 


The value ої 7 ах may be treated as the characteristic of the informational capability 
of the transmission channel. 
Slow Oscillation Approximation 


Iff =f ( f) * const — in this case for the Fourier components from the system 
(3)-(4) one has: 


| -іох(о) = k,8(@) - k,x(o) - k,x(@)* у(о) (15) 
“fay (o) = f (o) - kix(o)» у (в) | (16) 


Let us use in equations (15)-(16) the Taylor approximation for the x(o) and 
у(о) respectively. The system (15)-(16) may be simplified to the next one: 


—юх(о) = k,8(@) - k,x(@) - kx(o)y, - k,x"(o)y, (17) 
Hey (0) = /(2) - хьз (0) Вх," (о) (18) 


were x, = KOLO ys | y(@)da, х, -5 fx(o)o%da y; - | y(o)o*do. From 
the last system one has: 
io -ks 
iy 19054, 19 


Á k, * ky, 
/ (9) (20) 


d EE - бо е 
у (с) "I kx, У(о) kx, 


7 k (19) 


From the system (19)-(20) one can get: 
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k, 4 К, +К У 1 P 
xlo) = —— —— |, 21 
(о) ky; ky; xm 20 
ңа En І | 
= x, kx, (22) 
»(e) Kos ? 
where 
3 
eof etr) 
CGB) =—— e (23) 
Е | | 
ЗВ В 
and 


R(t, o, p) = [Rej (24) 


k k, +k 1 
One en e» 
iE Хо, “| | 
ie mm m 
If n << |f(1)|then R(t,a,B) = f(t)G(t; o. B)so the source function f(t)may be 


interpreted as the modulating one. 


Correlation Functions And Activation Zone 


According to [1] synaptic transmission may be considered as the process 
which is isomorphic to the critical phenomena in binary liquid mixtures near the 
critical mixing point. Using some speculative approximations (namely 
plane-parallel geometry of synaptic clefts, location of receptors оп the surface in- 
side the layer, etc.) one may apply the results of investigation of binary liquids 
near the critical point for estimation of the linear size of the activation zone for 
AChR. 
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The main results ofthe theory of electromagnetic critical scattering and other 
critical phenomena in the case of liquid systems are based on the expression for 
the pair correlation function G, (r) of density (concentration) fluctuations [2-4]. 
The Munster method [3] could be used to determine the correlation characteris- 
tics of finite-size systems with geometry of a plane-parallel layer (-h<z <h, 
-оо < x, y « +0). The basic assumption of this method is concluded in solving the 
integral Ornstein-Zernike (OZ) equation with the delta-function as the zero ap- 
proximation for the direct correlation function f(r) and the Green function of the 
Helmholtz differential operator as the pair correlation function G, (r) under the 
zero boundary conditions. One can get the first iteration for the direct correlation 


function f Or) from the integral OZ equation using the zero approximation for 
the pair correlation function c? (r). Then the first iteration for the pair correla- 


tion function С Q can be determined by solving of the corresponding differen- 
tial OZ equation. Calculation of such characteristics as correlation length (acti- 
vation zone), shifts of the critical temperature and the density (concentration) 


can be studied using the expression for the pair correlation function со (г) in the 


first approximation. It is necessary to point out that the zero approximation for 
the pair correlation function was investigated in [1,4]. 


1) The direct and the pair correlation functions 
Integral OZ equation 


Gr) = f(r) +(P)[F (7) G(r - )ar (27) 


generates the corresponding differential OZ equation in the double-moment ap- 
proximation: 


(v? = к) -G,(r) = 0, (28) 


where: C, = (р) | f(r)dr is the zero moment of the direct correlation function, 


C= n | (r)r*dr is the second moment of the direct correlation function f(r), 
к? ES 1-C, ; 
С, 


jo (r) = a -8(r).For a plane-parallel layer of thickness d = 2hthe initial iteration 


(p) — is the average density. In the zero approximation 


for the pair correlation function [1,4] 


23) 
7 з 2(2т+1 2 
69 (y, z) Sui да" тату ner id n 
\ 


4h? 2h 


36 


Synaptic Transmission as the Cooperative Process 


wherep = үх? +у and К, (и) ів the cylindrical Macdonald function. The first iter- 
ation for the direct correlation function can be expressed in the following form: 


г (ов) = ÈA (9): co eres (30) 


Then one can obtain from equation (28) using formula (30) for f (r) the following 
expression: 


І 
" : 2m +1) т Й 
"9 збу e PANG te ня = |^ | D (31) 


Differential OZ equation (29) gives for the first iteration of G(r): 


| i5 E cti 
| «(т + 2) М + j^ 
Of: oy NS - 
6; (p.z) = 2hn ONE mo pe cm h? чо h 
; (32) 
fa ota || (eral 
wae E ZT 
K, || 2 Па e А WS EIS 
jee ES T^ Me Е 
| 
\ 


where € = Jx? +С, /С,. 


The main contribution to G, (р, 2 ) is given by terms with small m in (32). Therefore, 
one has 


feas [re] n e] LB) е 


Whenr — 0 in the zero approximation С, (r — o. The first iteration for the pair 
correlation function has a finite value even whenr = 0. So in the first approximation 
the pair correlation function has more realistic behavior then in the zero one. 


2) The correlation length of order parameter fluctuations 


It is natural to determine the correlation length R, of order parameter fluctu- 
ations according to the ра relation [1,4]: 


EC т? «dr 
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Fig. 2. Dependence Кс оп т. 


After integrating (34) one can determine such an expression for the correlation 
length Кс in the case of a plane-parallel layer: 


RO =h $ + 4 5 (1-4), 
QT 


ЕЕЕ ney? + Сь — 
4 са 


(35) 


According to [1,2] the critical behavior of the parameter к may be presented as fol- 
1 


lows: к - т" fo (л/т?) = si: ж) where у = 063, с = 192 and В «033 — is 


the critical exponents,t = (T -T )/T, and Ax = (x - хо )/x_ —is the deviations of 
the temperature and the concentration (density) from their critical values. The as- 
ymptotic behavior of the functions f, p and } (у) are the next: f, (y> о) ~ у, 

f(y > о) у" and h(»290)-5(yo 0) =K- Here x, is taken to be 
Ко = 5nm^ , JC, /С, = 1nmis the Debye length of relaxation, thickness й = 10 nm. 


Figure 2 demonstrates the dependence R, ont. 
3) The shifts of the critical parameters 


In the close vicinity of the critical isochore one has the following formula for 
the critical temperature: 


Тс 


аа Еа 
1 +(2/KA,)¥ 


(36) 


where Тс is the new critical temperature, T, is the bulk critical temperature is taken 


to be 300 К, A, = /1+ $ + 16 
п? п? +4h’C, /C, 


= 136. In the vicinity of the critical iso- 


1 
therm (т << Ax^) one obtains the following formula for the critical density x¿ (K), 


where x, is taken to be 300kg /m°: 
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40 r^ ЕТ = 


a (37) 


Formulae (36),(37) demonstrate the dependence of shifts of the critical temperature 
AT, =T} -Тс and the critical density Ах, =x¢ 7% on geometrical factor K: 
1 


АТ, ~ КУ for the critical temperature and Ax, ~ K 7-1 for the density where the geo- 


metric factor К = 2h/R, and the critical exponent ô = 4.5. Figures 3 and 4 demon- 
strate the dependence ДТ, and Ax, on К. 

It is important to stress that this results are in agreement with the experi- 
ments carried out by Lutz and co-authors and other investigations [5,6]. As it is 
seen, with increasing of the thickness d = 2hof a plane-parallel layer critical tem- 
perature and density are approaching the bulk values Т. and x. 


The Restricted Cylinder System 


In order to investigate the synaptic transmission it is naturally to consider the 
system with the geometry of the restricted cylinder: 0 <p < and -h < z < h. 
Similar to the previous case the boundary conditions are zero to satisfy the 


| K 


| 

| 

| 
0 zd | | RA 
10 15 20 25 


Fig. 4. Dependence Ax, on K. 
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Fig. 5. The pair correlation function. The dependence on р 


limiting transition to the spatially infinite system. Then the pair correlation func- 
tion may be found in such form: 


G, (р, 2) E Y > бань г) а (38) 


пед m=0 


where J, (x — is the Bessel function and и’ — is the solution of the equation 
J, (x) =0. Taking into account that 


[s Js n(n +05)z 
"а h (39) 


Tah momi IEEE 7 
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Fig. 6. The pair correlation function. The dependence on z 
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Fig. 7. The dependence of the correlation length on the temperature 


MIN 
"O УЛ) , плн oa (40) 
к” + + 
4h? а? 


where V = лай — is the volume of the synaptic cleft. Thus for the pair correlation 
function one has the expression (J, (u) — is the Bessel function): 


i БЕ й 24 


1 оо оо 1 
G (pz) ==: 2 | ic 
‚(е ) V 22 т Hm) : z^(2n +1). u? ( ) 
к“ + PE 


4h? а? 


Figures 5 and 6 demonstrate the behavior of the pair correlation function for 
the system with the parameters a=100nm and h=10nm in the vicinity of the 
bulk critical point. 

The correlation length may be found from the formula (34) and after integrat- 
ing one has the next expression: 


2a) 


sy а Te +h? eS m 
2209 . | о is д? (2п +1) 


L 


(42) 


3 С 


$ (22 *1)u,Ji (ни Ки 


n^ (2n«1) Vs 

4h? a 

The temperature dependence of the correlation length is shown in figure 7 
(for the restricted cylinder system к, is taken to be x, = 10" и”). 


where кі, = к? + 
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Table 1. The sifts of the critical parameters 


[x] x | AT, | Axo |K] K | Ат, Ax, 
| 125 | 97888 | 132.141 | | 25 | 18.771 | 87131 
|| 50 | 96976 | 131803 | s | 15.796 | 83.444 
5 | 75 | 96806 131140 | 20 | 75 | 15.222 | 82670 
| | 100 | 96746 131.718} | 100 | 15020 | 82392 
| 125 | 96719 131708| 125 | 14925 | 82261 
25 | 43258 | 107150 | | 25 | 14925 | 82.261 
| 50 | 41308 105944 | 50 | 11645 | 77247 
10 | 75 | 40946 | 105713 | 25 | 15 | 11004 | 76.137 
| 100 | 40836 105.631, 100 | 10776 | 75731 | 
125 | 40.777 | 105.593 | 125 10670 | 75540 | 
| 25 | 26163 | 94626 | 
© | 50 23.609 | 92.248 | 
5 75 2312 | 91773 | 
| 100 | 22951 | 91603 | 
| 125 | 22.871 | 91.524 | 


The shift of the critical temperature in the bulk critical point is 


AT. T? 
aey азу 
K? Kj 
and the sift of the critical concentration 
x 
Axe = 8 > 
A я H и (44) 
(ук Ki 


The shifts of the critical parameters depend on the geometrical factors K, = 2hx, 
and К, = ак,. The Table 1 demonstrates the shifts of the critical parameters for the 
Тс =300K and x, =300kg/m’. It is seen that these shifts may be significant. 
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Conclusions 


1) As it’s shown above the values of stationary points of the concentration of 
ACh and R and the informational capability of the transmission channel depend 
on the kinetic coefficients. Therefore changes of temperature and influence of 
different external fields affect the kinetic coefficients and as the result the station- 
ary points and the informational capability of the channel. 

2) The source function in the case of its slow temporal changes modulates the 
temporal evolution of the AChR concentration. 

3) The correlative characteristics of the system mediator-receptor may be 
described with the help of the first approximation of the pair correlation function 
that has more realistic behavior then the zero one. 

4) The correlation length determines the zone of reagents activation. It de- 
pends on the geometric factor of the cleft and on the temperature as well. So 
conformal changes of the geometrical form of the synaptic cleft and shifts of the 
temperature affects on the activation zone. 

5) As the result of the space limitations shifts of the critical temperature and 
the density take place. Their dependence on the geometric factor must be taken 
in considering when studying the critical phenomena in such a system. 


СИНАПТИЧНИЙ ПЕРЕХІД ЯК КООПЕРАТИВНИЙ ПРОЦЕС 
О.В. ЧАЛИЙ, О.М. ВАСИЛЬЄВ 


Сучасна теорія міжмембранних взаємодій базується на уявленнях про хімічні реагенти що за- 
безпечують зв'язок між двома нейронами в синаптичній щілині чи між нейроном та м'язовим волокном. 
Близько десяти мільйонів молекул-реагентів синхронно вивільнюється під впливом нервового імпульсу 
в результаті чого має місце хімічна реакція. Для опису процесів, що відбуваються B синапсі, використано 
спрощену схему, причому відповідні кінетичні рівняння мають нелінійну форму. Розглядаючи синапти- 
чний перехід як процес, що ізоморфний критичній поведінці бінарних сумішей, можна розрахувати ко- 
реляційну функцію флуктуацій концентрації (густини) і визначити розміри зони активації та інші 
критичні параметри. 


СИНАПТИЧЕСКИЙ ПЕРЕХОД КАК КООПЕРАТИВНИЙ ПРОЦЕСС 
AB. ЧАЛЫЙ, А.Н. ВАСИЛЬЕВ 


Современная теория межмембранных взаимодействий базируется на представлениях о химиче- 
ских реагентах, которые обеспечивают связь между двумя нейронами в синаптической щели или между 
нейроном и мышечным волокном. Около десяти миллионов молекул-реагентов синхронно высвобож- 
дается под действием нервного импульса в результате чего имеет место химическая реакция. Для описа- 
ния процессов, которые происходят в синапсе, использована упрощенная схема, причем 
соответствующие кинетические уравнения имеют нелинейную форму. Рассматривая синаптический 
переход как процесс, изоморфный критическому поведению бинарных смесей, можно рассчитать кор- 
реляционную функцию флуктуаций концентрации (плотности) и определить размеры зоны активации 
а также другие критические параметры. 
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Abstract. In smooth muscle, Ca”*-dependent ion channels are one of the effector mol- 
ecules involved in physiologically important regulatory processes employing Ca” as a 
second messenger. These channels couple changes of | Са"), in the vicinity of the cell 
plasma membrane to changes in membrane potential thus affecting voltage-gated ion 
channels and as a result the excitability of the cell. In this study, by combining fluores- 
cence confocal imaging of |Са""), with simultaneous recording of whole-cell currents 
in voltage-clamped vascular myocytes we demonstrate that depending on the mem- 
brane voltage calcium sparks may evoke a burst of synchronous opening of either 
Ca’*-activated СГ channels ог Ca’*-activated K* channels thus leading to the dis- 
charge of spontaneous transient inward currents (STICs) or spontaneous transient 
outward current (STOCs), respectively. We show that comparative analysis of the 
amplitudes of calcium sparks and corresponding STOCs may provide information 
about colocalisation between Ca’*-release channels (ryanodine receptors, RyRs) and 
clusters ої Ca’*-activated К" channels. Propagating waves of elevated | Са""), observed 
during membrane depolarization seem to arise from spatio-temporal recruitment of 
local Ca’*-release events. Similar to visceral myocytes, in vascular myocytes the spa- 
tial nonuniformity of sarcoplasmic reticulum and RyRs distribution within the cell 
may account for the existence of ‘frequent discharge sites’ and the wide variation in the 
Ca^ wave propagation velocities observed. 

Keywords: calcium spark, calcium wave, spontaneous transient outward currents, 
spontaneous transient inward currents, patch clamp, smooth muscle cell, laser scan- 
ning confocal microscopy 


Introduction 


Changes in intracellular free Са"" concentration ([Са?*]) play а central role in 
the control of diverse activities of the cell by regulating various enzymes, the con- 
tractile proteins and membrane ion channels. The latter, Ca’*-dependent ion 
channels, are of special interest since they couple changes іп | Са""|, to changes in 
membrane potential thus affecting voltage-gated ion channels and as a result the 
excitability of the cell. Many of these Ca^-activated channels are also volt- 
age-sensitive so that, depending on the membrane potential, different types may 
become responsive to changes in [Ca^]. Further Ca” signal specificity arises as 
the Са" signal is highly structured both in space and time (e.g. Са" “sparks”, 
“blips”, “puffs”, waves, oscillations [1]). Important biological information is en- 
coded in the various spatial and temporal characteristics which is translated by 
Ca’*-regulated effectors into specific responses. For example, 


44 “Physics of the Alive”, Vol. 8, No. 1, 2000 


Modulation of the Membrane Cunductance by Calcium Release Events in Vascular Myocytes 


calmodulin-dependent protein kinase II has been recently identified as a Ca” os- 
cillation frequency decoder [2]. 

More than ten years ago spontaneous transient outward currents (STOCs) 
arising form the synchronous opening of up to hundred Са" -асбіуатед potassium 
channels were detected in visceral and vascular myocytes and it was suggested 
that STOCs are triggered by periodic release of small amounts of Са" from the 
sarcoplasmic reticulum (SR) [3]. Later spontaneous transient inward current 
(STIC) was identified in myocytes from trachea [4] and portal vein [5] as a cur- 
rent carried through Са" -асбіуатед chloride channels. More recently spontane- 
ous Ca’*-release events were identified in both cardiac [6] and skeletal muscle [7] 
by means of Са" sensitive dyes introduced into cells. These were termed ‘sparks’ 
on account of the flash of fluorescent light they produced [6]. STOCs observed in 
single smooth muscle cells were attributed to the occurrence of Ca^ sparks when 
later the latter were first recorded in vascular myocytes, although simultaneous 
recordings of STOCs and sparks were not at that time made in the same cell [8]. 
We have more recently described spontaneous Са" release events and waves of 
elevated [Ca^], in single smooth muscle cells (SMCs) from the longitudinal layer 
of the small intestine [9]. In this paper we describe the results of experiments 
where membrane current was recorded in voltage-clamped single smooth muscle 
cells from small mesenteric arteries while changes in | Са" |, were monitored si- 
multaneously by confocal imaging of fluo-3 fluorescence. 


Materials and Methods 


Experimental procedures were generally the same as previously described 
[9]. Briefly, experiments were performed at room temperature on single smooth 
muscle cells isolated from small branches (diameter about 200 um) ofthe anterior 
mesenteric artery of male guinea-pigs (350-400 g) after collagenase and elastase 
(0.1% and 0.01%, w/v) treatment at 36°C for 25-35 min. The myocytes were 
loaded with the fluorescent Ca” -sensitive indicator fluo-3 by dialysis of the cell 
through the patch pipette with solution containing (in mM): KCl, 96.5; NaCl, 4; 
MgCl, 5; Ма, АТР, 1; creatine, 5; glucose, 20; Hepes, 10; fluo-3 (pentapotassium 
salt), 0.1; pH adjusted to 7.4 with KOH.. To ensure equilibration between the pi- 
pette solution and cytosol, electrical recordings and fluorescence imaging were 
commenced at least 5 min after establishment of the whole-cell configuration. 
The composition of the bath solution was (in mM): NaCl, 120; KCl, 6; CaCl, 2.5; 
MgCl, 1.2; glucose, 12; Hepes 10; pH adjusted to 7.4 with NaOH. 

Experimental chambers containing cells were placed on the stage of an 
Axiovert 100 inverted microscope attached to a LSM 410 laser-scanning unit 
(Zeiss, Oberkochen, Germany). Cells were imaged using a Zeiss 
plan-Apochromat 63x1.4 N.A. oil-immersion objective. The pinhole was set to 
provide an х-у resolution of 04um and a confocal optical section (z resolution) of 
0.8-1.0 um (measured with 0.2 um fluorescent beads). 

Image processing was carried out using an Indy workstation (Silicon Graph- 
ics, Inc., Mountain View, California, USA) with custom routines written in IDL 
(Research Systems, Inc., Boulder, Colorado, USA). Where appropriate, data are 
expressed as mean values + S.D. for the number of cells (п) analysed. 
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The contribution of the zero-frequency component of the power spectra of a 
line scan image was removed by subtraction of the average normalised fluores- 
cence from the signal at each pixel along time profiles crossing sites of interest. 
For each site , 3 power spectra (from 3 adjacent pixels along the scan line) were 
calculated from the Fourier transform of the pixel data and averaged, and dis- 
played as a low-frequency (0-20 Hz) noise spectrum. 

Membrane currents were measured using conventional whole-cell patch 
clamp technique [10]. Fire-polished patch pipettes filled with pipette solution 
(see above) were connected to the head stage of an Axopatch-1D (Axon Instru- 
ments) input amplifier and had free-tip resistance of 2-4 MQ. Liquid junction po- 
tentials were nulled with an offset circuit of the amplifier before seal formation. 
After formation of tight seal between pipette and cell membrane, fast capacitive 
transients were minimised using electrode capacitance compensation controls. 
Following formation of whole-cell configuration input resistance was evaluated 
from the change in steady current elicited by voltage stepping from -70mV to 
-60 mV. If input resistance exceeded 2 GQ further electrical recordings and fluo- 
rescence imaging were performed. No electronic compensation for cell capaci- 
tance, series resistance and leakage conductance was introduced. Electrical sig- 
nal was filtered at 1 kHz (-3 dB frequency) by four-pole low-pass Bessel filter. 
Voltage protocols were generated and electrical signals were digitised at 5 kHz 
using DigiData 1200 hosted by 486 PC running pCLAMP 6.0 software (Axon In- 
struments). To synchronise line-scan imaging with electrical recordings the ex- 
perimental voltage protocol was initiated by trigger pulse taken from line clock 
output of the LSM 410 laser-scanning unit. Accuracy of synchronisation was 
tested by detecting а 1 ms flash of light emitted froma xenon arc lamp triggered si- 
multaneously with voltage protocol. . 

The chemicals and drugs used in this study: collagenase (Type 1A), elastase 
(Type 1), soybean trypsin inhibitor (Type II-S), bovine serum albumin, 
N-2-hydroxyethylpiperazine-N'-2-ethanesulphonic acid (Hepes), dimethyl 
sulphoxide (DMSO), and cadmium chloride were obtained from Sigma Chemi- 
cal Co., Poole, Dorset, UK; fluo-3 acetoxymethy] ester, fluo-3 pentapotassium 
salt and Pluronic F-127 were obtained from Molecular Probes, Eugene, Oregon, 
USA; ryanodine and thapsigargin were obtained from Calbiochem, Beeston, 
Nottingham, UK. 


Results 


Freshly isolated smooth muscle cells from small mesenteric arteries were 
50-80 um long and 3-8 um in diameter. As we have previously demonstrated, 
these myocytes posses Ca^ -activated potassium channels and when [Ca^] is not 
clamped (0.3 mM EGTA in the patch pipette) exhibit spontaneous transient out- 
ward currents (STOCs) which are abolished in the presence of caffeine (4 mM) 
and blocked by ТЕА? [11, 12]. In the experiments described here the Ca’*-sensi- 
tive indicator fluo-3 was introduced into the cell by dialysis through the patch pi- 
pette with solution containing 0.1 mM offluo-3 (pentapotassium salt). The mem- 
brane potential was routinely clamped at -70 mV and stepped to the levels be- 
tween -30 mV and +10 mV for periods 2-3 s while membrane currents were 
recorded; at the same time confocal imaging of fluo-3 fluorescence was 
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performed by laser scanning unit along а line parallel and close to (within 1-2 шт) 
the plasma membrane. Successive fluorescent images of the scan line were then 
aligned from left to right so that the horizontal dimension of the line-scan images 
reflects time (increasing from left to right) and the vertical dimension shows the 
position along the scan line. 

Upon stepping to -30 mV or 0 mV STOCS were seen (Fig. 14) in 23 cells of 42 
myocytes examined. When a line-scan image was obtained during the period 
which included the depolarizing step, Са" -release events were usually seen 
(n8). In the example shown in Fig. 1 the scan line was set close (within 1-2 um) 
and approximately parallel to the cell plasma membrane (panel B) and scanned 
every 2 ms. The corresponding line-scan image is shown in Fig. 1C. As can beseen 
from the image, local Ca’*-release events were extremely variable in size and du- 
ration as in the case of intestinal smooth muscle cells [9]. The time course of nor- 
malised fluorescence at site iii (panel C) is shown in panel D. Corresponding trace 
of the whole-cell current is shown in panel E. It is apparent that each calcium 
spark (Fig.1C and D) was associated with a STOC (Fig. 1E). However, some 
STOCs were not associated with sparks or other Са" release events at site iii, but 
were likely to be associated with calcium release events elsewhere in the cell; for 
example, two small STOCs sandwiched between two large ones (Fig. 1E) corre- 
spond to calcium events just below the site iii (Fig. 1C) and appearing as small ele- 
vations of fluorescence in trace c (Fig. 1D). In such small cells (usually less than 
5 um in diameter) there was a high correlation observed between STOCS and the 
discharge of Ca’*-release events. In 23 myocytes 151 STOCs were observed of 
which 51 were associated with calcium sparks. Sparks without STOCs were not 
observed. In some cases, with membrane depolarization to 0 mV unitary calcium 
sparks could not be seen in line-scan image due to the loss of the image contrast 
that develops as [Ca^], increases [13]. However STOCS were seen in the corre- 
sponding records of whole-cell current (see Fig. 5). When STOCs were observed 
without a calcium release event in the line-scan record then it is likely that the 
scan line did not cross the site where the corresponding event was discharged. In 
the example shown in Fig. 1 sparks outlasted the STOCS associated with them. 

Besides well resolved calcium sparks, much smaller and faster local fluctua- 
tions in normalised fluorescence intensity are seen in the line-scan image shown 
in Fig. 2C (for example at site ii). Unlike microsparks detected in ileal myocytes 
(see [9]) these fluctuations occurred much more frequently (however they had 
similar amplitude, spatial size and kinetics) which complicates their resolution in 
the temporal domain (Fig. 2B). These fluctuations, however, were clearly distin- 
guishable from the background noise (Fig. 24) the power spectrum of which was 
essentially flat (see inset in panel А) as would be expected if the noise in the data 
arose from photon collection statistics (or ‘shot noise"). In contrast the power 
spectrum of the signal at site ii (Fig. 1C) revealed an excess noise in the 
low-frequency range (see inset in panel B, Fig. 2), suggesting that small local fluc- 
tuations in | Са"), were occurring (cf. [9]). Thus, similar to longitudinal intestinal 
smooth muscle cells, it seems possible that several microsparks discharging 
closely in time may give rise to a spark and that sparks were not randomly dis- 
charged throughout the cell but occurred at certain preferred sites or frequent 
discharge sites (FDSs). This is supported by the observation that most STOCs 
(Fig. 1E) were accompanied by a spark (Fig. 1C and D). Because STOCs were re- 
corded as currents through the whole of the cell membrane but the fluorescence 
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Fig. 1.Imaging of [Ca] with simultaneous recordings of net membrane current in voltage-clamped myocyte 
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from mesenteric artery. 

A, net membrane current recorded in single myocyte. The voltage protocol is shown on the top. Note 
spontaneous transient outward currents (STOCs) in whole-cell current record during voltage steps to 
—30 mV and 0 mV; potassium equilibrium potential was E к = -72 mV (see Methods for composition 
of the pipette and bath solutions). B, transmitted-light image of the same cell (patch pipette is seen on 
the right hand site). Black line in the cell image indicates the position of the line repetitively scanned 
during line-scan imaging. The transmitted-light image was taken after line-scan protocol was per- 
formed. C, confocal line-scan image formed by aligning (from left to right) the successive images of 
flou-3 fluorescence along the scan line (indicated in panel В), acquired at 500 Hz. Scale bars: x, 6 um; 
t, 300 ms. The image was obtained during voltage step to —30 mV. The fluorescence intensity was nor- 
malised to the average fluorescence intensity during the first 67 lines and grey-scale coded (bar below 
the image). D, the time course of the fluorescent signal from site iii (depicted in panel C). E, corre- 
sponding record of the net membrane current (at the membrane voltage of —30 mV). The correlation 
between the occurrence of calcium sparks and STOCs is emphasised by grey dashed lines. 
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Fig. 2.Fluctuations of the fluorescent signal at different sites of the line-scan image. 
The time course of the normalised fluorescence at two sites depicted as (i) and (ii) in Fig, 1 Cis shown 
in panels A and B, respectively. After subtraction of the average normalised fluorescence from the sig- 
nal at each of two sites, the power spectrum of the fluorescent signal was calculated from the Fourier 
transform of the line-scan data (see Methods). The insets show the low-frequency (0—25 Hz) noise 
spectrum at these two sites, respectively. Note an excess noise in the low-frequency domain at site ii 
(panel В). 


signal was measured only along one line within the cytosol, this finding suggests 
that during the observation period in this case the majority of sparks were dis- 
charged within the region crossed by the scan line. 

It should be noted that there was a very close correlation between the relative 
increase in | Са""|, and associated STOC amplitude when sparks occurred at the 
same spot along the scan line (e.g. Fig. 1). In the example shown in Fig. 3 the situa- 
tion is different. Local [Ca^], transients are seen during membrane depolariza- 
tion to -30 mV (Fig. 34). These trigger STOCS as can be seen from the corre- 
sponding whole-cell current record (Fig. 3C). When the temporal profiles of the 
normalised fluorescent signal at three sites in the line-scan image (bars a, b and c, 
panel A) were plotted (b, c and d in panel B, respectively) it reveals that the sec- 
ond Ca” spark consisted of two asynchronous Ca^ release events. It is notable 
that two of the three peaks of the corresponding STOC (panel C) occurred syn- 
chronously with the peaks in the time course of the fluorescent signal shown in 
panels b and c (as emphasised by arrows in Fig. 3). Taking into account that the 
scan line was aligned parallel to and within 1-2 um of the cell membrane it is 
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Fig. 3. Spontaneous transient outward currents (STOCs) triggered by Са?* sparks at membrane voltage of 
—30 mV. 
A, normalised line-scan image acquired at 500 Hz revealed several localised Ca?*-release events (Ca?* 
sparks). The scan line was positioned parallel to the surface of the cell. Scale bars: x=5 um, 7—300 ms. 
B, plots of the temporal profile of fluorescent signal at three sites labelled as a, b and c in panel A are 
shown in panels a, b and c, respectively. C, corresponding record of whole-cell membrane current. 
The correlation between the occurrence of STOCs and calcium sparks is emphasised by dotted lines. 
Note that asynchronous multiple Ca?*-release during second spark (panel A and plots c and d in panel 
B) is ‘sensed’ by Ca”*-activated potassium channels giving rise to *multi-peak' in the second STOC 
(panel D). The myocyte was dialysed with solution containing 0.1 mM fluo-3 (see Methods for com- 
position of pipette and bath solutions); E, = -72 mV. 


interesting that the first peak in the second STOC is much smaller than the sec- 
ond peak while the amplitude of both underlying Са"" release events was almost 
the same (compare b and c іп Fig. 3B). ). The most likely explanation for this find- 
ing is that the cluster of Са" -асіїуатед potassium channels giving rise to the 
STOC lies close to a SR Са" release site depicted by Багз а in the line-scan image 
(Fig. 34), while the site indicated by bars c is displaced from that cluster by more 
than 3 um. 

When the membrane voltage was held at -70 mV, spontaneous transient in- 
ward currents, STICs, (due to the opening of Ca’*-sensitive chloride channels) 
were observed in 7 cells of 42 myocytes examined. In 8 of 15 cases a STIC was 
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Fig. 4.Spontaneous transient inward currents (STICs) triggered by Са?* sparks at membrane voltage of 
—70 mV. 
A, line-scan image formed by stacking from left to right the sequential fluorescence images of the scan 
line acquired at 500 Hz and normalised to the average fluorescence during the first 100 lines. The scan 
line was positioned parallel to the surface of the cell. Two localised Ca -release events (calcium sparks) 
are seen on the image. Scale bars: х=5 um, 300 ms. В, the temporal profile of calcium sparks is illus- 
trated by a plot of the time course of the fluorescence change averaged over the region indicated by 
black bars in panel A. C, corresponding record of whole-cell membrane current showed two STICs. 
The composition of intra- and extracellular solutions was the same as in previous figure. Calculated 
equilibrium potential for chloride ions (E,,) was —5 mV. 


associated with а Ca^ spark (while no spark without a STIC was observed) an ex- 
ample of which із shown in Fig. 4. In this cell three Са’ sparks were observed over 
a period of 3 s of which two are shown on the image (panel A). The time course of 
normalised fluorescence at site indicated by bars in the line-scan image is shown 
in panel B and the corresponding trace of the whole-cell current is presented in 
panel C. It is apparent that each Са" spark is associated with a STIC. It should be 

noted that in this experiment membrane voltage (-70 mV) was held close to the 
potassium equilibrium potential (Е, (-72 mV) which completely eliminates any 
possible contribution of K* to the whole-cell current. As can be seen (compare B 
and C, Fig.4) Ca^ sparks outlasted STICs associated with them. Furthermore we 
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found that Ca^ sparks discharged at the same site (within the limits of our resolu- 
tion) may activate either a STIC or a STOC depending on the membrane voltage 
(n —3, data not shown) suggesting that both Ca?*-activates potassium and chlo- 
ride channels may lie in close vicinity to the same SR Ca’*-release site. Thus the 
membrane potential by altering the balance of Cl and K channel sensitivities to 
[Ca^] regulates which type of channel will open. 

We have recently demonstrated that in ileal SMCs spatio-temporal recruit- 
ment of ‘elementary’ calcium sparks may give rise to whole cell [Ca^] transient 
and/or propagating calcium wave [9]. In SMCs from mesenteric artery the depo- 
larizing step to 0 mV in some cells gave rise to a calcium wave which spread 
throughout the myocyte. It was possible to visualise the initiation site of such 
waves on occasions by line-scan imaging. In the example shown, the record of 
membrane current revealed numerous STOCS (Fig. 5/4) but associated sparks 
could not be distinguished clearly due to the large increase in global fluorescence 
(Fig. 5B) which lead to the loss of image contrast that develops as | Са""| increases 
[13]. Examination of the rising phase of [Ca’*],.on an expanded time scale 
(Fig. 5C) showed that initially (at very beginning of the voltage step) there was a 
relatively uniform increase in fluorescence along the scan line which culminated 
after about 80 ms in a Ca’*-release event close to the middle of the scan line (indi- 
cated by arrow in panels C and D). From the initiating Ca’*-release event a cal- 
cium wave propagated at a velocity of about 60 um.s and then more slowly. The 
relatively uniform rise in | Са""|, at early times may reflect calcium entering the 
cell through voltage-dependent channels. This apparently initiates a spark-like 
calcium release at a single point along the scan line. However, it is curious that the 
calcium wave propagates from this point at a high velocity which implies that cal- 
cium-induced calcium release (CICR) is involved, although CICR did not occur 
from the sites involved in the calcium wave during the initial period when calcium 
was entering through voltage-dependent channels. Thus the CICR sites triggered 
by calcium from the calcium store were not, apparently, triggerable by calcium 
entering through the plasmalemma. A shaded surface plot of the image from 
panel С (Fig. 5D) illustrates that fluctuations in the fluorescence signal during the 
calcium wave are significantly larger than those before the voltage step was ap- 
plied and those at sites where the wave had not invaded. It seems likely that the 
calcium wave consisted of numerous, local calcium release events varying in am- 
plitude and size (Fig. 5D). 


Discussion 


Localised spontaneous |Са" |, transients were first predicted in smooth mus- 
cle cells more than ten years ago when spontaneous transient outward currents 
(STOCs) reflecting the clustered opening of Ca’*-activated potassium channels 
were detected in single smooth muscle cells under voltage clamp [3]. It was there- 
fore of interest to explore whether calcium sparks are responsible for STOCs acti- 
vation in smooth muscle by simultaneous recording of whole-cell current and im- 
aging of [Ca^]. We have found that in the small vascular myocytes Ca^ spark dis- 
charge was well correlated with STOC occurrence (see Fig. 1). It is perhaps 
fortuitous that such correlation was obtained since it would appear that, with few 
exceptions, most of the Ca’*-release events in the cell during the depolarizing step 
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Fig. 5. Net membrane current and [Ca**] wave elicited by voltage step to 0 mV. 


A, net membrane current in response to voltage step to 0 mV (voltage protocol is shown on the top). 

B, corresponding line-scan image. The scan line was oriented parallel to the plasmalemma. The fluo- 
rescence intensity was normalised to the average fluorescence intensity during the first 1024 lines. C, 
the leading edge of the calcium wave (from panel B) is presented on an expanded time scale, enhanced 
magnification and different coding of the fluorescence intensity (grey-scale bar on the left). Corre- 
sponding voltage protocol is shown at the bottom. Propagation velocities estimated from the slopes of 
lines drawn along the wave front are: vl = 60 um · s; v2 = 10 um · 8... For panels В and C scale bars 
are: x, 5 um; г, 400 ms. D, image from panel C is shown as shaded surface plot. Dashed line on the im- 
age indicates the time when the cell membrane voltage was stepped to 0 mV. Note 'spark-like' event at 
the wave front (indicated by arrows in panels C and D). 
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were detected along the line of scan. This also implies that most frequent dis- 
charge sites (FDSs) in this cell were overlain by the scan line. It is unlikely, how- 
ever, that there are only a few, large clumps of SR in a cell which are capable of 
discharging calcium; perhaps each FDS is multiple and capable of entraining 
other, closely adjacent, sites on occasions even though these may not give rise to 
Ca**-release events spontaneously during the period of observation. Neverthe- 
less the non-uniform distribution of FDSs in smooth muscle (compared to stri- 
ated muscles) may reflect the irregular distribution of SR from which calcium is 
released. It is apparent (Fig. 1) that Ca’*-activated potassium channels may be 
considered as a very sensitive endogenous detector of the rise іп | Са""|, (see [9]). 
However, such a good correlation between STOCs and calcium sparks was not 
found in experiments on the smooth muscle cells from the cat oesophagus and 
gastric smooth muscle cells from Bufo Marinus [14, 15]. It should be noted, how- 
ever, that in myocytes from mesenteric artery correlation between the relative in- 
crease in [Ca^], during Ca^ spark discharge and associated STOC amplitude was 
not always as good as in the case when all sparks where occurred at the same spot 
along the scan line (Fig. 1) and this is illustrated by Fig. 2. Indeed, the STOC cor- 
responding to the peak signal shown in panel Bc is much smaller whereas the cor- 
relation between STOCs amplitude and the peak of Са?" signal shown in panels 
Ba and Bb is much better. It should be noted that the scan line was positioned 
close (1-2 um) and parallel to the membrane such that the distance for Са?* diffu- 
sion perpendicular to the membrane is roughly the same. As it was mentioned 
above, this discrepancy could be explained by variations in the distance between 
the cluster of Ca’*-activated potassium channels and the SR Ca’*-release site. 
Since the sparks we observed were abolished by thapsigargin and ryanodine (data 
not shown) it is clear that they are due to spontaneous SR calcium release via 
RyR activation. Similar to visceral myocytes, Ca’*-release events smaller and 
faster than calcium spark (which we termed ‘microsparks’ [9]) were observed in 
vascular SMCs (see Fig. 1). Evidence for the existence of Ca’*-release events 
smaller than sparks have been also obtained in cardiac myocytes [16] and skeletal 
muscles [17]. It is possible that a spark represents the discharge close in time of 
several microsparks and that continuous microspark discharge also underlies the 
much longer-lasting rises in [Ca^], which occurred (Fig. 1C). Itis unclear at pres- 
ent whether sparks and microsparks are two distinct populations of elemental 
events or whether they simply represent two extremes of a continuum of localised 
Ca^-release events in the cell. Although near the limit of resolution, microsparks 
could be clearly distinguished from ‘shot noise’ in the signal by the appearance of 
excess noise in low frequency domain in the power spectrum of the fluorescent 
signal (Fig. 2). 

When Са" sparks where first detected in vascular SMCs [8] it was proposed 
that STOCs triggered by local [Ca’’], transients may produce hyperpolarization 
of the cell membrane which might be sufficient to turn of Са?" ‘window’ current 
thus leading to decrease in global |Са""), (which, however, has not been demon- 
strated) and myocyte relaxation. This proposal implies that a specific physiologi- 
cal role of Са" sparks is vasodilation. Our finding that Са" spark may trigger 
spontaneous transient inward current (STIC, [4, 5]) reflects that situation is more 
complex. Indeed, it is widely appreciated that chloride channels are a key deter- 
minant of membrane potential in variety of SMCs and their opening may lead to 
membrane depolarization and contraction of the myocyte [18, 19, 20]. Since the 
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same Са?" spark may activate either Ca’*-activated chloride channels or Са" -ас- 
tivated potassium channels (depending on the membrane potential) sparks 
would rather stabilise the membrane potential on certain level than lead to mem- 
brane hyperpolarization. Indeed, we found that at —70 mV STICs (but not 
STOCs) were triggered by Са?" sparks since this potential was close to the EK of 
—72 mV. With membrane depolarization STOCs became dominant since the 
electrochemical gradient was diminished for СГ but increased for К”; moreover 
Ca’*-activated potassium channels are known to become more sensitive to Са" 
with depolarization. 

STICs were found to decline much faster than corresponding Ca’*-release 
event (Fig. 4). It has been suggested that STIC decay is due to the intrinsic gating 
properties ої Ca’*-activated СГ channels [21, 22]. Our results show that the de- 
cline in [Ca"']i is slower than the decline of the STIC. This raises the question of 
why the СГ channels close while [Ca"'], still exceeds the threshold [Ca^], at which 
they opened, and why they are capable of a second opening within a few hundred 
milliseconds while [Ca^] is still elevated (Fig. 4). One possible explanation might 
be that Ca?*-activated СГ channels undergo desensitisation to Са" after their 
opening and recover their sensitivity after closing. This possibility should be fur- 
ther explored in the future experiments. 

We have previously demonstrated in ileal myocytes that spatio-temporal re- 
cruitment of ‘elementary’ Са" sparks may give rise to a calcium wave leadingtoa 
cell-wide elevation in | Са""), [9]. When a calcium wave was initiated by depolar- 
ization of a vascular myocyte (Fig. 5) fluorescence rose almost uniformly in 
subplasmalemmal area before a initiating spark occurred; the calcium wave then 
spread to what are, presumably, other adjacent discharge sites which had been 
previously refractory to discharge by entering Ca”. At later times ‘unitary’ cal- 
cium releases could not be well resolved within the calcium wave by fluo-3 emis- 
sion, but nevertheless the trace of net membrane current reveals STOCs which 
must reflect variations over time in subplasmalemmal [Ca^ ] (Fig. 5.4). These re- 
sults imply that there is a hierarchy of calcium discharge sites: some (‘initiator 
sites’) readily discharge spontaneously particularly if the cell is depolarized; of 
these some or all can be triggered by entering calcium, or calcium released from 
adjacent sites, to initiate a calcium wave; other calcium discharge sites (‘dormant 
sites’) are only revealed as contributors to a calcium wave. It is possible that initia- 
tor and dormant sites may reverse roles over longer periods of observation. 
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МОДУЛЯЦІЯ ПРОВІДНОСТІ МЕМБРАНИ ПОДІЯМИ ВИВІЛЬНЕННЯ КАЛЬЦІЮ У 
МІОЦИТАХ КРОВОНОСНИХ СУДИН 


Д.В. ГОРДІЄНКО 


В гладеньких м'язах Са" -залежні іонні канали є ефекторними молекулами які приймають участь 
у фізіологічно важливих регуляторних процесах, використуючих Са" як вторинний посередник. Ці Ka- 
нали трансформують зміну ІСа?"), поблизу плазматичної мембрани y зміну мембранного потенціалу, 
впливаючи таким чином на потенціалзалежні іонні канали i, як наслідок, на збудливіть клітини. Викори- 
стовуючи флуоресцентну конфокальну мікроскопію з одночасною реєстрацією трансмембранних іон- 
них струмів судинної гладеньком'язової клітини в умовах фіксації напруги ми встановили, що в 
залежності від рівня мембранного потенціалу кальцієві "спарки" можуть викликати сінхронне відкри- 
вання або Са""-активуємих СІ каналів, або Са'"-активуємих К“ каналів, наводячи відповідно до BHHHK- 


55 


D.V. GORDIENKO 


нення або спонтанних тимчасових вхідних струмів (STICs), або спонтанних тимчасових вихідних 
струмів (STOCs). Ми також демонструємо, що порівняльний аналіз амплітуд кальцієвих спарків i відпо- 
відних STOCs може бути використовано для отримання інформації про колокалізацію каналів вивіль- 
нення Са" (ріанодінвих рецеторів, RyRs) і кластерів Са""-активуємих К" каналів. Кальцієві хвилі, які 
розповсюджувалися по клітині у відповідь на деполярізацію мембрани, вірогідно виникали як наслідок 
просторово-часового інтегрування локальних подій вивільнення Са?*. Подібо до вісцеральних міоцитів, 
нявність "місць частого вивільнення" Са" (МЧВ) та варіабільність швидкоті розповсюдження каль- 
цієвих хвиль в міоцитах ковоносних судин вірогідно є наслідком неоднорідності розподілу SR та RyRs y 
клітини. 


МОДУЛЯЦИЯ ПРОВОДИМОСТИ МЕМБРАНЫ СОБЫТИЯМИ ВЫСВОБОЖДЕНИЯ 
КАЛЬЦИЯ В МИОЦИТАХ КРОВЕНОСНЫХ СОСУДОВ 


Д.В. ГОРДИЕНКО 


Вгладкихмышцах Са?*-зависимые ионные каналы являются эффекторными молекулами вовле- 
ченными в физиологически важные регуляторные процессы, исползующие Са“ как вторичный посред- 
ник. Эти каналы трансформируют изменения [Ca^], вблизи плазматической мембраны клетки в 
изменения мембранного потенциала, влияя таким образом на потенциалактивируемые ионные каналы 
и, врезультате, на возбудимость клетки. Используя флуоресцентную конфокальную микроскопию C од- 
новременной регистрацией трансмембранных ионных токов сосудистой гладкомышечной клетки в ус- 
ловиях фиксации напряжения мы демонстрируем, что в зависимости от величины мембранного 
потенциала кальциевые “спарки” могут вызывать синхронное открывание либо Са?*-активируемых СГ 
каналов, либо либо Са””-активируемых К“ каналов, приводя соответственно к возникновению либо 
спонтанных преходящих входящих токов (STICs), либо спонтанных преходящих выходящих токов 
(STOCs). Мы также демонстрируем, что информация о взаиморасположении каналов высвобождения 
Са" (рианодиновых рецепторов, RyRs) и кластеров Са""-активируемьхх K^ каналов может быть получе- 
на путем сравнительного анализа амплитуд кальциевых спарков и соответствующих STOCS. Распро- 
страняющиеся кальциевые волны инициируемые деполяризацией мембраны по-видимому являются 
результатом пространственно-временного интегрирования локальных событий высвобождения Ca”. 
Подобно висцеральным миоцитам, наличие “мест частого высвобождения” Са"" (MYB) и варибиль- 
ность скорости распространения кальциевых волн в миоцитах ковеносных сосудов вероято являются 
результатом неравномерного распределения саркоплазматического ретикулума и RyRs в клетке. 
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Abstract. 500 MHz NMR spectroscopy has been used to investigate the complexation 
of the anthracycline antibiotic daunomycin (DAU) with self-complementary 
deoxyhexanucleotides, 5'-d(CGTACG) and 5’-d(CGCGCG)), of different base se- 
quence in aqueous salt solution. 2D homonuclear 'H NMR spectroscopy (TOCSY 
and NOESY) and heteronuclear 'H-"P NMR spectroscopy (HMBC) have been used 
for complete assignment of the non-exchangeable protons and the phosphorus reso- 
nance signals, respectively, and for a qualitative determination of the preferred bind- 
ing sites of the drug. Analysis shows that DAU intercalates preferentially into the 
terminal d(CpG) site of each of the hexanucleotides and that the aminosugar of the 
antibiotic is situated in the minor groove of the hexamer duplex, partly eclipsing the 
third base pair. 

A quantitative determination of the complexation of DAU with deoxyhexanucleoti- 
des has been made using the experimental concentration and temperature 
dependences of the drug proton chemical shifts, which have been analysed in terms of 
the equilibrium reaction constants, limiting proton chemical shifts and 
thermodynamical parameters (enthalpies, AH, entropies, AS) of different drug-DNA 
complexes (1:1, 1:2, 2:1, 2:2) in aqueous solution. It is found that there is an apparent 
site-specificity on drug binding with a triplet nucleotide sequence containing two adja- 
cent C-G base pairs flanked by an A-T base pair in the 5’-d(CGTACG) hexamer over a 
triplet with three successive C-G base pairs as in 5’-d(CGCGCG). The major reason 
for the preferential binding of DAU is found to stem from hydrophobic interactions of 
its aminosugar moiety in the minor groove of the double helix. 

The most favourable structures of the 1:2 complexes have been derived from the in- 
duced limiting proton chemical shifts of the drug in the intercalated complexes with 
the hexanucleotide duplex, in conjunction with 2D NOE data, which are in good 
agreement with the X-ray results [28-30]. 

Keywords: Daunomycin, drug-DNA complexation, self-complementary deoxyhexa- 
nucleotides, base sequence, 'H/'P NMR, thermodynamical and structural parame- 
ters. 


Introduction 


Daunomycin (DAU), is a prominent member of the anthracycline antibiot- 
ics, a class of compounds that has well documented anti-cancer properties [1-3]. 
The chemotherapeutic efficacy of DAU has been attributed to its intercalative 
binding with DNA, resulting in inhibition of DNA replication and RNA tran- 
scription [1,3]. In contrast to some other intercalating drugs, it is thought that 
DAU does not exhibit significant sequence-specificity on interaction with DNA 
[4]. Data about the selective binding of DAU with synthetic deoxyoligonucleoti- 
des and native DNA are contradictory as well as theoretical investigations [5-12]. 

The use of native DNA to investigate the selectivity of binding of a drug is lim- 
ited because of the complexity of the structure and conformational variability of 
DNA molecules and the large variety of drug binding sites in a polymer chain. On 


58 "Physics of the Alive", Vol. 8, No. 1, 2000 


NMR Investigation of Sequence-Specificity of Daunomycin Binding with Deoxyhexanucleotides 


the other hand, it has been shown [13-16] that the distinctive structural features of 
drug-oligonucleotide complexes and the specificity of binding of aromatic lig- 
ands, such as proflavine, ethidium bromide and actinomycin D, with nucleic acids 
is manifested in short, self-complementary oligonucleotides. A more complete 
understanding of the nature of the physical interactions responsible for molecu- 
lar complexation can be achieved by consideration of the thermodynamics of mo- 
lecular complexation in aqueous solution [14]; in particular, for systems investi- 
gated under the same solution conditions, conclusions about the selective binding 
of low molecular weight aromatic compounds with defined sites in the 
oligonucleotide sequence can be drawn by comparison of the equilibrium con- 
stants of ligand complexation with oligonucleotides of different length and base 
sequence in the chain [14,17]. Using short oligonucleotides in solution, there is a 
complex dynamic equilibrium of single-stranded and double-stranded 
conformational states of the oligonucleotide chains and their complexes with 
drugs [13,14]. 

In order to investigate quantitatively the sequence specific binding of the an- 
thracycline antibiotic with double-stranded oligonucleotides, the complexation 
of DAU to the self-complementary deoxyhexanucleotides, 5’-d(CGTACG) and 
5’-d(CGCGCG), in aqueous salt solution, has been studied by 1р an 2D 
(500 MHz) NMR spectroscopy. The presence of identical d(CpG) terminal sites 
in the hexamers studied, provides an opportunity to determine the possible effect 
of A-T or C-G base pairs on the affinity of binding of DAU with the triplet binding 
sites in the DNA double helix. 


Materials and methods 


The structure of the DAU molecule with the numbering of all the proton posi- 
tions used in the NMR study is presented in Fig. 1. DAU was purchased from 
‘Fluka’ Chemical Company and used without further purification. The drug was 
lyophilised from D,O and re-dissolved in deuteriated 0.1 moll" phosphate buffer, 

рр 7.1, containing 10" moll’ EDTA. The 

concentration of antibiotic stock solution 

was measured spectrophotometrically 

O OH O using the molar absorption coefficient 
T tere €=115001mol™ cm” (A=480 пт) [18,19]. 
р Be Al OH The deoxyhexanucleoside penta- 

3 i i 8 phosphates, 5’-d(CCTACG), and 
5’-d(CGCGCG), were synthesised by the 
Me оо OH ‘Oswel DNA Service’ (University of 
о Southampton, UK), lyophilised from 

99.95% D,O and re-dissolved in 0.1 mol Г 
phosphate buffer, pD 7.1, containing 


T ere, 10" mol Г EDTA. 
ber до 1D and 2D NMR measurements were 
NH, made on a 500 MHz Bruker DRX spec- 
OH trometer with the residual water peak satu- 


rated during relaxation. Proton chemical 


Fig. 1. Structural formula of Daunomycin (DAU). shifts were measured relative to an internal 
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Table 1. "Р Chemical shifts of 5’-d(CGTACG) and 5’-d(CGCGCG) (8,, ppm)?, 
Daunomycin+5’-d(CGTACG) and Daunomycin+5’-d(CGCGCG) (5,*) and changes 
іп УР chemical shifts on binding of Daunomycin to the hexamers (Дб, = 8,* -6,) 


| ROO OE 4. Ro ИЕ Е | B 

| 5’-d(CGTACG)+Daunomycin 

|, | 261 | as | 27 | 26 | am | 
| a | an | 254 | 242 | 283 | 155 | 
| as, | -005 | 4032 | 4033 | -015 ^ 4097 | 
| 5’-d(CGCGCG) + Daunomycin | 
| & | 257 | 2732 | 260 | 2616 | 2513 | 
| зи | 259 | 25 257 | 252 | 220 | 
| Ab, | 90056 | 90142 | 40038 | 40096 | 40313 | 


а) "Р chemical shift referenced with respect to 0.1 mol I“ phosphate buffer signal. 


reference, TMA (tetramethylammonium bromide). "P NMR chemical shifts 
were referenced with respect to the "Р signal of the 0.1 mol I” phosphate buffer. 
Methods for preparation of solutions and the NMR experimental procedures are 
described previously [14,20]. 2D TOCSY and NOESY experiments were carried 
out at 7=303 or 308 К using standard pulse sequences described previously 
[14,17]. HMBC (heteronuclear multiple bond coherence) spectra [21] have been 
used to analyse the "Н-"Р spin system, using 8K data points in the F2 dimension 
(^P) and 256 scans in the F1 dimension (H); to improve spectral quality, a sinu- 
soidal apodisation function was applied to both dimensions before Fourier 
transformation. 


Results 


NMR assignment of all the non-exchangeable protons of the deoxyhexanu- 
cleotides, 5’-d(CGCGCG) and 5’-d(CGTACG) was made using homonuclear 
2D TOCSY and 2D NOESY measurements. The equilibrium parameters of 
self-association of 5’-d(CGCGCG) and 5’-d(CGTACG) in aqueous solution 
were determined earlier [22] under the same experimental conditions. 

Analysis of the NOE spectra of the mixed solutions shows a qualitatively simi- 
lar DNA conformation (B-form) when complexed with DAU as with the 
tetranucleotides [23]. Changes in the "Р spectra (Table 1, especially to the termi- 
nal phosphate Р.), as well as observations of intermolecular NOE contacts 
(Table2 , especially to the terminal G6, C5 and C1), provide evidence that DAU 
intercalates into the terminal d(CpG) sites of the hexanucleotide duplexes. The 
same situation was observed for complexation of DAU with deoxytetranucleoti- 
des [23]. 
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Table 2. Intermolecular NOE contacts between non-exchangeable protons of 
Daunomycin and 5’-d(CGTACG) (D,=1.33 mmol Г', №=2.66 mmol Г") and 
5'-d(CGCGCG) (D,=0.96 mmol I!, №=2.63 mmol I), 7-323 К, рр 7.1 in 


Proton of Daunomycin | Proton of 5 -d(CGTACG) | Proton of 5 -d(CGCGCG) 
(Chemical Shift, ppm) 


(Chemical Shift, ppm) 


| (Chemical Shift, ppm) 


| 
| 
| 
| 


НГ (5.41) 


| H1 G6 (6.09) 

| HI C5 (5.55) 

| H4' G6 (4.21) 

| H5' G6 (4.04) 

| H5" G6 (3.85) 

| H2//H2" C5 (2.40) 


| H5' G6 (4.06) 
| H5" G6 (4.03) 
| H4' G6 (4.19) 


| 


| H7 (4.76) | H5” C5 (3.79) 
| НГ G6 (6.09) 
HS’ (4.38) НІ" ТЗ (5.63) 
H4 (3.81) | H2’/H2” C5 (2.40) 


| H3’ C5 (4.89) 


OCH, (3.86) 


HS’ C1 (3.58) 


| H5 C5 (5.57) 


НІ G2 (5.62) 


НӨ C5 (7.30) 


! H3 СІ (4.67) 
| H5 C5 (5.32) 


| H10e (3.04) 


| НІ" G2 (5.62) 


| H8a (2.00) 


| НУ" G6 (5.85) 


HS’ G6 (4.04) 


| H4'/H5' C5 (4.11) 
На G6 (4.21) 
НГ G6 (6.09) 
| НУ А4 (4.14) 


| H1’ G2 (5.85) 


H2' (2.10) 


H5" G6 (3.85) 
H5' G6 (4.04) 


| H4'/HS'C5 (4.11) 
| H4'G6 (4.21) 
| HI'G6 (6.09) 


НУАА (4.14) 


| H2" (1.87) 


H5" G6 (3.85) 


| H5’ G6 (4.04) 
| HA'/H5' C5 (4.11) 
НА G6 (4.21) 
| HI’ G6 (6.09) 


H5' A4 (4.14) 


| H2" C5 (2.34) 


| CH, (1.29) 


H2' C5 2.40) 


| H4' A4 (4.38) 

| HA'/H5' C5 (4.11) 
| НУ A4 (4.14) 

| H5" A4 (3.96) 
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Fig. 2. Experimental dependence of chemical shifts of some of the non-exchangeable protons of DAU in 
mixed solution with the deoxyhexanucleotides: a) with concentration of the hexamer; 
[5' - d(CGTACG) + DAU, D, = 0.19 mmol/l; 5: — d(CGCGCG) + DAU, D, = 0.41 mmol/l]; and 
b) with temperature; [5'- d(CGTACG)+ DAU, D, = 0.19 mmol/l, N, = 0.41 mmol/l; 
(5' - d(CGCGCG) + DAU, D, = 0.13 mmol/l, М, = 0.44 mmol/l]. 


A. qualitative comparison of the experimental curves for some of the 
non-exchangeable protons of DAU in the mixed solutions with deoxyhexanu- 
cleotides indicates a stronger binding of the drug with S’-d(CGTACG) than with 
5’-d(CGCGCG) as is seen from the slightly higher curvature of the experimental 
concentration and temperature dependences of proton chemical shifts of DAU 
in the mixture with 5’-d(CGTACG) (Fig. 2). In order to estimate quantitatively 
the complexation of DAU with deoxyhexanucleotides, a similar approach and 
model was used as for the complexation of DAU with tetranucleotides [23]. 

The following equilibrium reactions were considered: 


D + D e р, (a) N +N e——, М, (b) 
D + N < DN (c) D +N, «= DN, (d) (1) 
D + DN <» D,N (е) DN,-D«——5. D,N, (f) 
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Table 3. Calculated values of limiting proton chemical shifts (8, , ppm) and equilib- 


rium constants (К, 1 mol!) for complexation of daunomycin and the hexamers, 
5’-d(CGTACG) and 5’-d(CGCGCG) at 303 K. 


Proton| 8, | 8, | 8, | 8, | 8, |K, 1 | K, 10. | K, 10 | K,, 10 | 


m 


| Daunomycin + 5'-d(CGTACG) 
ні |771/7.78 |7.41| 744 | 1.83. | | 

H3 (680 7.29 [7.21 692 7.55 aiii pe d 
HI’ |4.82 5.07 5.45 479 5.52. | | 
| OCH, ЕЯ 3.89 3.79 3.31 | 4.02 | | | 
“Hilde 1234 294 2.89 230 3.05 

Н10а 2.60 2.56 2.51 2.52 2.81 
| А Daunomycin + 5'-d(CGCGCG) 
ні 7.44 7.55 |743 7.25 7.83 6849 22004280 3.1528 49+12 
H2 7.18 772 7.61 7.06 | 7.78 | | | | 
| H3 1698 7227 722 6.84 | 7.55 | 
OCH, 351 13.80 3.76 3.38 4.02 | | | | 
| Hl0e 1295 | 3.00 2.87 2.35 | 3.05 | | | | | 


| 
| 


where D and М are concentrations of the monomers of DAU and hexanucleotide, 
respectively, K, is the equilibrium constant for dimerization of the drug, Ky is the 
equilibrium constant for formation of the hexamer duplex, and K,-K, are equilib- 
rium constants for formation of 1:1 (DN), 1:2 (DN,), 2:1 (D,N) and 2:2 (D,N,) 
complexes, respectively. Under the conditions of the experiment, self-association of 
DAU is confined to dimerisation of molecules in which the dimerisation constants, 
K,, is equal to 290 and 188 І mol’ at Т=303 К and 313 K, respectively [24]. 


Using this scheme, the observed proton chemical shift of the antibiotic may 
be expressed in the form [14,17]: 


8 -(D/D)G, -2K,D8, + K,Nà; +KxK, N’ 8, + 


+2K,K,DN 5, +2K КК, DN? 5,), МУ 
where 5,, and, are the chemical shifts of the drug in the monomer and dimer and 
5, 79, are the limiting values of the proton chemical shifts of the antibiotic in the 1:1, 
1:2, 2:1 and 2:2 complexes; D, is the initial molar concentration of the drug. Minimi- 
sation of the discrepancy function has been carried out using eight parameters 
(К, - K, and6, -8,) employing the numerical procedure of determination of pa- 
rameters published previously [14,17]. The calculated values of equilibrium parame- 
ters,5, —ó, and К, — K,, are presented in Table 3 for the systems studied, together 
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Table 4. Thermodynamic parameters (AG°, AH’, kJ mol, AS?, J тої"! К of 
daunomycin complexation with the deoxyhexanucleotides 5'-d(CGTACG) and 
5’-d(CGCGCG). 


| Reaction | -AG° | -AH? | -AS? 

| Daunomycin + 5’-d(CGTACG) 

Іремеюьоу | 24 | 160520 | — 45056 
|D«N,—55DN, | 4044 | 5346 | 45+16 
|D«DN —5—D,N | — 28t5 | 941 | 210428 
[DN *D-E— DN, әбу | ааб | 50 

| Daunomycin + 5'-d(CGCGCG) 

РУМ ЪрМ | 306 | | 175324 | 484451 
|D«N,—5—DN, | — 38M — | и | 201533 | 
|D4DNC5SDN 21424 88:13 — | — 225434 | 
[DN +D EDN] agag | 45 | зом | 


L | M | 


with the errors resulting from averaging the data from nine protons. In the calcula- 
tions, the data on concentration (at two different temperatures) and temperature 
dependences of proton chemical shifts of DAU in the mixed solutions (Fig. 2) have 
been jointly processed. 

The calculated parameters for DAU complexation with the hexamers at 
T=303 K are presented in Table 3. 

Comparison of the calculated equilibrium constants, K, aK, for 1:1, 12,24 
and 2:2 complex formation of DAU with deoxyhexanucleotides (Table 3) are 
qualitatively similar for both systems studied. Thus, the probability of 1:2 com- 
plex formation of DAU with the hexamers (K,) 15 much higher than formation of 
any other drug-hexanucleotide complex in aqueous solution. The magnitudes of 
equilibrium constants K, for 2:2 complex formation are relatively low (Table 3), 
in spite of the potential existence of two triplet binding sites in the hexamer se- 
quence, i.e. complexation of the second DAU molecule with the hexamer duplex 
is obviously anti-cooperative in solution. This is probably due to the distinctive 
features of the binding of DAU to DNA, as it contains a bulky, positively charged 
aminosugar ring and hinders further binding of the antibiotic in solution. Com- 
parison of the calculated K, values for the two systems studied shows that a triplet 
containing two adjacent C-G base pairs flanked by an А-Т base pair in the 
5’-d(CGTACG) sequence is a more preferable site for DAU binding toa DNA 
duplex than a triplet with three successive C-G base pairs, as in 5-d(CGCGCG). 
It follows that DAU shows some site specificity when С.С is replaced by an АТ 
base pair in the triplet nucleotide sequence, in agreement with footprinting titra- 
tion experiments [9] and theoretical investigations [12]. 
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Thermodynamic Parameters of Complexation. Quantitative assessment of 
the thermodynamical parameters for each type of complex requires the determi- 
nation of the temperature dependences of their mole fractions [14,25]. This can 
be accomplished by analysis of the experimental dependences of proton chemical 
shifts of the interacting molecules (Fig. 2b), taking into account the additive con- 
tributions of all associated forms in solution [14,25]. 

The mean values of Gibb’s free energies, enthalpies and entropies deter- 
mined for complexation of DAU with deoxytetranucleotides are summarised in 
Table 4. 

The thermodynamic parameters, summarised in Table 4, show that all the re- 
actions of complex formation of DAU with the hexamers are exothermic. ДС", 
АН? and AS" values for the formation of 1:1, 2:1 and 2:2 complexes of DAU with 
the two hexamers are similar to each other, within experimental error, for 
5’-d(CGCGCG) and 5’-d(CGTACG). However, a large difference in the ther- 
modynamic parameters of 1:2 complex formation is observed, indicating appar- 
ent distinctions between the binding affinity in the double-helical triplet sites of 
the hexanucleotides studied. It is seen from Table 5 that the absolute values of 
enthalpy and entropy for complexation of DAU with 5’-d(CGTACG) are sub- 
stantially smaller than ДН" and 45° for 1:2 complex formation with 
5’-d(CGCGCG) (the differences in АН" and AS" for the two hexamers are ap- 
proximately two and four times, respectively). Thus, it may be concluded that hy- 
drophobic interactions, which lead to positive entropy contributions, play a 
much more significant role on drug binding to S’-d(CGTACG) compared to the 
5’-d(CGCGCG) oligomer. 


Fig. 3. Calculated structure of the 1:2 complex of DAU and 5’-d(CGTACG); the antibiotic chromophore is 
shaded. 
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Structures of the DAU-deoxyhexanucleotide complexes. Quantitative deter- 
mination of the structures of DAU with the deoxyhexanucleotide duplexes (the 
1:2 complex) has been made using X-PLOR software [26] taking into consider- 
ation intra- and inter-molecular NOE contacts in conjunction with the calculated 
values of limiting proton chemical shifts (Table 3). A modified model of the 
equivalent magnetic dipoles [14,17], based on approximation of quan- 
tum-mechanical iso-shielding curves of nucleic acid bases [27], has been used to 
calculate, as a function of the conformational parameters of the double helix, the 
shielding of the non-exchangeable protons of the DAU chromophore by adjacent 
base pairs in the intercalated 1:2 complex. 

The most favourable calculated structure of the fragment of the 1:2 complex 
of DAU with deoxyhexanucleotide, 5’-d(CGTACG) corresponding to intercala- 
tion of the drug into the 5’-d(CG) site of the hexamer is shown as an example in 
Fig. 3. 

The intercalated 1:2 complex of DAU with 5’-d(CGTACG) is characterised 
by the following helix parameters: D, =0.62 пт, О =34°, @=-59°, т =-64°, 
D, =009пт, D, =0.20 пт, р =3.6°, к= 20°. The 1:2 complex of DAU with 
5’-d(CGCGCG) is described by similar parameters of the double helix within er- 
ror limits. The calculations have shown that in the intercalated complex the anti- 
biotic chromophore is oriented at an angle v (у, = 120° for 5’-d(CGTACG) and 


у, = 130° for 5’-d(CGCGCG)) to the base pair axis in the double helix of the 


hexamer duplex, which is similar to that found in the crystal structures of the in- 
tercalated complexes [28-30]. The differences in the y values for intercalation of 
DAU into the hexamers studied may be due to the characteristics of the location 
of the amino sugar residue of DAU in the minor groove of the DNA double helix 
in different base sequences. At the same time it may be concluded that the struc- 
tural conformational parameters of the intercalated complexes of DAU with 
deoxyhexanucleotides, which have been found to be practically the same, within 
error limits for the two systems studied, are determined, in the main, by the base 
sequence (5'-d(CG)) at the intercalation site. 

In conclusion it should be emphasised that: 

(i) Quantitative analysis of the complexation of DAU with the 
deoxyhexanucleotides, 5'-d(CGCGCG) and 5’-d(CGTACG), shows an apparent 
site specificity on drug binding with a triplet nucleotide sequence containing two 
adjacent C-G base pairs flanked by an A-T base pair in the 5’-d(CGTACG) 
hexamer over a triplet with three successive C-G base pairs as in 5’-d(CGCGCG); 

(ii) DAU shows some sequence specificity on binding with d(CpG) sites in 
single-stranded deoxyhexanucleotides; 

(iii) Hydrophobic interactions of the aminosugar residue in the minor groove 
of the double helix are the major effects in the preferential binding of DAU with 
triplet sequences containing an A-T base pair in 5’-d(CGTACG) compared with 
three successive С.С base pairs in the 5’-d(CGCGCG) hexamer in aqueous 
solution. 
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ДОСЛІДЖЕННЯ СИКВЕНС-СПЕЦИФІЧНОСТІ ЗВ'ЯЗУВАННЯ ДАУНОМЩИНУ 
3 ДЕЗОКСИГЕКСАНУКЛЕОТИДАМИ МЕТОДОМ ЯМР-СПЕКТРОСКОПП 


ОН. ВЕСЕЛКОВ, В.І. ПАХОМОВ, Л.Н. ДИМАНТ, Р.ДЖ. ITOH, Д.Б. ДЕВІС 


Вивчено комплексоутворення антрациклинового антибіотика дауномщину (DAU) з самокомп- 
лементарними дезоксигексануклеотидами 5'-d(CGTACG) i 5’-d(CGCGCG) різної послідовності основ 
у ланцюгу у водно-сольовому розчині методом ЯМР-спектроскопії (500 МГц). Двовимірні гомоядерна 
'Н-ЯМР-спектроскопія (TOCSY i NOESY) та гетероядерна "Н-"Р ЯМР-спектроскопія (HMBC) вико- 
ристано для повного віднесення резонансних сигналів необмінних протонів та ядер фосфору, відповід- 
но, а також для якісного визначення місць переважної посадки антибіотика при зв'язуванні з 
гексамером. Встановлено, що DAU інтеркалює у термінальні d(CpG)-caitru гексануклеотидів, при цьо- 
му аміноцукрова група антибіотику розміщується в малій бороздці гексамерного дуплексу, частково пе- 
рекриваючи третю пару основ. 

Виміряно концентраційні і температурні залежності протонних хімічних зсувів молекул, на ос- 
нові котрих розраховані рівноважні константи реакції комплексоутворення DAU з дезоксигексанукле- 
отидами, граничні значення хімічних зсувів протонів дауноміцину і термодинамічні параметри 
(ентальпія AH і ентропія AS) для різних комплексів антибіотик-ДНК (1:1, 1:2, 2:112:2) уводному розчині. 
Встановлено, що має місце означена специфічність зв'язування антибіотику з триплетною нуклеотид- 
ною послідовністю, яка має дві суміжні С-С-пари основ, фланковані А-Т-парою основ y гексамері 
5’-d(CGTACG), в порівнянні з триплетом, який складається з трьох послідовних С.С-пар основ y гекса- 
мері 5’-d(CGCGCG). Аналіз отриманих результатів показав, що специфічність зв'язування дауноміци- 
ну з ДНК зумовлено гідрофобними взаємодіями аміноцукрової групи антибіотика у малій бороздці 
подвійно спіралі. Базуючись на розрахованих значеннях індукованих протонних хімічних зсувів анти- 
біотика у складі інтеркальованих комплексів та даних 2М-МОЕ побудовано найімовірніші просторові 
структури 1:2 комплексів, які добре узгоджуються з рентгеноструктурними даними [28-30]. 


ИССЛЕДОВАНИЕ СИКВЕНС-СПЕЦИФИЧНОСТИ СВЯЗЫВАНИЯ ДАУНОМИЦИНА 
С ДЕЗОКСИГЕКСАНУКЛЕОТИДАМИ МЕТОДОМ ЯМР-СПЕКТРОСКОПИИ 


А.Н. ВЕСЕЛКОВ, В.И. ПАХОМОВ, Л.Н. ДЫМАНТ, Р.ДЖ. HTOH, Д.Б. ДЭВИС 


Изучено комплексообразование антрациклинового антибиотика дауномицина (DAU) с camo- 
комплементарными дезоксигексануклеотидами S’-d(CGTACG) и 5’-d(CGCGCG) различной последо- 
вательности оснований в цепи в водно-солевом растворе методом ЯМР-спектроскопии (500 МГц). 
Двухмерные гомоядерная 'Н-ЯМР-спектроскопия (TOCSY и NOESY) и гетероядерная ‘H-"P 
ЯМР-спектроскопия (HMBC) использовались для полного отнесения резонансных сигналов необме- 
нивающихся протонов и ядер фосфора, соответственно, а также для качественного определения мест 
преимущественной посадки антибиотика при связывании с гексамером. Установлено, что DAU интер- 
калирует в терминальные 4(СрС)-сайть гексануклеотидов, при этом аминосахарная группа антибиоти- 
ка располагается в малой бороздке гексамерного дуплекса, частично перекрывая третью пару 
оснований. 

Измерены концентрационные и температурные зависимости протонных химических сдвигов 
молекул, на основании которых рассчитаны равновесные константы реакций комплексообразования 
DAU сдезоксигексануклеотидами, предельные значения химических сдвигов протонов дауномицина и 
термодинамические параметры (энтальпия AH и энтропия AS) для различных комплексов антибио- 
тик-ДНК (1:1, 1:2, 2:1 и 2:2) в водном растворе. Установлено, что имеет место определенная специфич- 
ность связывания антибиотика с триплетной нуклеотидной последовательностью, содержащей две 
смежные С-С-парь оснований, фланкированные А.Т-парой оснований в гексамере 5’-d(CGTACG), по 
сравнению с триплетом, состоящим из трех последовательных С.Сб-пар оснований в гексамере 
5'-4(СССССС). Анализ полученных результатов показал, что специфичность связывания дауномици- 
на с ДНК обусловлена гидрофобными взаимодействиями аминосахарной группы антибиотика в малой 
бороздке двойной спирали. На основе расчетных значений индуцированных протонных химических 
сдвигов антибиотика в составе интеркалированных комплексов и данных 2M-NOE построены наиболее 
вероятные пространственные структуры 1:2 комплексов, которые хорошо согласуются с рентгеност- 
руктурными данными [28-30]. 
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Abstract. Protein folding is an intriguing problem for both experimental and theoreti- 
cal studies. Protein molecules are synthesized as linear polymer chains. They adopt 
their three-dimensional structures by a stochastic process of self-assembly. The infor- 
mation on the protein three-dimensional structure is imprinted in a sequence com- 
posed of 20 amino acids. This “folding code” is degenerate and distributed over the 
whole sequence. In the present report we overview the two basic approaches to pro- 
tein folding problem — “pathway view” and “landscape view”. Next we present the 
analysis of protein folding problem based Boltzmann distribution function and show 
that considering the folded state being in equilibrium with the magnitude of unfolded 
and partially folded states results in complete disagreement with experiment, a para- 
dox. Any motion along any trajectory from native to denatured state will not be re- 
versed, it will be lost among the trajectories connecting astronomical number of 
denatured states. The ideas to resolve this paradox are discussed. 


Keywords: protein, native states, folding, activation energy, transition states 


Introduction 


Proteins are unique macromolecules. They are synthesized as linear chains of 
amino acid residues and attain their three-dimensional structure as a spontane- 
ous process of folding. Folding is a process of realization of linear information 
contained in a sequence of made up by 20 different amino acids and its transfor- 
mation into a three-dimensional pattern. This process is spontaneous, it does not 
require neither supply of energy or information from the outside [1,2]. Folded 
protein structures are compact, highly ordered, saturated by intramolecular 
bonds (non-covalent electrostatic and hydrophobic interactions, hydrogen 
bonds). The folded state is separated from unfolded (denatured) state by a high 
energy barrier, which is responsible for maintaining protein structure and func- 
tion ina relatively broad range of external conditions. Protein folds because of ех» 
istence of free energy gradient between unfolded and folded states. Even under 
optimal conditions for the folded state, this state is only marginally more stable 
than the fully unfolded state [3]. Typical values of АС”, for small natural proteins 
are from -5 to -10 kcal/mol (-20 to -40 kJ/mol, roughly by 10-20 КТ).ТБіз small 
gradient is a result of compensation between much larger values, which deter- 
mine enthalpy and entropy of unfolded state largely as the interactions with the 
solvent and that of folded state — with neighboring groups of atoms in protein 
structure [4,5]. | 

The ability of proteins fold into unique three-dimensional structures sponta- 
neously, without external supply of information or energy, has attracted the 
minds of many researchers. The information about these structures are coded in 
amino acid sequences [6], but our abilities to predict the structures based on the 
sequences and to design new proteins which could adopt desired structures re- 
main very limited. Most obscure are the initial steps of folding, which because of 
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shorter duration and smaller effects produced are hidden under the transitions 
occurring at the latest steps. Does the folding occur as a random search between 
the huge number of possible conformations, or this search is directed and facili- 
tated by additional mechanisms? Because a variety of configurations exists for ev- 
ery amino acid residue in the chain, the whole chain can in principle adopt an as- 
tronomical number of configurations and the random search is absolutely incom- 
patible with real, often very fast, folding rates. 


How the protein folds? 


Presently two different methodologies are explored in protein folding re- 
search. One of them is based on the classical point of view on protein folding as a 
reversible process of a transition between two states, unfolded (U) and folded 
(F), the folding can be viewed as an ordinary reversible chemical reaction [7,8]: 


U eF (1) 


Despite its simplicity, this equation allows to describe a number of important 
features of folding/unfolding reaction observed in experiment. In particular, 
these are the ability to observe two U and F forms in equilibrium [2,7,8], exponen- 
tial rate function of acquisition this equilibrium [9], the Arrhenius temperature 
dependence of this reaction [10]. The process of folding often exhibits a more 
complicated behavior, which is usually described by the presence of reaction in- 
termediates [11-13] on the folding/unfolding pathways: 


Us..cloF (2) 


The major transition state is thought to be between the latest intermediate I 
and the F state. Sometimes the folding intermediates are observed directly by 
physical methods [14,15], and it is a matter of discussion if they belong to 
on-pathway or off-pathway routes [16]. Regarding the activation barrier, it is usu- 
ally thought that this barrier represents a highly ordered protein chain configura- 
tion very close to the folded state [6], which still lacks the necessary non-covalent 
bonds to stabilize this structure, and therefore this barrier contains both enthalpy 
and entropy contributions [4,17]. 

The other methodology, which is sometimes called the “new view” [19,20] 
has been developed recently based on computer modeling on protein latex mod- 
els. It eliminates the chemical reaction formalism and the pathway concepts [20]. 
Instead of considering in the course of folding of discrete U and F states, it consid- 
ers a landscape made up of enormous population of microscopically different 
sub-states and attempts to describe them by statistical mechanics models. The 
observable macroscopic conformational “states” are analyzed as the distribu- 
tions or ensembles of individual chain conformations. The U state represented by 
tremendous number of polypeptide chain conformations can funnel to a single 
folded state by diffusing along multiple routes in conformational space. The over- 
all reaction rate is an integral over many simultaneously occurring smaller pro- 
cesses. Within this concept, the reaction intermediates are just kinetic traps, 
which are not necessary to be followed by each member by the unfolded state 
distribution. | 
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Presently the classical pathway view and the landscape view are not only dif- 
ferent concepts, they are applied to different classes of problems. The classical 
approach in combination with solvent modulation of protein stability [21] and 
site-directed mutagenesis of protein sequences [22] became an extremely impor- 
tant tool in elucidation the sequential steps of formation of native protein struc- 
ture [8]. The landscape view in contrast, is intensively explored in computer mod- 
eling of stochastic processes in folding using the simplified models [23,24]. The di- 
versity of unfolded protein conformations, impossibility to fold a protein in one 
or several sequential steps and the necessity of parallel folding pathways were rec- 
ognized long ago [7], but there was no demand to develop these issues form the 
side of experiment. The experiments on protein folding usually start by rapid 
transferring the system to the folding permissive conditions (the conditions 
where the folded state is thermodynamically more stable). This reaction usually 
exhibits exponential time dependence, in amazing agreement within the descrip- 
tion of an activated chemical reaction, in which the equilibrium is shifted towards 
the product state. 


Why proteins fold so fast? Levinthal paradox 


The assumption оп а single pathway between U and Е states, however, results 
in a significant conceptual problem known as Levinthal [25] paradox. The para- 
dox is in the following. Since the unfolded state is a state composed of an astro- 
nomical number of conformations of polypeptide chain, it can not find the ther- 
modynamically stable native state without a random search between these con- 
formations, which requires astronomical times. Why then protein folds so fast? 

The attempts to resolve Levinthal paradox have stimulated the development 
of a “new view”. Instead of finding a transition from U to F state as a random 
search between astronomical number of conformations the system can converge 
to F state by many pathways very fast under the influence of permanent (down- 
hill) gradient of free energy [11]. The smooth (low frustrated [23] ) free energy 
gradient connecting the U and F states is needed for that. The process of folding is 
viewed as a motion on a multidimensional landscape of protein conformations, 
which adopts a shape of a funnel oriented downhill along the free energy gradi- 
ent. The motion along this gradient is associated with the formation of native 
bonds and the decrease of the number of conformations. This motion allows 
maintain favorable entropy-enthalpy balance, and the entropy decrease along 
the pathway, which is associated with the decrease of the number of chain confor- 
mations, is compensated by enthalpy changes [26]. Levinthal paradox was usually 
considered for a case of folding of unfolded nascent protein under folding permis- 
sive conditions, when the folded state is more stable, than an unfolded one. 


A new populational paradox 


There are conditions, however, when the folding state has no energetic pref- 
erence, and both folded and unfolded states exist together in rapidly established 
dynamic equilibrium (Fig. 1). So, how does the reaction of folding occur without 
downhill energy gradient? And ifit exists, is it sufficient to drive the folding path- 
way? We investigated this situation using a computer model based on Boltzman 
distribution of protein states at equilibrium. 
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U<=F 


РЕНН eee E EAN Pe EEEE 
Fig. 1 Schematic presentation of F and U states in configurational space. How a unique F state and an astro- 
nomical number of U states co-exist in equilibrium? 


Consider a system, which can exist inn states of equal or almost equal energy. 
Thermally activated transitions between the states occur easily. Let us select one 
state, j, between 1 andn and find the probability of finding the system in the state j 
after thermal equilibration. This probability, Р, will be determined by Boltzmann 
probability function [27,28]: 


Fee (3) 


ізі 


where E, is the free energy of conformation j, k is the Boltzmann constant and the 
summation is performed over all и states. Eq. (3) is a very general model- 
independent expression, which does not require any knowledge about particular 
system, and the only requirement is the thermodynamic equilibrium. Therefore any 
additive measurable parameter characterizing the whole population of states in the 
system (e.g. protein conformations) can be computed by summation of contribu- 
tions from individual states with the account of their Boltzmann probabilities. 

With no additional assumptions on the behavior of the system we can add one 
more state, М, which exists in equilibrium with all other п states. The state М will 
serve for folded (native) state in protein folding, and the other n states will repre- 
sent the states, which may exist in equilibrium with the native state. The probabil- 
ity of populating the N state will be: 


-Е у /КТ 
e^" 


е-ЕкТ + $ e ЕТ | (4) 
i=l 


Py = 


Using Eq. (4) we can study the behavior of our system with the increase of the 
number of states, п. 
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Fig. 2 Illustration of the model. The state N exists in thermodynamic equilibrium with п states. The states п 
are considered to be of equal energy, while the energy of state N is variable by the value Е. 


Let us first suggest that allistates, including the state М, are of the same stabil- 
ity: E,=E,. Then, P, will decay as a function of п as п, and already at n>100 it 
will become below 1%. This determines the limit on the size of the system in 
which the population of state N can be observed at a detectable level. 

Consider then the system with the increased stability of state N, which should 
possess lower energy well (Fig. 2). For convenience, this increase of stability is ex- 
pressed in КТ units: AE,=E-E, = акт. Using Eq. (4) we can analyze in detail the 
Boltzmann probability P, as a function of two parameters — AE, and n. 

Primarily, we ask the question, how P, will change with the increase of n at 


different AE, values (Fig. 3). With the energy difference of 1-2 kT, P, (n) decays 


very fast, with no essential difference with the case of AE v=0. The population of 
state N falls below 1% already at several decades of the change of n. With the in- 
crease of energy gap, the probability is concentrated іп N state in a larger extent, 
and the decrease of P, is observed at substantially larger n values. It is essential, 


to note, that with the increase of AE, up to 20 KT (which may be considered as the 
higher limit of stability of small globular proteins) Р, still sharply decreases on the 
scale of n of ten decades. To make it infinitely small no astronomical number of 
states is needed. 

Representative in this respect are the results obtained by assigning a particu- 
lar state probability P, (from 0.01 to 0.99) and finding ЛЕ, values, which are 
needed to achieve this probability at particular number of states (Fig. 4). These 
functions are linear if n is presented on logarithmic scale. To achieve a 99% prob- 
ability of N state, in a two-state equilibrium AE, of only 5 КТ is needed, with 10 
states it is about 7 KT, while with 107 states the excess stabilization energy of 18 КТ 
is required. It is clearly seen that with infinitely largen, AE, will also tend to infin- 
ity as a function AE, /kT = In(Pyn/(1-P, )). 

These calculations are made for the state of thermodynamic equilibrium. 
What happens on the approach to equilibrium? Any pathway from the state N 
leads to one ofn states, while any transition from any оби states to № state becomes 
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Fig. 3 Calculated values of Boltzmann state probability, Py, as a function of a number of states, п, at different 
AE, values. 
A — small number of states (п < 10), linear scale. B — large number of states (n < 10), logarithmic 
scale. AE, =0(1); АТ (2); 2kT (3); 5kT (4); 10kT (5); 15 kT (6) and 20 kT (7). 


extremely low probable because of the necessity of choice within a huge variety of 
n states. If the state N is populated initially, its population will vanish in time as a 
result of attaining the equilibrium. Thus, we easily derive a statement that a highly 
populated single state in the system can not exist in dynamic equilibrium with a 
relatively large population of other states of similar energy. 

The transitions between folded and unfolded protein states are highly coop- 
erative, they occur іп a relatively narrow range of external conditions, such as low 
pH, high temperature, high pressure or presence of chemical denaturant [29]. In 
the absence of off-pathway reactions, such as intermolecular aggregation of un- 
folded protein, they display a high level of reversibility [1,7,8,29]. Consider any 
transition of protein conformation between F and U states occurring under the 
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n 


Fig. 4 Calculated values of AE, as a function of a number of States, n, (logarithmic scale) at different 
Boltzmann state probabilities for the М state, Py. 
Py = 0.01 (1); 01(2); 05(3); 09(4); 099 (5). 


influence of changing the external conditions. Let the reaction be completely re- 
versible. If the folding is a two-state co-operative process, then in the transition 
range the folded and unfolded forms should exist in equilibrium [30]. In this case 
the equilibrium constant K, and free energy of folding ЛО", (the difference be- 


tween the free energies of F and U states, G^, and С.) can be evaluated from the 
equations: 


Ка =[Е [О] (5) 
АС° = С, -G°u = RTIn K a» (6) 


Within the transition range we can always find a point where К, =1, and the two 
forms are represented in equal concentrations (the half-transition point). This situa- 
tion is depicted in Fig. 1, which illustrates that the U and F forms are quite different 
physical entities. The F state is characterized by a very narrow distribution in multi- 
dimensional conformational space, while for U state the conformational space in 
which this state is distributed is of astronomical dimension. 

It is believed [30], that in the reaction of folding the transition state is close to 
the folded F state and the activational barrier contains both enthalpic and 
entropic components. It is a low-entropy conformation, which is almost folded 
but lacks necessary bonds to stabilize it. Can this state be easily achieved by diffu- 
sion from any part of a distribution over astronomical number of possible 
polypeptide chain configurations, which is the U state? (Fig. 5). 

Simple estimates demonstrate that in the absence of thermodynamic driving 
force this is not possible. Consider any trajectory leading from the F state over 
transition state to the U state. The unfolded product may adopt one of the 
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astronomical number of configurations. It is separated from F state by a relatively 
high activational barrier, while the barriers separating it from other U states are 
lower and the choice immense. The free energy gradient which could drive the 
system to the F state does not exist (we set K,, = land AG°r =0, see Eq. 6), so this 
product will be trapped among the configurations of the U state without any 
chances to return back to the F state. 

The small ДО", values of -5 to -10 kcal/mol result in the estimates of K,, val- 


ues (Eq. 5) ona level of 10"-10". This means, that even most folded proteins under 
the conditions of their optimal stability must exhibit spontaneous unfold- 
ing-folding transformations and spend some time in the unfolded state. This fact 
has been recognized long ago based on the finite rates of hydrogen-deuterium ex- 
change reaction exhibited by native proteins [31] and confirmed by numerous ex- 
perimental data [32,33]. In order to maintain the concentration of the folded state 
high, the rates of folding should exceed the rates of unfolding by many orders of 
magnitude. Thus, if the rate of total unfolding is 107-107" s, which corresponds to 
half-times from 2 h to 80 days, the rate of spontaneous refolding should be as high 
as 1 s” [30]. Experimental data [8] demonstrate the high rates of refolding after 
the spontaneous unfolding in the native state. Since the spontaneous unfoldings 
do occur, they should depopulate occasionally the folded state and populate one 
of the unfolded states. The way back can not easily be found (Fig. 5) and the nec- 
essary balance between unfolding and folding rates can not be maintained. If this 
can happen, the protein folded state can not exist at all. 

Let us then go back to Levinthal paradox, which 15. considered for unfolded 
protein under folding conditions. Taking typical values of АС’, from -5 to 
-10 kcal/mol (-20 to -40 kJ/mol, roughly by 10-20 КТ) and dividing by a number 
of amino acid residues in a polypeptide chain and the number of configurational 
degrees of freedom for every amino acid, we obtain the infinitesimally small values 
of 0.01—0.02 kcal/mol (assuming a protein with 200 residues with 3 degrees of free- 
dom per residue). They are even much lower than kT (~0.6 kcal/mol). How with 
this small energy difference between neighboring chain configurations one of them 
will become depopulated and the other — populated on the pathway of folding? 
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The nature of unfolded state 


Thus, we demonstrate that there can be no equilibrium between unique 
well-organized folded state and an astronomical number of unfolded states. If a 
molecule appears in the unfolded state, it gets an immense number of choices and 
no guide for its proper folding. The motion of this molecule will be lost among the 
trajectories connecting astronomical number of denatured states. Thermody- 
namic driving force is not sufficient to make the process of folding uni-directional 
and preserve the folded structures from spontaneous unfolding. This statement is 
in complete contradiction with а common practice and represents a new paradox. 

The clue to this paradox may also lie in specific properties of the unfolded 
state. What makes astronomical number of states to behave in the folding reac- 
tion like a single well-defined state? Why a formal description of folding/unfold- 
ing based on simple schemes of chemical kinetics does not originate a paradox? 

Under unfolding conditions, where the interactions between protein groups 
are weaker than the interactions with the solvent the unfolded state can be inter- 
preted as a thermodynamical phase — an ensemble of interconverting conforma- 
tions. Whether we attribute the term “state” to a whole system or to all members 
of the distribution individually, is of no principal importance, there is no possibil- 
ity to distinguish a particular member of the distribution or characterize its prop- 
erties. Even at strong unfolding conditions the polypeptide chain does not oc- 
cupy an astronomical number of conformations, there exists some forming and 
breaking transient structure [34,35]. When the system is transferred to folding 
conditions, not only a polypeptide chain as a whole, but individual groups of at- 
oms become unstable in their interaction with the solvent and tend to associate. 
This reduces substantially а number of populated conformations. It is a step along 
the folding pathway, but this step is neither specific nor cooperative. This behav- 
ior is analogous to the well-known contraction of polymer chains when moved 
into a poorer solvent [36]. The system is not at thermodynamic equilibrium, thus 
the Boltzmann distribution is not established. Then on the folding pathway the 
system converges to a state, which is much closer to the native state, since it con- 
tains a variety of native transient contacts (the pre-folded state, P [7]). The ob- 
served cooperative folding transition occurs essentially between pre-folded and 
folded states with relatively high activational barrier. This explains the frequently 
observed exponential folding kinetics and other properties, which do not require 
for explanation anything beyond Eqs. (1) and (2). The observations of dispersive 
folding kinetics are rare [37]. They may result from a variety of folding pathways 
when the activational barrier is unusually low. The state of U <> F equilibrium is 
actually the P <> F equilibrium with a low population of more unfolded states, 
which appear at higher energies. 

The intermediate states, described by Eq. (2) are essential parts of folding 
mechanisms of many proteins. In the “pathway view” formulation they decrease 
the height of an activational barrier, while in “landscape view” they allow to 
achieve the less frustrated folding route. In the folding kinetics they can appear 
transiently, while at equilibrium they can be populated according to Boltzmann 
distribution. How many of them can exist? Our estimates (see Fig. 2) show, that 
several millions, not more. But this is not an astronomical number, as it could be 
intuitively expected based on “landscape view”. 
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БІЛКОВИЙ ФОЛДІНГ AK ПРОБЛЕМА ТЕОРЕТИЧНОЇ БІОФІЗИКИ 
О.П. ДЕМЧЕНКО 


Білковий фолдінг — це інтригуюча проблема як для експериментальної, так і для теоретичної 
біофізики. Білкові молекули синтезуються як лінійні полімерні ланцюги. Вони набувають тривимірної 
конформації у стохастичному процесі самозборки. Інформація про тривимірну структуру білка закла- 
дена в послідовності, що складається з 20 амінокислот. Цей "код фолдінгу" вироджений і розповсюдже- 
ний на всю послідовність. 

В цій роботі ми розглядаємо два основні підходи до проблеми білкового фолдінгу, що були обра- 
зно названі "шлях" і “ландшафт”. Далі ми представимо аналіз проблеми білкового фолдінгу, що базуєть- 
ся на функції розподілу Больтсмана і покажемо, що якщо вважати, що нативний стан має бути в 
рівновазі з величезною кількістю денатурованих і проміжних станів, то такий стан не зможе існувати, а 
це являє собою парадокс. Це тому, що любий рух за траєкторією від нативного стану до астрономічної 
кількості денатурованих станів не зможе бути пройдений у зворотному напрямку, а загубиться проміж 
траєкторій, що з'єднують денатуровані стани. Обговорюються ідеї розв'язання цього парадоксу. 


ФОЛДИНГ БЕЛКА КАК ПРОБЛЕМА ТЕОРЕТИЧЕСКОЙ БИОФИЗИКИ 
А.П. ДЕМЧЕНКО 


Фолдинг белка — интригующая проблема как для зкспериментальньхх, так и для теоретических 
исследований. Белковые молекулы синтезируются как линейные полимерные цепочки. Они образуют 
трехмерную структуру в стохастическом процессе самосборки. Информация о трехмерной структуре 
белка заложена в последовательность, состоящую из набора из 20 аминокислот. Этот “код фолдинга” 
вырожден и распределен по всей последовательности. 

В настоящем сообщении мы рассмотрим два основных подхода к проблеме белкового фолдинга 
— “направленный путь” и “ландшафт”. Далее мы представляем анализ проблемы белкового фолдинга, 
основанный на функции распределения Больтсманна, и показываем, что рассмотрение свернутого CO- 
стояния в равновеии с большим количеством развернутых и частично свернутых состояний приводит к 
полному несоответствию эксперименту, парадоксу. Любое движение по любой траектории от нативно- 
го к денатурирванному состоянию не будет обратимым, оно будет потеряно между траекториями, свя- 
зывающими астрономическое число денатурированных состояний. Обсуждаются идеи разрешения 
этого вопроса. 
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Abstract. The influences of Zn, Cd, Cu, Ag, Au, Fe cyanide complexes on the 
electrokinetic and transmembrane potentials of cells of Pseudomonas fluorescens 
strain B-5040 has been studied and the correlation between the values under investiga- 
tion is found. It is investigated that degree of damage in a membrane by the cyanide 
complexes depends on the coordinating metal nature and is determined by the bio- 
chemical transformation rate. 


Ключевые слова: электрокинетический потенциал, трансмембранный 
потенциал, микробные клетки 


1. Введение 


Токсичность цианидов для живых организмов хорошо известна и свя- 
зывается с их повышенной способностью образовывать комплексы с метал- 
лами, входящими в состав клеточных ферментов [1]. Тем не менее некото- 
рые организмы способны не только выживать в присутствии цианида, но и 
выделяют его в процессе жизнедеятельности. Это уникальное и полезное 
свойство микроорганизмов существовать в окружении высокотоксичных 
ионов и одновременно их утилизировать используется при разработке эко- 
логически безопасных методов обезвреживания цианидов с помощью 
микроорганизмов. 

Предполагается несколько механизмов устойчивости клеток к циани- 
дам, включающих в том числе их биохимическую трансформацию в неток- 
сичные вещества с помощью ферментов (оксигеназ, цианаз и др.) [2], а так- 
же наличие у клеток альтернативных систем переноса электронов взамен 
поврежденных цитохромов [3]. Цитохромоксидаза, которая играет ключе- 
вую роль в клеточной энергетике, и цианидоксигеназа, которая катализи- 
рует окисление цианида в цианат, являются мембраннно-ассоциированны- 
ми ферментами. Поэтому при поступлении цианида в клетку и его возмож- 
ном каталитическом окислении могут изменяться как энергетика 
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клеточной мембраны, так и свойства двойного электрического слоя за счет 
изменения поверхностного заряда. Задача становится еще более интерес- 
ной и неоднозначной, если клетка находится в контакте не с простым циа- 
нидом, а с цианидным комплексом металла, который в процессе трансфор- 
мации меняет заряд и высвобождающийся металл также каким-то образом 
взаимодействует с клеточной стенкой и мембраной. 

К настоящему времени уже известны работы, в которых была экспери- 
ментально обнаружена взаимосвязь мембранного и поверхностного потен- 
циалов у различных микрообъектов [4]. Авторы [5] в результате теоретиче- 
ского анализа пришли к заключению, что взаимная зависимость указанных 
параметров определяется, во-первых, перераспределением трансмембран- 
ного потенциала между емкостями мембраны и двойного электрического 
слоя, а, во-вторых, изменением заряда клеточной поверхности. В результа- 
те таких изменений достигается стабилизация внутримембранного поля. 

Целью настоящей работы было исследование методами микроэлектро- 
фореза и проникающих ионов (ТФФ*) электроповерхностных и мембран- 
ных свойств культуры Pseudomonas fluorescens ВКПМ В-5040 в присутст- 
вии простого цианида натрия и комплексов с переходными металлами, что 
позволит оценить биоэнергетические процессы, протекающие в клетке при 
метаболических превращениях цианидов. 


2. Материалы и методы 


В опытах использовали штамм Pseudomonas fluorescens BKIIM B-5040 
из коллекции ИБКХ НАНУ. Культуру выращивали методом глубинного 
культивирования в питательной среде 5М [6]. Для измерения электрофоре- 
тической подвижности клетки дважды отмывали от питательной среды 
дистиллированной водой и концентрировали в фосфатно-цитратном бу- 
ферном растворе (pH 7.6). 

Электрофоретическую подвижность клеток измеряли методом микро- 
электрофореза в прямоугольной ячейке закрытого типа в ручном режиме 
регистрации клеток [7]. Измерительная ячейка была снабжена электродны- 
ми отделениями. Величину электрокинетического потенциала (ЭКП) pac- 
считывали по формуле Смолуховского [8]. 

Трансмембранный потенциал определяли методом электрофореза про- 
никающих ионов тетрафенилфосфоний бромида (ТФФ*). Ero концентра- 
цию измеряли с помощью ТФФ”-селективного электрода [9]. 

Поскольку основные системы трансформации энергии локализованы 
на внутренней мембране, то измерение трансмембранного потенциала у 
грамотрицательных бактерий с помощью проникающих ионов возможно 
лишь после пермеабилизации внешней мембраны трилоном Би трис-буфе- 
ром. Подготовка клеток P.fluorescens к измерению мембранного потенциа- 
ла состояла в следующем: путем центрифугирования при 5600xg клетки от- 
мывали 25 MM трис-НСІ-буферным раствором; содержащим 1% глюкозы, 
РН 7.6. Осадок затем ресуспендировали в 4 мл того же буферного раствора. 
В полученный объем добавляли 800 мкл 50 MM раствора этилендиа- 
минтетраацетата, доведенного с помощью 0.1н KOH до pH 7.6 и перемеши- 
вали 20 мин на магнитной мешалке при температуре 28°С. После этого 
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Рис. 1. Влияние цианида на поглощение ионов ТФФ* клетками Pseudomonas fluorescens B-5040 в pac- 
творах солей : 1 — NaCN, 2 — Na,Cd(CN),, 3 — Na,Zn(CN),, 4 — NaAg(CN),, 
5 — Na,Cu(CN),, 6 — NaAu(CN),, 7 — K;Fe(CN),. 


бактерии осаждали центрифугированием и снова ресуспендировали в 2 мл 
трис-НСі-буферного раствора с глюкозой, присутствие которой в буфер- 
ном растворе способствовало сохранению активности бактерий. Аликвоты 
суспензии подготовленных таким образом клеток вносили в модельные 
растворы, содержащие МаСМ или соответствующую комплексную соль, 
выдерживали не менее 10 мин и затем измеряли интенсивность поглоще- 
ния ТФФ*. Объем измерительной ячейки составлял 5 мл. Все эксперимен- 
ты проводились сериями на клетках одного урожая так, чтобы наблюдае- 
мый результат влияния агентов всегда можно было соотнести с 
контрольным. 


3. Результаты и обсуждение 


Ранее нами было показано, что культура P.fluorescens В-5040 эффектив- 
но разрушает простой цианид натрия, а также его комплексы с цинком, ме- 
дью, серебром и никелем до нетоксичных углекислого газа и аммиака [6]. 
Скорость и эффективность микробного разрушения цианидных комплек- 
сов металлов по сравнению со свободными цианидами заметно зависела от 
термодинамической и особенно кинетической прочности комплексов и 
природы координирующего металла [10]. 

Как видно из представленных на рис. 1 данных, цианиды способны из- 
менять мембранный потенциал изученной культуры. Состав комплекса 
оказывал заметное влияние на интенсивность поглощения ТФФ" клетками 
бактерий, преинкубированных с соответствующим цианидом. 
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Эффективность поглощения оценивали по величине отношения m/m,, где 
m, — концентрация TOO" в среде инкубации с необработанными клетка- 
ми, а т — концентрация в среде инкубации с клетками, предварительно об- 
работанными цианистой солью в течение 10 мин. Благодаря высокой про- 
никающей способности ТФФ" происходило выравнивание его электрохи- 
мических потенциалов по обе стороны клеточной мембраны. Мембранная 
разность потенциалов связана с концентрацией ТФФ* уравнением Нерн- 
ста: AP--RT/F(In (ТФФ'|, / [ТФФ\]), где индексу “i” соответствует KOH- 
центрация ТФФ` во внутреннем объеме клетки, индексу “0” — концентра- 
ция в среде инкубации. По изменению концентрации TOO" в среде инку- 
бации можно было судить об относительном изменении мембранного 
потенциала клетки. Поэтому мы полагали, что отношение m/m,, коррели- 
руя с величиной AY, отражало степень повреждения клеточной мембраны 
цианидами. Использование относительной величины m/m, позволило так- 
же вычесть пассивную составляющую, связанную с простой сорбцией 
ТФФ* полиэлектролитами (белки, полисахариды и др.), заполняющими 
периплазматическое пространство клетки. 

По характеру воздействия на приведенные на рис. 1 зависимости m/m, 
OT концентрации, цианиды можно условно разделить Ha три группы. B nep- 
вую вошли комплексы золота, серебра и меди. Их объединяет немонотон- 
ный характер кривых с максимумом при концентрациях 10-20 мг/л CN’. 
При этих концентрациях, а для дицианаурата при концентрации больше 
25 мг/л CN , обработанные клетки поглощали больше ТФФ*, чем необрабо- 
танные. Указанные комплексы относятся к числу достаточно стабильных, а 
координирующий их металл к группе так называемых электроположитель- 
ных металлов. Во вторую группу вошли комплексы кадмия и цинка. Здесь с 
ростом концентрации комплекса кривые характеризовались вначале плав- 
ным снижением сорбции TOO", а затем стабилизацией на уровне 30-40% 
от уровня необработанных клеток. Остаточный уровень энергизации со- 
хранялся даже в присутствии очень токсичного комплекса кадмия при кон- 
центрации 200 мг/л, хотя комплексы цинка и кадмия относительно непроч- 
ные и легко диссоциируют, высвобождая СМ и ион металла. К третьей 
группе можно отнести комплекс с железом. Клетки отреагировали на вне- 
сение гексацианоферрат-ионов небольшим снижением сорбционной емко- 
сти TOO", дальнейшее повышении концентрации этой соли не отразилось 
существенно на степени энергизации клетки. Данное соединение отличает- 
ся большой кинетической инертностью и остается в растворе в практиче- 
ски недиссоциированном состоянии. 

Известно, что основными генераторами трансмембранного потенциала 
клетки являются АТФаза и дыхательная цель, поэтому можна было ожи- 
дать, что максимальное ингибирующее действие будут оказывать соедине- 
ния, содержащие большую концентрацию свободного СМ”. Однако, как 
видно из рис. 1, преинкубация клеток Р. fluorescens B-5040 с NaCN сопрово- 
ждалась плавным повышением m/m, в изученном интервале концентраций 
10-200 мг/л, а наибольшее повреждающее воздействие на мембраны оказы- 
вали легкодиссоциированные непрочные цианокомплексы кадмия и цин- 
ка. Отсюда можно предположить, что снижение степени энергизации кле- 
ток, находящихся в контакте с цианид-ионами и их комплексами с металла- 
ми, определяется в первую очередь природой металла, входящего в состав 
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Рис. 2. Зависимость электрокинетического потенциала ($) клеток Pseudomonas fluorescens B-5040 от 
концентрации цианида в растворах : 1 — NaCN, 2 — Na,Cu(CN),, 3 — K,Fe(CN),, 
4 — NaAg(CN),, 5 — Na,Cd(CN),. 


комплекса, и обусловлено индивидуальными особенностями повреждения 
мембраны этими металлами. Двухзарядные катионы цинка и кадмия могут 
сорбироваться клеточной стенкой, а возможно проникать внутрь клетки 
путем диффузии, в то время как для комплексов меди и серебра можно ожи- 
дать энергозависимый характер их проникновения внутрь клетки по меха- 
низму антипорта с протонами, образующимися в результате гликолиза 
глюкозы, входящей в состав питательной среды, подобно тому, как это от- 
мечено для аккумуляции ионов меди клетками Escherichia coli в [11]. Дейст- 
вительно, для двух последних комплексов нами было отмечено существен- 
ное закисление инкубационной среды после их введения в суспензию бак- 
терий [6]. Отмеченная на рис. 1 при малых концентрациях активация, при 
более высоких концентрациях сменялась ингибированием, что должно 
привести к остановке транспорта меди и серебра в клетку, ограничивая тем 
самым повреждение мембраны. Характер представленных зависимостей 
свидетельствует в пользу высказанного предположения. 

Являясь мерой метаболической активности клетки, мембранный потен- 
циал определенным образом связан C ее электроповерхностными свойства- 
ми. Исходя из такой предпосылки можно было ожидать, что одновременно 
с изменением мембранного потенциала, контролируемого нами через из- 
менение отношения M/M, под влиянием цианидных соединений различно- 
го состава будут изменяться электроповерхностные свойства клетки, что в 
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Рис. 3. Зависимость электрофоретической подвижности (DOI) клеток от pH в буфферном растворе 
(1), а также в присутствии МаСМ (2), Na,Cu(CN), (3), и NaAg(CN), (4). Концентрация СМ 
20 мг/л. 


свою очередь может повлиять на (-потенциал, поскольку электрокинети- 
ческий потенциал изменяется симбатно поверхностному. 

Изменение электроповерхностных свойств может быть связано как со 
степенью энергизации клетки, так и с процессами биохимической транс- 
формации цианистых соединений, через их взаимодействие с продуктами 
(или полупродуктами) разложения. Приведенные на рис. 2 кривые отража- 
ютизменение ЭКП клеток под влиянием возрастающих концентраций циа- 
нидных соединений. Полученные данные склоняют к тому, что заряд на по- 
верхности активно метаболирующей клетки в цианидсодержащем раство- 
ре зависел не столько от сорбции комплексных анионов, сколько от 
результатов ферментативного окисления цианида. Рост отрицательных 
значений ЭКП клетки с повышением концентрации МаСМ мы связываем с 
образованием в качестве промежуточного продукта разложения цианида 
аспарагиновой кислоты и, следовательно, с возрастанием концентрации 
карбоксильных групп на поверхности клетки [6]. Отметим, что электроки- 
нетический потенциал убитых кипячением клеток не зависел OT концентра- 
ции цианида и был на 7-8 мВ ниже, чем у живых. 

Комплексные анионы Zn(CN); и Са(СМ)?` ввиду ux непрочности дис- 
социируют в растворе, а в слабокислых и нейтральных растворах эти метал- 
лы легко гидролизуют. Поэтому с клеткой взаимодействуют циан-ионы и 
коллоидные гидроксиды указанных металлов, что находит свое выражение 
в небольшом снижении ЭКП. Комплексы Си(СМ) ; M Ag(CN); значитель- 
но более устойчивы и очень слабо диссоциированы в растворе, поэтому они 
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взаимодействуют с клеткой целиком, оказывая на поверхностный заряд 
приблизительно одинаковое воздействие — в присутствии этих соедине- 
ний зависимость носит экстремальный характер. Многозарядный анион 
Fe(CN); изменил ЭКП не более, чем на 5-8 мВ. 

Для интерпретации полученных результатов полезными будут пред- 
ставленные на рис. 3 данные изменения электрофоретической подвижно- 
сти клеток под влиянием рН в растворах, содержащих цианиды различного 
состава. Как видно, характер зависимостей не отличался для необработан- 
ных и обработанных цианидами клеток — закисление раствора вело к сни- 
жению электрофоретической подвижности. Однако наблюдали отличие в 
абсолютных значениях указанных величин. 

Немонотонный характер зависимостей ЭКП от концентрации цианида 
(рис. 2) можно связать с одновременным протеканием нескольких процес- 
сов, вносящих вклад в изменение электрических характеристик поверхно- 
сти клетки. Сюда входят: изменение рН культуральной среды в результате 
накопления во внеклеточном растворе органических кислот и выброса 
протонов в результате гликолиза; изменение числа карбоксильных групп на 
поверхности клетки в результате ферментативной трансформации циани- 
да; взаимодействие клеточной стенки и мембраны с металлом, входящим в 
состав комплекса. 

Сопоставление результатов измерения мембранного потенциала клет- 
ки по поглощению ионов TOO” и ее электрокинетических свойств (рис. 1и 
рис. 2) обнаружило заметную корреляцию между этими величинами, что 
свидетельствовало об их взаимной зависимости. Оба процесса объединены 
процессом метаболического разрушения цианида, в результате которого 
меняется число ионногенных групп на поверхности и происходит ингиби- 
рование мембранных ферментов металлом, освобождающимся при биохи- 
мической деструкции цианида. Этим объясняется тот факт, что клетки, 
проявляющие резистентность к цианиду, остаются к нему достаточно чув- 
ствительными и это прослеживается на мембранных и электроповерхност- 
ных свойствах. Способность утилизировать цианид является по-видимому 
вынужденной защитной мерой, позволяющей клетке существовать в циа- 
нид-содержащей среде. 


ВПЛИВ ЦІАНІДНИХ КОМПЛЕКСІВ ПЕРЕХІДНИХ МЕТАЛІВ НА ЕНЕРГЕТИЧНІ 
ХАРАКТЕРИСТИКИ БАКТЕРІЙ 


В.І. ПОДОЛЬСЬКА, Л.М. ЯКУБЕНКО, 3.Р. УЛЬБЕРГ, Н.І. ГРИЩЕНКО, А.І. ЄРМОЛЕНКО 


Вивчено вплив ціанідних комплексів Zn, Cd, Cu, Ag, Au, Fe на електрокінетичний та трансмемб- 
ранний потенціали клітин бактерій Pseudomonas fluorescens B-5040 і помічена кореляція поміж вивче- 
ними характеристиками. Встановлено, що ступінь ушкодження мембрани ціанідними комплексами 
залежить від природи координуючого метала і визначається швидкістю біохімічної трансформації ком- 
плексу. 


ВЛИЯНИЕ ЦИАНИДНЫХ КОМПЛЕКСОВ ПЕРЕХОДНЫХ МЕТАЛЛОВ HA 
ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ БАКТЕРИЙ 


В.И. ПОДОЛЬСКАЯ, Л.Н. AKYBEHKO, 3.P. УЛЬБЕРГ, Н.И. ГРИЩЕНКО, А.И. ЕРМОЛЕНКО 


Изучено влияние цианидных комплексов Zn, Cd, Cu, Ag Au, Fe на электрокинетический и 
трансмембранный потенциалы клеток бактерий Pseudomonas fluc e: `епз В-5040 и отмечена корреля- 
ция между изученными характеристиками. Установлено, что степеє: овреждения мембраны цианид- 
ными комплексами зависит от природы координирующего ме т? и определяется скоростью 
биохимической трансформации комплекса. 
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КОНЦЕНТРАЦІЙНИЙ ПОТЕНЦІАЛ НЕРНСТА З 
ВРАХУВАННЯМ ТЕРМОДИФУЗІЙНИХ ЕФЕКТІВ 


СП. САМКО, О.В. ЧАЛИЙ 


Київський національний університет імені Тараса Шевченка 


Національний медичний університет імені О.О. Богомольця 


CONCENTRATION NERNST POTENTIAL WITH 
THERMAL DIFFUSION EFFECTS 


S.P. SAMKO, A.V. CHALY 


Abstract. The influence of thermal diffusion effects on the concentration Nernst po- 
tential is examined. The phenomenologic coefficients characterizing the appearance 
of flux of charged particles under temperature gradient were evaluated from experi- 
mental data. The proposed results are of a great importance with taking into account 
the application of hyper- and hypothermia in modern medical practice (cryotherapy, 
oncology and other areas). 


Ключові слова: потенціал Нернста, термодифузія, гіпотермія, гіпертермія 

Однією з найважливіших властивостей основного структурного еле- 
мента живих систем — клітини — є наявність різниці електричних потен- 
ціалів по обидві сторони цитоплазматичної мембрани. Як відомо (див., 
напр. (1-51), розрізняють: 1) потенціали спокою — потенціали в нормально 
функціонуючій незбудженій клітині, які існують між цитоплазмою ізовніш- 
нім середовищем за рахунок різниці концентрацій іонів, а також їх різних 
коефіцієнтів проникності; 2) потенціали дії, які виникають при збудженні 
клітини і пов'язані із зміною іонних проникностей в порівнянні з їх значен- 
нями в незбудженій клітині. 

Такі ж за своєю природою, як іна біологічних мембранах, різниці потен- 
ціалів будуть виникати, якщо, наприклад, два розчини одного і того ж елек- 
троліту з різними концентраціями будуть розділені деякою штучною мемб- 
раною, проникною лише для іонів одного знаку. 

В даній роботі буде вивчений вплив термодифузійних ефектів на зна- 
чення концентраційного потенціалу Нернста, який є одним із основних 
мембранних потенціалів спокою. Така постановка задачі, в першу чергу, 
маєтеоретичне значення, а також може бути цікавою у випадку застосування 
в медичній практиці (кріомедицині, онкології) методів гіпер- та гіпотермії. 


1. Концентраційний потенціал Нернста 


Рівновага в однокомпонентній двохфазній системі має місце, якщо ви- 
конуються умови рівності температури, тиску і хімічних потенціалів обох 
фаз, тобто T=const, P=const, це const. Будемо вважати різними фазами зов- 
нішню (е) і внутрішню (i) сторони біологічної мембрани. Нехай спочатку 
умови рівності температур Т,=Т, і тисків P, P, виконані. В цьому випадку 
потік іонів через мембрану буде описуватись відомим електродифузійним 
рівнянням Нернста-Планка, яке може бути перетворене до іншого співвід- 
ношення для потоку заряджених частинок — рівняння Теорелла: 
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S". Were 
I Pasa (НУ +ВТш С +2Еф). (1) 


Величина, що стоїть в дужках, називається електрохімічним потенціалом ц. 
Він визначає вільну енергію Гіббса в розрахунку на один моль при наявності 
розчиненої речовини і електричного поля. Хімічний потенціал чистого 
розчинника у, вважається однаковим по обидві сторони мембрани. 
Умовою іонної рівноваги є рівність електрохімічних потенціалів обох 
фаз. Звідси випливає, що 
RTI. C; 


Pm 79.79: = e і (2) 


Для оцінки величини концентраційного потенціалу Нернста візьмемо 
температуру Т=293 К, а валентність 2-1, тобто будем розглядати 
одновалентний іон, яким є, наприклад, іон калію. З врахуванням числових 
значень В. =8.31.10° Дж/кмоль:К i F=9.65-10" Кл/кмоль, а також використо- 
вуючи експериментальні дані (див. [1]) для конкретного об'єкта - - гігантсь- 
кого аксона кальмара, в якому концентрації іонів всередині та зовні аксона 
рівні відповідно [К]=392 мМ/л і [K,]=22 мМ/л, маємо таке значення 
мембранного потенціалу Нернста: 9, - 72.67 мВ. 


2. Термодифузійна зміна концентрацій 


Застосовуючи загальні положення термодинаміки необоротних проце- 
сів для опису термомеханічних явищ (див., наприклад, [5,7-10]), запишемо 
вираз для густини дифузійного потоку розчиненої речовини, тобто кількос- 
ті речовини через одиницю площі за одиницю часу, для однокомпонентної 
двохфазної системи, в якій існують відносно малі градієнти концентрації і 
температури: 


" sean CON Ier 
I- AY vc Ead ; (3) 
\ 


де D — коефіцієнт дифузії, а К, — термодифузійне відношення. 

В рамках термодинаміки необоротних процесів кінетичні коефіцієнти 
типу D i K, є феноменологічними параметрами, з'ясуванням мікроскопіч- 
ного змісту яких займається відповідна кінетична теорія. Наприклад, роз- 
в'язок кінетичного рівняння Больцмана в наближенні Лорентца (модель 
пружнього зіткнення малих і великих сферичних частинок) дає такий вираз 
для термодифузійного відношення: K,=n/2, де n — концентрація малих 
частинок [5,9]. В залежності від знаку К. розчинений компонент дифундує 
або в холодну область (при K,>0), або в нагріту (при K,<0). Як показують 
експерименти [11,12], у випадку розчинів електролітів концентрація 
розчиненого електроліту завжди збільшується в холодній області, тобто ве- 
личина К, є додатньою. 
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Для оцінки числового значення термодифузійного відношення K, CKO- 
ристаємось експериментальними даними по термодифузії в рідинах [13]. 
В [13] приведені експериментальні значення термодифузійних відношень 
A, в інтервалі температур Т 298-313 К для рідин, в яких мольна доля 
першого, дифундуючого, компонента становить 0.09 i 0.5. Термодифузійне 
відношення A, є безрозмірною величиною. 3 виразу для дифузійного 
потоку першого компонента в бінарному розчині при наявності термо- 
дифузії [13] можна показати, що величини K, i А, пов'язані наступним спів- 
відношенням: 


M 
К. = Ca (4) 


де C=C,+C, — молярна концентрація розчину, М, — молярна Maca розчин- 
ника, M= X, M ,+x,M, — середня молярна маса розчину, a X, i x, — відповідно 
мольні долі розчиненої речовини i розчинника. 

Приведені значення змінюються в достатньо широких межах, а саме: 
A,=0.06-1.23. Для одержання числового значення термодифузійного 
відношення K, візьмемо мінімальне значення A,=0.06. Величину молярної 
концентрації розчину С оцінимо через молярну концентрацію чистого роз- 
чинника (будемо вважати, що це вода), тобто С, ЕС, 7996 М/л, i середню 
молярну концентрацію калію, тобто С, -207 мМ/л. Таким чином, молярна 
концентрація розчину C=C,+C, =55. 8 М/л, що за формулою (4) дає 
наступну оцінку для величини термодифузійного відношення: K,=3.348 М/л. 

В біологічних клітинах між зовнішнім і внутрішнім середовищем за ра- 
хунок активного транспорту підтримуються необхідні іонні градієнти. От- 
же, в такій системі явище термодифузії приведе до деякого перерозподілу 
концентрації, а отже, і до зміни існуючих градієнтів. 

Нехай температура змінюється всередині клітини. Припустимо Т,>Т,. 
В цьому випадку за рахунок термодифузії частина іонів калію повинна пе- 
рейти із клітини в зовнішнє середовище, тобто в більш холодну область. В 
стаціонарному стані концентрації іонів по обидві сторони мембрани дещо 
зміняться в порівнянні з своїми початковими значеннями. Вважатимемо, 
що перерозподіл концентрацій по обидві сторони мембрани є симетричним 
(це показують експерименти [14]), тобто концентрації іонів калію зовні і 
всередині клітини зміняться на одну і ту ж величину AC. Тоді С/= С-АС i 
C/=C,+AC, де C/ і C/ — концентрації іонів всередині і зовні клітини, які 
встановилися в стаціонарному стані. Для оцінки величини АС використа- 
ємо умову стаціонарності 


II ee o4 5 
dx T dx 2 
Інтегрування рівняння (5) y лівій частині дасть різницю концентрацій 
між двома фазами з температурами Т; 1Т,, яка виникла за рахунок термоди- 


фузії, тобто 2AC. Вважаючи K, =const при деякій середній температурі, 
після інтегрування рівняння (5) отримуємо для АС наступний вираз: 
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ln ~. 
T (6) 


При чисельому значенні термодифузійного відношення K,—3.348 М/л 
значення величин AC, С; i С," в залежності від різниці температур приведені 
в таблиці 1. 


Таблиця 1 
T, K | AT, K | AC, мМіл | C, , мМ/л C, , мМ/л 
293.01 | 0.01 | 0.057 | 391.943 | 22.057 
2931 | — 01 | 0571 | 391429 22.571 
| 294 | 1 | 57 | 3863 | 277 
296 | 3 | 17 | 375 | 39 
28 | 5 | 28 | зе 50 | 


Для окремої клітини, яка поміщена в деякий фізіологічний розчин, що 
має об'єм набагато більший за об'єм клітини, слід вважати, що зменшення 
концентрації іонів всередині клітини суттєвим чином не вплине на їх кон- 
центрацію зовні і тому С/=С-АС1С/=С.. 


3. Розрахунок концентраційного потенціалу при наявності градієнта 
температури 


Потік заряджених частинок (в нашому випадку - - одновалентних іонів 
калію) через цитоплазматичну мембрану при наявності градієнтів концен- 
трації, потенціалу електричного поля і температури записується наступним 
чином: 


2 Е ome 
I = -D| VC +— VT | - CbVo, 
| T Ф (7) 


де b — рухливість іонів. Формула (7) являє собою узагальнення електродифу- 
зійного рівняння Нернста-Планка при наявності градієнта температури. 

Перетворемо це рівняння до відповідного узагальненого рівняння Тео- 
релла. Очевидно, що рівняння (7) повністю еквівалентне наступному 
рівнянню (8): 


Е = Dis = Јат zl 
ls «se ELE Es -InTV З) (8) 
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Тепер розкладемо функцію Т(х), яка описує розподіл температури 
всередині мембрани (товщина мембрани L), в ряд Тейлора в околі точки 
хе 0, розміщеної посередині мембрани: 


че-то) (A) ee ө 


Тоді в наближенні постійного градієнта температури, тобто при збере- 
женні перших двох доданків у розкладі (9), розподіл температури в мембра- 
ні буде описуватись наступним чином: 


T(x) = 4 +1, xd = x, T(0) = 


e, 


(10) 


Із (10) безпосередньо випливає, що температура Ha зовнішній і 
внутрішній границях мембрани рівна відповідно 


І, dE. L a (ha ae 
(-£]-19 57^ = Te (2-10) ^^ -7. (11) 
Враховуючи, що відношення (T;-T,)/(T,* T.) є набагато менше одиниці, 
для InT(x) отримуєм наступний вираз: 


4. EUM 2(Т, =T,) X 
InT(x) dU xD ur EU (12) 


Чисельні оцінки обох доданків в правій частині формули (12) 
показують, що при температурах Т,=293 К i Т, -298 К чисельне значення 
другого доданку є на три порядки меншим від першого. Це дозволяє з 
високою точністю зберегти в (12) лише перший доданок, а саме: 


InT(x) «In A = const. (13) 


Таким чином, вихідне рівняння (7) для потоку частинок з врахуванням 
виразів (8), (13) і співвідношення Ейнштейна b=D(zF/RT) зводиться до 
узагальненого рівняння Теорелла наступного вигляду: 


= DS ZED zFD FD K 2T 
I - 29 + ZFo + НУ se cole PE. Te) (14) 


де сума перших трьох доданків B фігурних дужках являє собою звичайний 
електрохімічний потенціал у, а останній доданок описує внесок від градієнта 
температури. Величину в фігурних дужках рівняння Теорелла (14) для 
системи з градієнтом температури природньо назвати термоелектрохімічним 
потенціалом |’. 

З умови рівності термоелектрохімічних потенціалів по обидві сторони 
мембрани, тобто и; = и, отримуєм наступний вираз для концентраційного 
потенціалу Нернста з врахуванням термодифузійних ефектів фу: 
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qe ue ME ONE TE: 
gps dur inet c quc e |, 
TEM | у (е ит i] d 


1 1 


Для різниці температур Т, - Т, << T(0) легко показати, що 
27, Т-Т Q 25 -LOT 


In ; - 
Ir ТАС Xm dem 


В результат! концентраційний потенціал Нернста на мембрані з 
врахуванням градієнта температури буде описуватись наступною фор- 


мулою: 
EU EA 
= c In 1 K 1 e t Я 
Е | Ur T G = n2 


1 e 


Tyr T, — середня температура в системі. Отриманий вираз для різниці 
електричних потенціалів складається з двох доданків, де перший член являє 
собою звичайний концентраційний потенціал (з врахуванням перерозподілу 
концентрацій), а другий обумовлений безпосередньо градієнтом температу- 
ри і при Т,- T, перетворюється в нуль. 

Скориставшись чисельними оцінками з 11.1 i 2, а також вважаючи, що 
С; =C, - AC i C; =C,, можна оцінити потенціал за формулою (16). Значен- 
няфи B залежності від різниці температур АТ=Т-Т, приведені в таблиці 2. 


Таблиця 2 
AT, K | Pu, MB | 
0.01 | 72.73 
0.1 | 073.34 | 
І | 79.34 | 
3 | 92.73 
5 | 106.16 | 


4. Коефіцієнт Cope і розрахунок концентраційного потенціалу при 
наявності градієнта температури 


При розгляді термодифузії в розчинах |7,11,12) густину потоку розчине- 
ної речовини записують в наступному вигляді: 


I, = -nD(Vx, +0XX,VT) j (17) 
де n — повне число молів (розчиненої речовини i розчинника) Ha одиницю 


об'єму; х, і х, — відповідно мольні долі розчиненої речовини і розчинника; 
D — коефіцієнт концентраційної дифузії; о — коефіцієнт Соре. 
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Для оцінки величини коефіцієнта Соре і термодифузійної зміни кон- 
центрацій скористаємось експериментальними даними по термодифузії в 
рідинах і розчинах електролітів (14). В [14] приводиться експериментальне 
значення коефіцієнта Соре с для розчину Н,О-КСІ з концентрацією розчи- 
неного електроліту 0.2 М/л при середній температурі 293 K: o= 1.18-10° K”. 
Аналогічно до того, як це було зроблено B п.2, можна оцінити зміну KOH- 
центрацій іонів по обидві сторони мембрани під впливом градієнта темпе- 
ратури. В таблиці 3 приведені значення величини AC для гігантського ак- 
сона кальмара при Т, 293 К. 


Таблиця 3 
Т, К | C, MM/n 
0.01 | 0.002 
| 0.1 | 0.022 | 
| 1 | 0.218 
| 3 0.655 | 
| NM XN _ 1.091 


Тепер запишемо густину потоку заряджених частинок при наявності 
градієнта концентрації, потенціалу електричного поля і температури, в 
наступному вигляді: 


I= -b(vc й CoVT) - CbVo. (18) 


При цьому вважається, що мольна доля розчинника прямує до одиниці, тобто 
у формулі (17) x, — 1 Перетворемо це рівняння до вигляду рівняння 
Теорелла. З цією метою спочатку винесемо множник (-Cb/zF) разом з опера- 
тором градієнта, що дає і 
Е Cb = zFD 
ТЕ оао +> inc SED ar |. 
Z 


(19) 
Величину, що стоїть в дужках, по аналогії з попередніми викладками, назвемо 
термоелектрохімічним потенціалом. 

Умовою іонної рівноваги є рівність термоелектрохімічних потенціалів 
по обидві сторони мембрани. Звідси знаходимо вираз для різниці електрич- 
них потенціалів з врахуванням градієнта температури: 


х RT, C; 
Ou = ET {in с. -+ o(T, = то) (20) 
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Отримана формула (20) при Т,-Т, переходить y вираз (2) для рів- 
новажного концентраційного потенціалу Нернста. Перший член являє со- 
бою звичайний дифузійний потенціал. Другий член можна назвати “термо- 
дифузійним потенціалом". В початковий момент, коли ще відсутній 
перерозподіл концентрації, зумовлений термодифузією, має місце 
термодифузійний потенціал і звичайний рівноважний потенціал Нернста. 
Зміна дифузійного потенціалу виникає поступово і обумовлена перерозпо- 
ділом концентрації іонів по обидві сторони мембрани. Значення потенціалу 
для окремої клітини, розраховані за формулою (20), приведені в таблиці 4. 


Таблиця 4 
| AT, К фи, MB 
| 0.01 72.67 
| 0.1 227 
| _ І | 72.8 
| 3 | 73.1 
| 5 | eo 7734 


Таким чином, чисельні оцінки величини концентраційного потенціалу 
Нернста з врахуванням термодифузійних ефектів, отримані першим 
(табл. 2) і другим способом (табл. 4), виявляються не співпадаючими, але, як 
видно з отриманих результатів, в обох випадках при T,» T, потенціал Bce- 
редині клітини стає більш від'ємний, Відмінність результатів в обох випад- 
ках може пояснюватись, насамперед, тим, що величина К. досить сильно 
залежить від концентрації дифундуючого компонента. При середній 
концентрації іонів C=207 мМ/л їх мольна доля в розчині (розчинником 
вважається вода) становить 0.004. Це є на один-два порядки менше ніж 
мольна доля дифундуючого компонента в експериментах по термодифузії 
[13], з яких оцінювалось термодифузійне відношення К... 

Здається, що перевагу слід віддати результатам, які були отримані через 
оцінки коефіцієнта Соре, а саме, другим способом, але однозначно 
стверджувати цього не можна. Це пов'язано з тим, що експеримент по тер- 
модифузії (14| проводився в умовах, коли градієнт температури становив 
3 град/мм. Коефіцієнт Cope зростає при збільшенні градієнта температури 
(див., напр., [11,12]), а тому при більших градієнтах можна чекати зростаю- 
чого впливу термодифузійних ефектів. 

Подальше дослідження цієї проблеми потребує постановки прецизій- 
них експериментів, які могли б кількісно описати вплив термодифузійних 
ефектів на значення концентраційного потенціалу Нернста. Важливу роль 
при цьому грає експериментальне визначення таких термодифузійних 
характеристик як термодифузійне відношення і коефіцієнт Соре. 
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КОНЦЕНТРАЦІЙНИЙ ПОТЕНЦІАЛ НЕРНСТА З ВРАХУВАННЯМ ТЕРМОДИФУЗІЙНИХ 
ЕФЕКТІВ 


С.П. САМКО, О.В. ЧАЛИЙ 


Досліджено вплив термодифузійних ефектів на значення концентраційного потенціалу Нернс- 
та. Феноменологічні коефіцієнти, що характеризують виникнення потоку заряджених частинок під 
дією градієнта температури, оцінювалися із експериментальних даних. Запропоновані результати є ду- 
же важливими, враховуючи застосування в сучасній медичній практиці (кріотерапії, онкології та ІНШИХ 
областях) методів гіпер- та гіпотермії. 


КОНЦЕНТРАЦИОННЫЙ ПОТЕНЦИАЛ НЕРНСТА С УЧЕТОМ ТЕРМОДИФФУЗИОННЫХ 
ЭФФЕКТОВ 


С.П. САМКО, A.B. ЧАЛЫЙ 


Исследовано влияние термодиффузионных эффектов на значение концентрационного потен- 
циала Нернста. Феноменологические коэффициенты, которые характеризуют появление потока заря- 
женных частиц под воздействием градиента температуры, оценивались из эксперементальных даных. 
Предлагаемые результаты являются очень важными, учитывая применение в современной медицин- 
ской практике (криотерапии, онкологии и других областях) методов гипер- и гипотермии. 
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EXPERIMENTAL MODELLING ОЕ NONLINEARITY 
EFFECTS ALTERNATIVE MAGNETIC FIELDS AND 
LIVING ORGANISMS INTERACTION 

М.У. NETSVETOV, Р.К. KHIZHENKOV, У.І. AVSYANKIN 


Abstract. We propose the experimental model of interaction of polarised molecules 
and varying magnetic fields based on the symmetry of the Maxwell laws and the scales 
invariantly principles. In this paper showed the results of research of 
magnetopolarisating particles condensation from ferrofluidisating layer at the ampli- 
tude reduction on frequencies 1; 2,5; 5; 10; 30 Hz. It is showed a non-linear character of 
condenses forming with various particle numbers. Authors investigated a character of 
particles precipitating and overcoming of potential barrier rate. Obtained dependen- 
cies have at first case 2 maximums and at second — 1. The paper has a biological inter- 
pretation of the results. 

Ключевые слова: поляризованные биомолекулы, переменное магнитное поле, 
экспериментальные модели 


Ранее в ряде работ был показан нелинейный, волнообразный характер 
частотной зависимости влияния переменных магнитных полей H, Ha био- 
логические объекты [1, 4-7, 12, 14]. Известно также, что в низко- и сверхниз- 
кочастотном диапазоне Н, генерирует в проводящих материалах ЭДС и To- 
ки проводимости [11, 13]. Однако объяснить разнонаправленность дейст- 
вия H, на течение биологических процессов (например, регенерации) 
только с позиций индуцированных полем ЭДС и токов проводимости в жи- 
вых клетках и тканях не представляется возможным. Согласно данным 
[11, 13] в изотонических растворах хлорида натрия ЭДС и токи проводимо- 
сти прямопропорциональны частоте и амплитуде приложенного поля и 
имеют тот же порядок величин, что и у регистрируемых в живых тканях — 
от единиц до десятков мкА. 

Целью настоящей работы является создание и исследование экспери- 
ментальной модели, характеризующей взаимодействие поляризованных 
молекул (цвиттерионов) с переменными токами различной частоты и ам- 
плитуды. С точки зрения физики основу модели составляет симметрия за- 
конов Максвелла и принципы масштабной инвариантности. Поскольку 
взаимодействие магнитополяризованных частиц с магнитными полями и 
электрополяризованных частиц с электрическими внешне ничем не отли- 
чаются, использование подобного принципа моделирования можно счи- 
тать вполне оправданным. Такие свойства ферровзвесей как конденсация 
частиц в постоянных, переменных или вращающихся магнитных полях при 
надкритических значениях частот характерны всем масштабным уровням 
— от милли- до нанометровых размеров частиц [3]. Намагниченные и 
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меченые по полюсам сферы диаметром от 3-5. до 10 мм многократно ис- 
пользовали для моделирования различных процессов происходящих в на- 
норазмерных системах [2, 8-10]. 

В той же мере масштабная инвариантность свойственна и электрически 
поляризованным частицам. Таким образом, по результатам взаимодейст- 
вия намагниченных частиц с переменными магнитными полями мы можем 
судить о характере взаимодействия поляризованных частиц или молекул 
(например, аминокислот) с переменными токами. 

Рассмотрим магнитоожиженный (МО) слой, на который действуют од- 
новременно два переменных магнитных поля H', иН’, с частотами f, nf, Ec- 
ли Af достаточно велика (f, > 2f,), эффективность магнитоожижения будет 


такой же, как при действии одного поля с частотой f,. При малых значениях 
ДЕМО слой будет испытывать пульсации: синфазность вершин волн образу- 
ет пучности (суммирование амплитуд H,’ и Н,"), а противофазность — y3- 
ль c Н, =0. Частоту появления пучностей и узлов можно регулировать COOT- 
ношением величин f, и f, например, при f,=8 Гц и f,=7 Гц, она будет равна 
1 Гц, при 49 и 42 Гц соответственно 7 Гц. В узлах интенсивность магнито- 
ожижения резко падает вплоть до полной конденсации, а на пучностях 
слой вновь “вскипает”. Если f,=f,, но существует сдвиг во времени на поло- 
вину периода колебания Н, суммарное поле равно нулю, и эффект магнито- 
ожижения не наблюдается. 

В условиях f, «f, B зависимости от абсолютных величин частот обра- 
зующихся в узлах конденсат частиц будет отличаться по структуре. 
Поскольку чередование “кипящей” и конденсированной фазы является 
быстротекущим процессом, для детального рассмотрения структур ис- 
пользовали метод отключения поля. Шарики в конденсате делились на обо- 
собленные («мономеры») и агрегированные («полимеры»), а соотношение 
агрегатов с разным числом частиц менялось в зависимости от исходной час- 
тоты поля. Количество димеров и тримеров с ростом частоты уменьшалось. 
Число тетрамеров имело максимум при f=10 Гц. Пентамеры при 2,5 Гц 


0,1 1 10 2030 


Рис. 1. Интенсивность “выхватывания” частиц постоянным магнитным полем из магнитоожижен- 
ного слоя в зависимости от частоты переменного поля при амплитудахН, = 80(°)и 120 З (e). 
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проходили через минимум, и далее, с ростом частоты, их количество росло. 
Гексамерам также характерен минимум при 2,5 Гц, a при #=20 Гц ихколиче- 
ство проходило через максимум и далее вновь уменьшалось. Гептамеры 
при f<5 Гц не образовывались и число их максимально при f=10 Гц. 
Восьми- и девятизвенные агрегаты появлялись, начиная с f=20 Гц, а 10-ти 
=12-тис 30 Ги. 

В следующем эксперименте исследовали интенсивность осаждения 
частиц на активный субстрат — постоянный магнит (H,=2 K9), располо- 
женный на крышке контейнера. Чтобы оказаться в зоне притяжения поля 
магнита, частицам необходимо подняться от основания контейнера на вы- 
соту не менее 15-20 мм. На рис. 1 показана интенсивность «выхватывания» 
частиц полем H, из магнитоожиженного слоя. Общее число частиц — 100; 
М-средниє, по пяти измерениям для каждой частоты, количества частиц, 
осевших на субстрат за 1 минуту. При амплитуде Н,=80 Э при увеличении f 
наблюдалось сначала медленное увеличение М, затем быстрый рост, прохо- 
ждение через максимум в интервале Гот 0,5 до 2,0 Гц и минимум при 
f=4-5 Гц, вновь небольшое увеличение (при 6-25 Гц) и далее резкий спад до 
нуля при f=35 Гц. Увеличение амплитуды Н, до 120 Э в целом мало влияет 
на качественную картину процесса: основной максимум эффективности 
приходится Ha те же частоты, однако последующий за ним минимум смеща- 
ется с 4—5 на 2 Гц, а второй максимум возрастает по интенсивности. Кроме 
того, появляется еще один максимум на частоте 0,3 Гц. 

Для исследования частотной зависимости скорости У преодоления час- 
тицами потенциального барьера контейнер снабжали перегородкой высо- 
той 10 мм, а все частицы (40 шт.) изначально помещали по одну сторону. 
У — среднее по результатам пяти опытов число частиц, преодолевших 
барьер за 1 минуту на каждой частоте. Зависимость V(f) приведена на 
рис. 2. Рост { сначала сопровождался резким увеличением V, которая про- 
ходила через максимум при {=8-20 Гц, и последующим постепенным CHH- 
жением до нуля. 


М.мин" 
20 
| 
15 
10 
5 
f, Гц 
1 10 50 


Рис. 2. Зависимость скорости преодоления частицами потенциального барьера от частоты. 
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Из приведенных экспериментальных результатов видно, что все эффек- 
ты взаимодействия частиц между собой и c H, в зависимости от частоты HO- 
сят сложный нелинейный характер. Все описанные эффекты — преоблада- 
ниетого или иного типа агрегатов при конденсации, интенсивность взаимо- 
действия с активным субстратом (рис. 1), преодоление потенциального 
барьера (рис. 2) — имеют один, два и более максимумов (или минимумов), 
приходящихся на разные частотные интервалы. 

Перенеся полученные на моделях результаты на биосистему, можно 
дать некоторую качественную интерпретацию механизмов, лежащих в OC- 
нове разнонаправленности действия поля на биообъект при различных 
частотах. 

1. Если рассмотреть магнитоожижение при действии одной частоты f, 
как модель теплового движения молекул в растворе, то дополнительное 
воздействие второй частотой Ё, будет периодически усиливать и ослаблять 
интенсивность этого движения. Это в первом случае, т.е. в пучностях, долж- 
но способствовать более интенсивному перемешиванию частиц (молекул) 
в растворе и более частому попаданию их B нужное место и в нужное время, 
а во втором случае, в узлах — облегчать включение биомолекул (например, 
аминокислот) в процесс синтеза белков. Совпадение частот по величине 
(f,=f,), но при наличии сдвига во времени 0,5 периода колебания движение 
частиц прекратится, т.е. создадутся условия, препятствующие синтезу 
биополимеров. 

2. С позиций влияния индуцированного полем переменного тока на 
скорость синтеза белка можно рассматривать и экспериментальную зави- 
симость интенсивности осаждения частиц на активный субстрат от часто- 
ты. Здесь наблюдается два и даже три максимума, т.е. переменный ток мо- 
жет и способствовать, и препятствовать синтезу белка на матричной РНК. 

3. Процессы синтеза биополимеров невозможны без осуществления ак- 
тивного транспорта веществ, который, как правило, связан с преодолением 
одного или нескольких потенциальных барьеров, в первую очередь мем- 
бран. Здесь, согласно полученной модели, в живом организме также дол- 
жен наблюдаться максимум скорости трансмембранного переноса молекул 
приходящийся на определенный интервал частот наведенного переменно- 
го тока. 

Приведенные в настоящей работе экспериментальные модели возмож- 
ного механизма действия переменного электромагнитного поля на биоло- 
гические объекты и происходящие в них процессы дают лишь качественное 
объяснение разнонаправленности эффективности при различных часто- 
тах. Подтвердить или опровергнуть правильность предложенной модели 
можно только проведением дополнительных исследований. В первую оче- 
редь это физическое моделирование взаимодействия переменных электри- 
ческого и магнитного полей с цвиттерионами и получение количественных 
оценок. И второе — это биологическое моделирование. Если эксперимен- 
тальная модель отвечает действительности, то, во-первых, одновременное 
действие на живой организм двух магнитных полей с различными частота- 
ми должно быть более эффективным, чем одного, и, во-вторых, интервалы 
биологически активных частот должны сдвигаться в сторону больших или 
меньших значений в зависимости от температуры. 
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ЕКСПЕРИМЕНТАЛЬНЕ МОДЕЛЮВАННЯ НЕЛІНІЙНИХ ЕФЕКТІВ ВЗАЄМОДІЇ ЗМІННИХ 
МАГНІТНИХ ПОЛІВ З ЖИВИМИ ОРГАНІЗМАМИ 


М.В. НЕЦВЕТОВ, П.К. ХИЖЕНКОВ, В.І. АВСЯНКІН 


Запропонована експериментальна модель взаємодії поляризованих молекул із змінними магніт- 
ними полями. За основу моделі прийнята симетрія законів Максвела. Подані результати дослідів кон- 
денсації магнітополяризованих часток з магніторозрідненого шару при зменшенні амплітуди коливань 
на частотах 1; 1,5; 5; 10; 20 та 30 Гц. Показано нелінійний характер створення конденсатів з різним чис- 
лом часток. Досліджено процес осадження часток на активний субстрат та швидкість подолання ними 
потенційного бар'єру. Отримані залежності мають у першому випадку 2 максимуми, а у другому — 1. 
Наведено біологічну інтерпретацію отриманих результатів. 


ЭКСПЕРИМЕНТАЛЬНОЕ МОДЕЛИРОВАНИЕ НЕЛИНЕЙНЫХ ЭФФЕКТОВ 
ВЗАИМОДЕЙСТВИЯ ПЕРЕМЕННЫХ МАГНИТНЫХ ПОЛЕЙ С ЖИВЫМИ ОРГАНИЗМАМИ 


М.В. НЕЦВЕТОВ, П.К. ХИЖЕНКОВ, В.И. АВСЯНКИН 


Предложена экспериментальная модель взаимодействия поляразованных молекул с перемен- 
ными магнитными полями. В основу модели положена симметрия законов Максвелла и принципы мас- 
штабной инвариантности. Приведены результаты исследования конденсации магнитополяризованных 
частиц из магнитоожиженного слоя при уменьшении амплитуды колебаний на частотах 1 ;2,5;5;10;20 и 
30 Гц. Показан нелинейный характер образования конденсатов с различным количеством частиц. Ис- 
следован процесс осаждения частиц на активный субстрат и скорость преодоления ими потенциально- 
го барьера. Полученные зависимости имеют в первом случае 2 максимума, во втором — 1. Приведена 
биологическая интерпретация полученных результатов. 
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Abstract. In this paper we used the method Sit'ko-MRT for preventive maintenance 
and treatment of paresis digestive tract after operations on the colon. It is shown, that 
after the course treatment the faster and proof normalization contractile function in- 
testine was observed in comparison with traditional therapy. The medical effect was 
accompanied by restoring of clinical and biochemical parameters of blood, that ap- 
peared in cupping of endotoxicosis symptoms. The method Sit'ko-MRT was easily 
tolerated by the patients, did not cause side-effects and complications during observa- 
tion. 

Ключевые слова: послеоперационный парез кишечника, микроволновая 
резонансная терапия, эндотоксикоз 


Вступление 


Послеоперационные парезы и параличи дигестивного тракта остаются 
актуальной проблемой абдоминальной хирургии, поскольку являются вто- 
рым по частоте осложнением послеоперационного периода и развиваются 
в 30-40% случаев [1]. 

Из множества факторов, имеющих тормозящее влияние на моторику 
ЖКТ в условиях анестезии и хирургической травмы, наиболее значимым 
является гиперактивность симпатического звена вегетативной нервной 
системы, активизация которой во время операции угнетает желудочно- 
кишечную пропульсию и способствует сокращению кишечных сфинктеров 
[7]. Исследования послеоперационного состояния моторики желудка, 
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MPT в профилактике и лечении парезов желудочно-кищечного тракта после операций на толстой кишке 


двенадцатиперстной и толстой кишки показали, что наиболее стойко тор- 
мозилась электрофизиологическая и моторная активность левой половины 
толстой кишки (до 72 часов), что задерживает опорожнение кишечника до 
3-4 суток и более [8]. Следовательно, послеоперационная атония пищева- 
рительного тракта в первую очередь является проблемой толстой кишки. 

Экстрамуральная иннервация толстой кишки осуществляется симпати- 
ческими и парасимпатическими отделами нервной системы, причем пара" 
симпатические ветви активируют моторику кишки, а симпатические — yr- 
нетают. Исследованиями установлено, что местные рефлексы, возникаю- 
щие при воздействии на точки акупунктуры (ТА), являются активными 
стимуляторами моторной деятельности кишечника, если они возбуждают 
парасимпатическое звено [4]. В последнее десятилетие интерес клиници" 
стов различных специальностей привлекает метод микроволновой резо- 
нансной терапии (МРТ), основанный на облучении ТА низкоинтенсивным 
электромагнитным излучением миллиметрового диапазона. В литературе 
приведены положительные результаты влияния МРТ на процессы санорге- 
неза в различных областях медицины [6]. 

Для обоснования рекомендаций использования МРТ в лечении боль- 
ных хирургического профиля целесообразно привести сведения по некото- 
рым сторонам механизма его действия. Прежде всего, принципиальной 
особенностью современных технологий квантовой медицины является 
крайне низкий уровень внешних воздействий на организм пациента [5]. В 
частности, в микроволновой резонансной терапии устойчивый терапевти- 
ческий эффект наблюдался при интенсивностях внешнего электромагнит- 
ного облучения соответствующего участка тела пациента, сравнимых с та- 
ковыми для тепловых шумов окружающей среды, спектральная плотность 
которых в терапевтическом (миллиметровом) диапазоне длин волн при 
комнатных температурах составляет величину порядка 107" Вт/Ги cM. 

Для адекватности используемого воздействия необходимо оценить CO- 
стояния меридианов, соответствующих облучаемым ТА. Для этого предло- 
жен метод пульсовой диагностики, в котором определяются характеристи- 
ки пульсовой волны на участке цунь-коу с использованием компьютерного 
пульсодиагностического комплекса [3]. 

Целью работы явилось повышение эффективности профилактики и ле- 
чения парезов кишечника у больных после операций на толстой кишке пу- 
тем использования в до- и послеоперационном периоде метода МРТ в соче- 
тании с компьютерной пульсовой диагностикой. 


Объект и методы исследования 


В исследование включены результаты лечения 114 больных, опериро- 
ванных в отделе абдоминальной онкологии УНИИОР по поводу различных 
патологических процессов ободочной и прямой кишки, из них 77 пациента 
составили исследуемую группу. Им с целью предупреждения и лечения по- 
слеоперационной атонии кишечника использовали МРТ. Контрольная 
группа включала 37 больных, которым B XO- и послеоперационном периоде 
не проводили подобные мероприятия. В обеих группах больные были ран- 
домизированы по полу, возрасту, нозологическим формам методом 
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лотерейного отбора. В табл. 1 представлено распределение больных по ви- 
дам оперативных вмешательств. 


Таблица 1. Распределение больных по видам оперативного вмешательства 
Г T Т ] 


В перативного вмешательства | Исследуемая | Контрольная 
| nA onep | группа группа 
| Правосторонняя гемиколэктомия | 15 | 7 
ji | 
| | | 
| Левосторонняя гемиколэктомия | 5 2 
| 
Резекция сигмовидной кишки 9 4 
Передняя резекция прямой кишки | 11 5 
Брюшно-анальная резекция прямой | - 13 6 
| кишки | 
f Т 2 
| Операция Кеню-Майлса | 12 | 6 
| ро - i 
| | | 
| Реконструктивная колопластика | 12 | 7 
| | 
| | 
Всего | 77 | 37 
| i = 


С целью стимуляции моторики ЖКТ назначали МРТ, начиная с 2-су- 
точного послеоперационного периода с профилактической целью (10 еже- 
дневных сеансов). 

МРТ проводили 26 больным исследуемой группы в предоперационном 
периоде: 12 больным с функционирующими колостомами, 14 больным с 
субкомпенсированной формой хронического колостаза на фоне обтюри- 
рующей опухоли толстой кишки. 

Источником миллиметрового излучения служили генераторы электро- 
магнитного излучения АМРТ-01, АМРТ-02. Продолжительность воздейст- 
вия на ТА составляла 5 минут. За один сеанс облучали от двух до четырех TO- 
чек акупунктуры, в зависимости от энергетического состояния меридиа- 
нов. Энергетику меридианов оценивали с использованием компьютерного 
пульсодиагностического комплекса. Курс лечения — 8-10 процедур. 

Степень эндотоксикоза определяли по величине лейкоцитарного ин- 
декса интоксикации (ЛИИ) [2], показателям общего белка, мочевины, креа- 
тинина, электролитов в сыворотке крови. Полученные данные обрабатыва- 
ли методом вариационной статистики с использованием і-критерия 
Стьюдента. 


Результаты и их обсуждение 


В послеоперационном периоде у больных исследуемой группы приме- 
нение МРТ уже после двух процедур способстьовало купированию явле- 
ний гастростаза, снижению симптомов эндоген 1 интоксикации (субфеб- 
рилитет, тошнота, слабость), а полное вг ·,~ -овление перистальтики 
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наблюдали через 3,6=1,7 суток после операции. У больных контрольной 
группы, получавших стандартное лечения (инфузионная терапия, прозе- 
рин), перистальтика кишечника восстановилась через 5,5+2,1 суток после 
оперативного вмешательства (Р» 005). Динамика ЛИИ в исследуемой и 
контрольной группах свидетельствует о более быстром восстановлении 
этого показателя у больных, прошедших курс МРТ (табл. 2). 


Таблица 2. Динамика показателей ЛИИ у больных в послеоперационном пе- 
риоде 


Время после операции | Исследуемая группа | Контрольная группа 
| | 
| 5.451,20 | 5195112 | 
| 4-е сутки | 2,50+0,36 | 4,750,61* | 


1-е сутки 


* _Р<0,05 по сравнению с исследуемой группой на 4-е сутки. 


Биохимические показатели крови в динамике (креатинин, мочевина и 
электролиты плазмы) у больных исследуемой группы также имели более 
выраженную тенденцию к нормализации по сравнению с контрольной 
группой. 

Негативных побочных реакций и осложнений при использовании МРТ 
у больных исследуемой группы не отмечено, переносимость метода была 
оценена как удовлетворительная всеми пациентами. 

Исходя из полученных результатов можно заключить, что включение 
МРТ в послеоперационное комплексное лечение больных, перенесших 
вмешательства на толстой кишке, позволяет эффективно купировать ато- 
нический синдром со стороны ЖКТ, положительно влияет на востановле- 
ние клинических и биохимических показателей крови, проявляя тенден- 
цию к более быстрому снижению эндотоксикоза. 


Авторы выражают благодарность проф. С.П. Ситько за поддержку и интереск работе. 


МІКРОХВИЛЬОВА РЕЗОНАНСНА ТЕРАПІЯ У ПРОФІЛАКТИЦІ І ЛІКУВАННІ ПАРЕЗІВ 
ШЛУНКОВО-КИШКОВОГО ТРАКТУ ПІСЛЯ ОПЕРАЦІЙ НА ТОВСТІЙ КИШЦІ 


О.П. КУЗЬМЕНКО, І.Є. СОЛОВИОВ, А.В. ТОФАН 


У запропонованій клінічній роботі використаний метод Сітько-МРТ при профілактиці і ліку- 
ванні парезів дігестивного тракту після операцій на товстій кишці. Показано, що під впливом курсового 
лікування спостерігалася більш швидка і стійка нормалізація скорочувальної і пропульсивної функції 
кишечника порівняно з традиційною терапією. Лікувальний ефект супроводжувався відновленням 
клінічних і біохімічних показників крові, що виявлялося в подоланні симптомів ендотоксикозу. Метод 
Сітько-МРТ легко переносився пацієнтами, не викликав побічних ефектів і ускладнень y період спосте- 
реження. 
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МИКРОВОЛНОВАЯ РЕЗОНАНСНАЯ ТЕРАПИЯ В ПРОФИЛАКТИКЕ И ЛЕЧЕНИИ 
ПАРЕЗОВ ЖЕЛУДОЧНО-КИШЕЧНОГО ТРАКТА ПОСЛЕ ОПЕРАЦИЙ НА ТОЛСТОЙ 
КИШКЕ 


А.П. КУЗЬМЕНКО, И.Е. СОЛОВЬЕВ, A.B. ТОФАН 


В представленной клинической работе использован метод Ситько-МРТ при профилактикеиле- 


чении парезов дигестивного тракта после операций Ha толстой кишке. Показано, что под влиянием кур- 
сового лечения наблюдалось более быстрая и стойкая нормализация сократительной и пропульсивной 
функции кишечника по сравнению с традиционной терапией. Лечебный эффект сопровождался вос- 
становлением клинических и биохимических показателей крови, что проявлялось в купировании сим- 
птомов эндотоксикоза. Метод Ситько-МРТ легко переносился пациентами, не вызывал побочных 
эффектов и осложнений в период наблюдения. 
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НЕКОТОРЫЕ АСПЕКТЫ НАРУШЕНИЯ 
ГОМЕОСТАЗА ПРИ КРИТИЧЕСКИХ СОСТОЯНИЯХ 
И ИХ КОРРЕКЦИЯ С ПОМОЩЬЮ КВЧ-ТЕРАПИИ И 
КСЕНОИМПЛАНТАЦИИ ЭМБРИОНАЛЬНЫХ 
МОЗГОСПЕЦИФИЧЕСКИХ ТКАНЕЙ 


А.В. МАЛЫХИН 


Украинский НИИ клинической и экспериментальной неврологии и 
психиатрии 


SOME ASPECTS OF HOMEOSTASIS VIOLATION 
UNDER HEAVY CONDITIONS AND ITS CORRECTION 
BY THE UHF-THERAPY AND XENOIMPLANTATION 
OF EMBRIONIC BRAIN-SPECIFIC TISSUE 

A.V. MALIHIN 


Abstract. The problems of correlation of a carbohydrate, protein and fatty exchange 
and heavy condition of the patients are investigated. We have shown that the heavy 
condition is stipulated by violations of oxygen utilization in a system. EHF-therapy 
and cerebral structure of a tissue change an ion composition of plasma that is associ- 
ated with to an index of oxygen utilization. 

Ключевые слова: метаболизм, тяжесть состояния, ионный состав, транспорт 
кислорода, фагоцитоз 


Нарушения метаболизма, достигающие угрожающих параметров жиз- 
недеятельности больных являются основной проблемой теоретической и 
практической биологии и медицины [1,2,3]. Методологически проблема 
разделена на изучение этиологических факторов воздействия (первичных) 
и регуляторных особенностей организма с характерной для них нелиней- 
ной фрактальностью функций [3,4,5]. При этом, выявляемые различными 
методами исследований, изменения параметров жизнедеятельности орга- 
низма и их коррекции до физиологической нормы не решают проблемы, 
так как отсутствует единый патогенетический подход ко всему спектру ис- 
следуемых параметров. Перспективным в решении этой проблемы может 
быть подход с позиций диссипативных структур за счет так называемого от- 
рицательного прироста энтропии извне — КВЧ [3,4-7], компенсирующего 
обычную положительную энтропию. В такой ситуации особое значение 
приобретают методы, позволяющие не только выявить закономерности из- 
менения показателей, но и предвидеть возможные перестройки фазовых 
пространств гипотетических моделей рассматриваемых систем [3,5,8-11]. 

Целью работы явилось изучение нейробиологических адаптационных 
механизмов КВЧ-воздействий и ксеноимплантации, включающих исследо- 
вание взаимоотношений между структурой клетки и изучением адаптив- 
ных колебаний ее функционирования, которые взаимосвязаны с динами- 
кой диссоциации комплекса эритроцит-гемоглобин и газовым составом 
атмосферы. | 

Задачи работы: изучение структурно-функциональных изменений ак- 
тивности ферментативных систем; водно-электролитного обмена при 
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возникновении и регрессе синдромов эндогенной интоксикации, обуслов- 
ленной травмой ЦНС и алкоголизацией; клинических проявлений выра- 
женности — гемодинамических (регуляторно-интегративных) и метаболи- 
ческих (трофически защитных) расстройств с формированием соматонев- 
рологических синдромов. 


Материал и методы исследований 


Исследовано 150 больных в возрасте от 20 до 60 лет с различной степе- 
нью выраженности синдрома эндогенной интоксикации. Контрольная 
группа — 50 практически здоровых лиц (I группа). Проведен комплекс ин- 
струментальных (ЭЭГ, ЭКГ, РЗГ, КТ-головного мозга), иммунологических 
и биохимических исследований. 

Обследуемые больные были разделены на две группы в зависимости от 
этиологического фактора и выраженности клиники соматоневрологиче- 
ских проявлений синдрома эндогенной интоксикации с учетом изменений 
периферического состава крови, ионного состава, активности кининаз. 
Морфометрические показатели головного мозга оценивали по данным 
компьютерной томографии. 


Результаты исследований 


Показатели компьютерных исследований у обследуемых I группы ис- 
пользованы в качестве нормы для размеров ликворосодержащих путей: ин- 
декс передних рогов (25,41 +1,23)%; индекс центральных отделов 
(24,31=0,52)%; индекс задних рогов (37,41 = 0,19)%; ширина третьего желу- 
дочка (4,13+0,13)см. 

П группа: индекс передних рогов (28,23+0,71)%; индекс центральных 
отделов (32,52+1,53)%; индекс задних рогов (40,21+1,52)%; ширина 
третьего желудочка (5,62 —0,31)см. 

Ш группа: индекс передних рогов (33,2+1,4)%; индекс задних рогов 
(44,14+1,83)%; ширина третьего желудочка (7,12+1,25) см. 

Общий кининоген І группы (4,1::0,2) мкг/мл, кининазная активность 
(1,0+0,2) мкг/мл, зстеразная активность каликреинов — (7,8+3,0) мкг/мл. 
Содержание холестерина — 5+0.3 ммоль/л, содержание фосфолипидов — 
4.1=0.2 ммоль/л. | 

П группа: содержание холестерина — 6+0.4 ммоль/л, содержание boc- 
фолипидов — 4+0.3 ммоль/л., кининазная активность (-1.5=0.1), эстераз- 
ная активность каликреинов — 9+3 мкг/мл. 

Ш группа: содержание холестерина — 3.5+1.3ммоль/л, содержание 
фосфолипидов — 3+0.3ммоль/л., кининазная активность — 1.7+ 0.1, зсте- 
разная активность каликреинов — 12-3 мкг/мл. 

Ударный объем крови в среднем по I группе (95,1+3,1) мл; по 
П (89,2=5,2) мл; по Ш (87,1+6,1) мл. Ударный индекс в среднем по І группе 
(51,4=2,1) war; по II (50,8=2,5) мл/л; по Ш (47,2+3,1) мл/л?. 

Содержание натрия в крови в среднем по I группе (141,26=0,51) 
ммоль/л; по II (138,44+2,51) ммоль/л; по Ш (149,82+2,14) ммоль/л. Содер- 
жание натрия в моче по І группе (65,94+ 15,3) ммоль/л; по II (110,13+8,55) 
ммоль/л; по Ш (120,52+15,46) ммоль/л. 

Концентрация калия в моче больных I группы в среднем составляла 
(35,63+4,99) ммоль/л; П (35,94+5,09) ммоль/л; у Ш (27,63=3,97) ммоль/л. 
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Некоторые аспекты нарушения гомеостаза при критических состояниях и их коррекция 


Отношение содержания калия мочи к содержанию натрия плазмы состав- 
ляло соответственно 0.26, 0.25, 0.21. Удельный вес мочи в среднем по группе 
составлял: для 1 (100642) ma/cm’; II (1009+2) мл/см'; Ш (1011=5) мл/см.. 

Содержание креатинина плазмы по группам было следующим: в 
I (73,3=3,1) мкмоль/л, II (92,3+5,6) мкмоль/л, Ш (51,10+3,76) мкмоль/л; в 
моче соответственно по группам: (0,09+0.02) мкмоль/л, (0,14=0.02) 
мкмоль/л, (0.16+0.02) мкмоль/л. 

Изучение клинической картины больных всех групп показало, что кли- 
ника синдрома эндогенной интоксикации взаимосвязана с лейкоцитарной 
формулой крови, которая определяет образование перекисей жирных ки- 
слот. Формула крови регулирует процессы окисления-восстановления, 
протекающие в печени, почках, скелетной и сердечной мышцах. Процесс 
окисления взаимосвязан с рядом ферментативных реакций, протекающих 
в тромбоцитах и лейкоцитах. Развитие синдромов эндогенной интоксика- 
ции с нарушениями центральных и периферических звеньев регуляции 
кровообращения (лимбико-ретикулярные, ренин-ангиотензиновые, кар- 
диальные и сосудистые механизмы) связаны с активацией полипептидов, 
действующих на уровне периферических звеньев кровообращения, необ- 
ходимых в обеспечении доставки кислорода к клеткам паренхимы, а также 
в удалении из тканей метаболитов. Скорость тканевых процессов соответ- 
ствует электрохимическим процессам на уровне микроциркуляции и взаи- 
мосвязана с энтропией энергии движущейся крови для приведения в дейст- 
вие фильтрационного механизма. Временные различия обусловлены изме- 
нениями состава периферической крови к пролиферативным процессам, 
происходящим в костном мозге — 490” [13,14]. Клеточные мембраны име- 
ют частоты собственных колебаний в диапазоне 10"-10" Гц. 

Теоретически рассчитанные значения энтропии организма на основе 
стохастического математического анализа логических функций показыва- 
ют, что данные энергетических изменений в изучаемой биосреде можно 
описать с помощью теоремы Умова-Пойтинга: 


Jpeenps - J(5-Eyar - fs. E)av +H Zay -0, 


предположив, что энергетические процессы связаны C морфометрической xa- 
рактеристикой капиллярного русла, числом плотности капилляров, объемом 
капилляра, концентрацией глюкозы, холестерина, фосфолипидов и их изме- 
нениями вединицу времени. При этом можно предположить, что энергия по- 
глощается внутри объема биосреды, которая включает капилляр, интерсти- 
циальное пространство, с наличием ограниченной замкнутой поверхности $ 
(эритроциты). При таком условии скорость кровотока в системе микроцирку- 
ляции должна соответствовать: 


Su = p, V, ‘Su: YO а 0,063 мл / сек, 
(CA, -D,) «Ро "os :60,,. 


где P, — атм. давл.; У, — объем капиллярной сети; $, — сечение капиллярной 
сети; СА, — артериальное давление; D, — диастолическое давление; У, — 
объем зритроцита; 5,, — площадь эритроцита и 60,,, — временной интервал, 
УО — ударный объем. В таком случае первый интеграл выражения представ- 
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ляет собой мощность, переносимую биологической средой (жидкостью) в ре- 
зультате сокращения миокарда и продвижения ударного объема черезедини- 
цу поверхности биосреды. Второй интеграл — это скорость электрических 
процессов, обусловленная концентрацией ионов в микроциркуляторной сис- 
теме (сечение капиллярной сети 2,98 -10° - 3,13 -10°; среднее значение дли- 
ны капилляра 0,057 см, диаметр 0,0086 см и объем капилляра 11.10). Cym- 
марным выражением всех процессов является единица электромагнитного 
поля, аналогичная единице вещества — молю. Один моль содержит 6.02 .102 
частиц и составляет 399 кДж/моль или 0,399 Вт/сек [12]. Он определяет напря- 
жение газов в крови (О, и CO,) и ход температурных реакций. Изменения 
монитируемых показателей содержания ионов в единице объема биосреды 
связаны со скоростью диссоциации и восстановления комплекса эритро- 
цит-гемоглобин, которые обусловлены ходом биохимических свободно-ра- 
дикальных реакций, протекающих в 8 этапов. Ход реакций регулируется ско- 
ростью электрохимических процессов, концентрацией ионов, уровнем би- 
карбоната, концентрацией холестерина, фосфолипидов и белков, которые 
регулируют проницаемость мембран эритроцита и сосудистой стенки микро- 
циркуляторного звена кровообращения в зависимости от количественного и 
качественного состава лейкограммы. В общем виде скорость диссоциации 
комплекса эритроцит-гемоглобин (СДКЭ) может быть представлена следую- 
щим образом [15]: 


С, бы -Р, 1006 _ OCD(Na + NaHCO, ) 


т YO.(P+S)% | My, 


=50 +3 сек, 


где Му, — молекулярный вес натрия; P+S — содержание палочко- и 
сегментоядерных нейтрофилов, OCD — осмотическое давление, С, — клеточ- 
ная пролиферация. 

Коэффициент утилизации кислорода составил по группам: 
I (31,141,5)%; II (28,0+3)%; Ш (23,2+3)%. При этом во П группе увеличе- 
ние морфометрических размеров мозга отмечалось в 37% случаев, a B III — 
в 43%. Tun развития гемодинамических расстройств взаимосвязан CO CKO- 
ростями синтеза и распада комплекса, концентрации ионов натрия и калия, 
определяющих ход температурной реакции [15]. | 

С. целью коррекции выявленных нарушений системы адаптации на 
больных воздействовали КВЧ с частотным характером 51-65 Ггц в первые 
5 суток и проводили ксенотрансплантацию активных эмбриональных MO3- 
госпецифических тканей, состав которых выбирался в зависимости от ини- 
циирующего и доминирующего звена в системе развития механизма деза- 
даптации. Использовались эмбриональная ткань печени, поджелудочной 
‚ железы в сочетании или отдельно с нейрогенными мозгоспецифическими 
~ (locus coeruleus) тканями, продуцирующими эндогенные биологические 
регуляторы. Доза материала подбиралась индивидуально в зависимости от 
скорости образования и распада белка (концентрации креатинина, плаз- 
менного белка и холестерина, объема эритроцита, удельного веса мочи и 
концентрации липопротеидов). Согласно полученным данным оптималь- 
ная терапевтическая доза трансплантата составила 0,4-0.5 мг на килограмм 
массы. Частотный характер КВЧ-воздействий на первом этапе заболева- 
ния с последующей ксеноимплантацией были оптимальными при контроле 
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активности глютатионовой системы, активности Т-лимфоцитарной систе- 
мы и потребления кислорода. 

Общей биологической закономерностью активации адаптационных 
механизмов для всех групп больных при благоприятном течении синдрома 
эндогенной интоксикации было увеличение количества лейкоцитов при 
нормализации нейтрофилов за счет снижения палочкоядерных форми уве- 
личениє количества лимфоцитов (процентного содержания). Диапазон 
представительства тромбоцитов сужался от 120000 до 250000 ед. Уровень 
кининогена повышался до 4 мкг/мл при понижении кининазной активно- 
сти меньше 1 мкг/мл и эстеразной активности каликреинов до (6=2) мкг/мл 
в минуту. Коэффициент множественной корреляции (КМК) указывает на 
наличие достоверной связи сегментоядерных нейтрофилов со всей сово- 
купностью параметров кислотно-щелочного баланса, гемодинамики 
(КМК=0.757, Е=6.23) посредством связи с числовой величиной комплекса 
эритроцит-гемоглобин. 

Анализ исследований показал, что соотношения клеточного состава 
крови и коэффициента растворимости кислорода реализуются B многопро- 
фильные изменения ферментативных систем на органных уровнях (почек, 
печени, кардиореспираторной системы) гемоконцентрационных показате- 
лей и в изменения содержания электролитов плазмы и мочи, которые обу- 
словлены изменением объема циркулирующей крови, объема циркулирую- 
щей плазмы и минутного объема кровообращения. 

Надо полагать, что изменение этих параметров происходит при участии 
всех составляющих клеток крови с определенной их специализацией, опре- 
деляющих иммунологическую и биохимическую направленность указан- 
ных реакций на системном и органном уровнях в зависимости от наличия 
во внутренних органах метаболических субстратов, регулирующих akTHB- 
ность ферментов но принципу обратной связи. В основе этих обратных свя- 
зей, как свидетельствуют результаты биохимических инструментальных и 
иммунологических исследований находятся механизмы свободноради- 
кальных процессов фагоцитоза. | 


Обсуждение результатов 


Биологические механизмы КВЧ-воздействия на организм и ксе- 
нотрансплантации включают общие вопросы адаптации, основанной на 
принципе саморегуляции. Не вызывает сомнений тот факт, что в адаптив- 
ных реакциях особая роль принадлежит механизмам кровоснабжения, ре- 
гулируемым ЦНС, периферическими звеньями кровообращения, а также 
микроциркуляторным звеном кровообращения внутренних органов, обес- 
печивающих выживаемость организма в экстремальных условиях. 

Как показали наши исследования, адаптационный системный процесс 
разворачивался отсрочено во времени на ионном и молекулярном инфор- 
мационном уровнях, взаимообусловленных первичным информацион- 
но-колебательным алгоритмом, возникающим посредством взаимосвязи 
миотического деления эмбриональных клеток со скоростью образования и 
потребления кислорода, скоростью диссоциации комплекса эритро- 
цит-гемоглобин, регулируемого составом крови и газовым составом ат- 
мосферы, определяющим функционирование специализированных Op- 
ганных ферментативных систем и возникающих при этом условий 
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функционирования периферических гемодинамических систем. Отклоне- 
ния функционирования органов или систем органов определяют измене- 
ния активности лимбико-ретикулярного комплекса, ренин-ангеотензино- 
вой и каликреин-кининовой систем. 


Выводы 


Системные адаптационные механизмы КВЧ воздействий охватывают 
все регуляторные звенья и через определенный временной интервал опре- 
деляют адекватную систему потребления кислорода. 

Ксеноимплантация эмбриональных мозгоспецифических тканей воз- 
действует на адаптационные системы отсрочено (через 6-7 суток). 

Взаимосвязь между количественными изменениями водных сред 
организма и, как следствие, возникновение изменений концентрации белка 
креатинина направлены на адаптацию транспорта и потребления 
кислорода. 

Начальные механизмы адаптации включают изменения перифериче- 
ского состава крови, вызывающие изменения ионного состава плазмы, уро- 
вень регуляции в системе микроциркуляторного механизма кровообраще- 
ния с включением ферментативной активности почек и печени. 

Нарушение транспорта кислорода определяет степень выраженности 
клиники критических состояний, которая зависит от формирующихся из- 
менений кровообращения внутренних органов. Биохимический этап адап- 
тации достигается посредством изменения колебательно-информацион- 
ной характеристики системы, объединяемой нервной и гуморальной регу- 
ляцией, которые участвуют в формировании значений величины 
сердечного выброса, ударного объема, параметров водно-солевого обмена, 
общего сосудистого сопротивления и тем самым влияют на количественное 
перераспределение объема циркулирующей крови по органам: печень, 
почки, сердце, головной мозг, мышцы, остальные органы. Это подтвержда- 
ет существующие представления, что изменения кровоснабжения внутрен- 
них органов служат триггерным механизмом включения системы калик- 
реинов, простогландинов, взаимосвязанных с образованием и потреблени- 
ем кислорода. 

Спектр биохимических реакций предопределяет регуляцию окисли- 
тельно-восстановительных процессов посредством активации иммуноло- 
гической системы, неразрывно связанной с колебаниями гуморальных 
сред. Все эти процессы осуществляются через стромально-паренхиматоз- 
ные взаимодействия, посредником при которых является периферическая 
кровь (носитель активации ферментов специализированных органов). По- 
средством изменения частотного характера мембраны эритроцита, пере- 
дающегося путем сокращения белковых миофибрилл через базальные мем- 
браны микроциркулирующего русла на интерстициальные пространства, 
достигается биохимическая и биофизическая сбалансированность функ- 
ций внутренних органов в их взаимоотношениях с интерстициальным про- 
странством и внешней средой. Конечным результатом вышеуказанного 
процесса является сохранение постоянства внешней и внутренней среды, 
основу которого составляет ионное равновесие. 
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ДЕЯКІ АСПЕКТИ ПОРУШЕННЯ ГОМЕОСТАЗУ ПРИ КРИТИЧНИХ СТАНАХ ТА їхня 
КОРЕКЦІЯ ЗА ДОПОМОГОЮ КВЧ-ТЕРАШІ I КСЕНОІМПЛАНТАЦІЇ ЕМБРІОНАЛЬНИХ 
МОЗКОСПЕЦИФІЧНИХ ТКАНИН 


А.В. МАЛИХІН 


Досліджені питання взаємозв'язку вуглецевого, білкового, жирового обміну з тяжкістю стану 
хворого. Показано, що тяжкість стану обумовлена порушеннями, що виникають у системі утилізації 
кисню. КВЧ-терапія та мозкоспецифічні тканини змінюють іонний склад плазми, що зв'язаний з показ- 
ником утилізації кисню. 


НЕКОТОРЫЕ АСПЕКТЫ НАРУШЕНИЯ ГОМЕОСТАЗА ПРИ КРИТИЧЕСКИХ 
СОСТОЯНИЯХ И ИХ КОРРЕКЦИЯ С ПОМОЩЬЮ КВЧ-ТЕРАПИИ И 
КСЕНОИМПЛАНТАЦИИ ЭМБРИОНАЛЬНЫХ МОЗГОСПЕЦИФИЧЕСКИХ ТКАНЕЙ 


А.В. МАЛЫХИН 


Исследованы вопросы взаимосвязи углеводного, белкового, жирового обмена с тяжестью со- 
стояния больного. Показано, что тяжесть состояния обусловлена нарушениями в системе утилизации 
кислорода. КВЧ-терапия H мозгоспецифические ткани изменяют ионный состав плазмы, который взаи- 
мосвязан с показателем утилизации кислорода. 


ЛИТЕРАТУРА 


1. Е.О. Андреєв, М.У. Білий, С.П. Сітько Проявлення власних характеристичних частот 
організму // Доп. AHYPCP.— СБ. - 1984.— Ne10.— С. 56-59. 

2. П.К. Анохин Очерки по физиологии функциональных систем — M., 1975. 

3. В.П. Казначеев, Л.П. Михайлова Сверхслабые изменения в межклеточных взаимодейст- 
виях. 1981. 

4. Л.В. Недорезов, Б.А. Недорезова [| Автометрия. - 1993.— №2 С. 5-10. 

5. В.И. Петросян и др. Физика взаимодействия мм-волн с объектами различной природы 
// Радиотехника. — №9.— 1999.— С. 20-31. 

6. В.А. Нагорнев Кинетика клеточных элементов сосудистой стенки при атеросклерозе // 
Арх. IIar.— 1928.— №10.— C.— 86-95. 

7. B.3. Ланкин, Л.М. Вихерт Роль перекисного окисления липидов в этиологии и патоге- 
незе атеросклероза // Вопр. Мед. Химии - 1979.— №10.— С. 69-72. 

8 А.В. Малыхин Клинико-компьютерно-реоэнцифалографическая характеристика лик- 
ворно-гипертензионного синдрома у больных с изолированной и сочетательной ЧМТ // 
Український вісник психоневрології. T. 6, випуск 3(18).— 1998.— С. 45-47. 

9. Дж.Е. Коттрел Защита мозга. Анестезия и реанимация. - 1996.— №2.— C.81-85. 

10. Л.Б. Лихтерман, В.Н. Корниенко, A.A. Потапов, В.А. Кузьменко, В.И. Горбунов, 

В.М. Трошин Черепно-мозговая травма: прогноз течения и исходов.— М: Книга ЛТД, 
1993.— 208 с. 

11. В.К. Siesjo Cerebral circulation and metabolism / j. Neurosung., 1984.— Ne10(60).— 
P. 883-908. 

12. В.Е. Илларионов Основы лазерной терапии. M., 1992.— С. 16-17. 

13. Н.Л. Федоров Регуляция пролиферации кроветворных клеток. М., Медицина.— 1977, 
C273: | 

14. Е А. McCullok Stem sells in normal and leukemic hemopoesis.— 1983.— V.62.— N1.— 
P. 1-13. 

15. A.B. Малыхин Метод скрининговой оценки функционального и физиологического со- 
стояния человека в экстремальных ситуациях: Материалы Междунар. Конф. “Доноры 
для охраны здоровья Украины”. — Киев, 1999.— С. 71. 


115 


ИЗУЧЕНИЕ ПЕРОКСИДАЗНОЙ АКТИВНОСТИ 
КРОВИ У БОЛЬНЫХ ЦЕРЕБРАЛЬНЫМ 
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STUDY OF THE BLOOD PEROXIDAZIS ACTIVITY FOR 
PATIENTS WITH CEREBRAL ATHEROSCLEROSIS 
DURING MICROWAVE RESONANCE THERAPY 

V.M. KUZMENKO 


Abstract. We studied the blood peroxidazis activity during the microwave resonance 
therapy (MRT) on 41 patients with an early cerebrovascular pathologic form caused 
by atherosclerosis. Analysing these results we showed that the positive MRT effect re- 
lieved the peroxidation syndrome in patients with the brain blood circulation defi- 
ciency. Thus the blood peroxidazis activity indexes attests that the given method 
provides a useful estimation of treatment efficiency as well as prognosis of MRT suc- 
cess. 

Ключевые слова: церебральный атеросклероз, пероксидазная активность 
крови, микроволновая резонансная терапия 


Известно, что этиология и патогенез многих заболеваний связаны с па- 
тологическими изменениями свободнорадикальных реакций [4,9]. 

Одним из показателей, характеризующих состояние свободно-ради- 
кальных реакций, является пероксидазная ферментативная система, встре- 
чающаяся во всех клетках аэробных организмов [1,5,10]. 

Существует также мнение, что пероксидазная активность является од- 
ним из показателей имунной резистентности организма [3,6,8,11]. 

Рядом исследований показано, что этиология и патогенез цереброваску- 
лярных заболеваний тесно связаны с изменениями иммунной резистентно- 
сти, липидного и жирнокислотного состава крови, что непосредственным 
образом влияет на характер свободно-радикальных реакций [2,3,7,8]. 

Учитывая вышеизложенное, в настоящей работе изучалась пероксидаз- 
ная активность крови и ее изменения в процессе микроволновой резонанс- 
ной терапии (МРТ) у больных церебральным атеросклерозом. 


Материалы и методы исследования 


В данном сообщении представлены результаты определения перокси- 
дазной активности крови у 41 больного цереброваскулярной патологией, 
обусловленной атеросклерозом. Среди обследованных было 11 мужчин и 
30 женщин в возрасте от 30 до 60 лет (средний возраст 40,1+2,4 лет). 

Диагноз цереброваскулярного заболевания устанавливался на ос- 
новании сопоставления данных клинического осмотра с результатами 
параклинических исследований (изучение состояния глазного дна, реоэн- 
цефалография, ЭКГ, термография, допплерография, биохимические 
исследования). 


116 “Physics of the Alive”, Vol. 8, No. 1, 2000 


Изучение пероксидазной активности крови y больных церебральным атеросклерозом в процессе МРТ 


Нами выделены 2 группы 
больных: а) больные с началь- 
ными признаками недостаточ- 
ности мозгового кровообраще- 
ния (НИНКМ) — 30 человек; 
6) группа больных с дисцирку- 
ляторной энцефалопатией I ста- 
дии (11 наблюдений). 

Микроволновая резонанс- 
ная терапия проводилась по 
разработанной нами методике. 
При этом курс лечения состав- 
лял 5-10 сеансов 20 минутного 
воздействия электромагнитно- 
го излучения, мощностью 
1-3 мкВт/см’. 

Пероксидазную активность 
крови определяли до, в процес- 
се и после окончания курса 
МРТ. Исследования проводили 
нахемилюминометре ХЛМЦ-01 
в системе ЭДТА (0,5 mM), фос- 
фатный буфер (рн=6,2) — 
кровь (0,1 мл) по модифициро- 
ванному методу De Toledo [5]. 
Реакция инициировалась введе- 
нием 1,5:107"Мм Н.О.. 

Концентрации всех раство- 
ров приведены до конечного 
объема пробы 4 мл. Гепаринизи- 
рованную кровь разводили B CO- =~ z menep мн 
отношении 1:100 0,9%-ного Puc. 1. Характеристические кривые хемилюминисцен- 
NaCl; для исследования брали ции с низкими (1, 2) и высокими (3, 4) показате- 

" лями пероксидазной активности крови. Шкала 
0,1 мл разведенной крови. Из- самописца 25 · 10? имп/сек. 
меряли светосумму свечения за 
5 минути его интенсивность I, и Ї, По интенсивности свечения І, данной pe- 
акции представляется возможным определить концентрацию пероксидазы 
хрена от 10" M и выше. Для сопоставления с данными, полученными спек- 
трофотометрическим методом, содержание фермента определяли по ка- 
либровочной кривой, а его активность — в пересчете на единицу активно- 
сти в реакции с О-дианизидином. 


Результаты и их обсуждение 


В результате определения пероксидазной активности крови у больных 
цереброваскулярньми заболеваниями установлен ряд особенностей. Так, 
исходные кривые хемилюминесценции имели медленную (L) и быструю (1,) 
вспышки. При этом вспышки (I, и L) отличались различной величиной ам- 
плитуд, их соотношением, а также значением светосуммы (рис. 1). 
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По-видимому, это обусловлено индивидуальными особенностями функ- 
ционирования антиокислительных систем организма, преобладанием фер- 
ментативного или неферментативного биологического окисления, а также 
реактивностью организма в данных условиях. Этим можно объяснить зна- 
чительное варьирование хемилюминесцентных показателей крови. Так, 
интенсивность свечения (L) изменялась независимо от зтиологического 
фактора от 748,1+116,7 имп/сек до 1669,3+134,4 имп/сек, а 5 в пределах от 
96345,7+ 13177 до 242606+ 12552,3 (имп. 5 мин.). В пересчете по реакции с 
О-дианизидином это соответствует, в среднем, 220-660 единиц. Согласно 
данным литературы активность пероксидазы крови практически здоровых 
людей колеблется от 257,5 до 305,7 единиц. 

Кинетика хемилюминесцентной реакции при низкой интенсивности 
свечения имела замедленное развитие. Достижение максимальных значе- 
ний І, происходило за 2-3 мин. (рис. 1, кривые 1, 2), тогда как в группе боль- 
ных с высокими показателями за 0,8-1,2 мин. (рис. 1, кривые 3, 4). 

После 5 сеансов МРТ отмечалась динамика таких хемилюминес- 
центных показателей как кинетика свечения, его уровень и светосумма за 
5 минут. Так, в группе больных с положительным лечебным эффектом 
(улучшение общего состояния, регресс головных болей, головокружения, 
очаговой неврологической микросимптоматики) интенсивность свечения 
увеличивалась в среднем с 1300+ 134,8 имп/сек до 2750+290,7 имп/сек, cBe- 
тосумма с 176248+14178 имп. 5 минут до 212256+15552 имп. 5 минут 
(р<0,05). 

В случаях, когда в процессе МРТ клиническое состояние больных не из- 
менялось (8 наблюдений) не отмечено также и динамики хемилюми- 
несцентных показателей пероксидазной активности крови. 

При ухудшении клинического состояния (9 наблюдений) после 5 сеан- 
сов МРТ отчетливо проявлялось изменение параметров свечения и кинети- 
ки хемилюминесцентной реакции. При этом значительно уменьшалась ам- 
плитуда свечения, светосумма и увеличивалось время достижения макси- 
мального значения (L). 

Исследования пероксидазной активности крови показали, что при по- 
ложительном лечебном эффекте МРТ у больных с недостаточностью моз- 
гового кровообращения, очевидно, снимается синдром пероксидации. Об 
этом свидетельствует восстановление ферментативной антиоксидантной 
системы крови. Очевидно при этом в мембранах клетки активируются по- 
лиферментные комплексы, образующие цепи переноса электронов с вос- 
становлением кислорода. Следует отметить, что у больных с низкой перок- 
сидазной активностью и невысокой скоростью реакции наблюдается более 
высокий процент улучшения состояния после курса МРТ. 

Таким образом, анализ хемилюминесцентной информации, отра- 
жающей активностью пероксидазной системы, свидетельствует о том, что 
данный метод может быть применен для оценки эффективности лечения и 
в качестве прогностического теста при микроволновой резонансной 
терапии. 
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Изучение пероксидазной активности крови у больных церебральным атеросклерозом в процессе МРТ 


ВИВЧЕННЯ ПЕРОКСИДАЗНОЇ АКТИВНОСТІ КРОВІ У ХВОРИХ НА ЦЕРЕБРАЛЬНИЙ 
АТЕРОСКЛЕРОЗ В ПРОЦЕСІ МІКРОХВИЛЬОВОЇ РЕЗОНАНСНОЇ ТЕРАПП 


В.М. КУЗЬМЕНКО 


Проаналізовано результати вивчення пероксидазної активності крові в процесі мікрохвильової 
резонансної терапії у 41 хворого з ранніми формами цереброваскулярної патології обумовленої атеро- 
склерозом. Встановлено, що при позитивному ефекті МРТ у хворих з недостатністю мозкового кровоо- 
бігу знімається синдром пероксидації. Таким чином, визначення показників пероксидазної активності 
крові свідчить про те, що даний метод може бути застосований для оцінки ефективності лікування i B 
якості прогностичного тесту при мікрохвильовій резонансній терапії. - 


ИЗУЧЕНИЕ ПЕРОКСИДАЗНОЙ АКТИВНОСТИ КРОВИ У БОЛЬНЫХ ЦЕРЕБРАЛЬНЬМ 
АТЕРОСКЛЕРОЗОМ В ПРОЦЕССЕ МИКРОВОЛНОВОЙ РЕЗОНАНСНОЙ ТЕРАПИЙ 


В.М. КУЗЬМЕНКО 


Проведен анализ результатов изучения пероксидазной активности крови в процессе микровол- 
новой резонансной терапии у 41 больного сранними формами цереброваскулярной патологии обуслов- 
ленной атеросклерозом. Установлено, что при положительном эффекте МРТ у больных с 
недостаточностью мозгового кровообращения устраняется синдром пероксидации. Таким образом, 
изучение показателей пероксидазной активности крови свидетельствует о том, что данный метод мож- 
но применять для оценки эффективности лечения, а также в качестве прогностичнеского теста при 
микроволновой резонансной терапии. 


ЛИТЕРАТУРА 


1. Биохемилюминесценция в медицине и сельском хозяйстве.— Ташкент. — 1988.— 110 с. 

2. И.Ш. Весельский, А.В. Сонник Применение корректоров процессов перекисного окис- 
ления липидов и гиостаза в комплексном лечении больных с цереброваскулярными рас- 
стройствами. // Ж. невропатологии и психиатрии им. С.С. Корсакова.— 1997.— T. 97.2.— 
С. 51-54. 

3. Л.Ф. Дмитриев Биохимические аспекты атерогенеза: Роль антиоксидантов // Терапев- 
тический архив. - 1995.— T. 67, Ne12.— С. 73-77. 

4. В.М. Кузьменко, Б.П. Грубник, В.В. Кузьменко Микроволновая резонансная терапия 
больных с начальными проявлениями недостаточности мозгового кровообращения, ра- 
ботающих на химическом производстве / Тез. докл. междун. н/п конф. “Достижения в 
применении методов микроволновой резонансной терапии в Украине и за рубежом”. — 
Донецк.— 1996.— С. 14-15. 

5. В.М. Кузьменко, Б.П. Грубник Коррекция метаболических нарушений у больных с доин- 
сультными формами сосудистой патологии головного мозга при помощи микроволно- 
вой резонансной терапии // Физика живого.— 1997.— Т. 5, №1.— С. 115-119. 

6. Н.А. Мітряєва, М.А. Пиханова, Т.С. Бакай та ін. Адаптивні системи регуляції у ліквіда- 
торів наслідків аварії на ЧАЕС (за даними 7-річного спостереження) // Укр. радіологіч- 
ний журнал. - 1995.— T. 3.— С. 213-218. 

7. В.А. Нагорнев Атерогенез и иммунное воспаление // Архив патологии.— 1995.— T. 57.— 
№3.— С. 6-14. 

8. Я.В. Пишель, О.Е. Прокофьева, Н.Н. Дубенко и др. Значение перекисного окисления ли- 
пидов в патогенезе нарушений мозгового кровообращения при атеросклерозе. Актуаль- 
ные вопросы невропатологии, психиатрии и нейрохирургии / Тезисы докл. 1987.— 
С. 44-45. | 

9. Р.Р. Шакарашвили Свободнорадикальное окисление иантиокислительная активность в 
спинномозговой жидкости больных с острыми цереброваскулярными заболеваниями. 
Некоторые актуальные вопросы ангионеврологии. - Тбилиси. - 1984.— С. 98-107. 

10. Oxygen Radicals Chem. and Biol. Proc. 3 Int. Conf., Neuherberg, July, 10-15, 1983.— Berlin, 
New York, 1984.— 820 р. 

11. Free Radicals Ziver Injury. Proc. Int. Meet, Turin, June 27-29, 1985.— Oxford, Washington, 
D.C., 1985, 1020 p. 


119 


ІНФОРМАЦІЯ ДЛЯ АВТОРІВ 


Головним призначенням журналу є публікація статей та повідомлень, які демон- 
струють тенденцію переходу від традиційної біофізики до «фізики живого». Цей 
перехід грунтується на застосуванні синергетичних та квантових принципів в до- 
слідженні живого на всіх рівнях його ієрархічної організації. 

Журнал складається з трьох рубрикацій. У рубриці «Синергетика та квантова 
фізика живого» будуть друкуватися статті загального характеру з фундаменталь- 
них проблем синергетики, квантової механіки та їх можливого взаємозв'язку. Руб- 
рика «Біофізика» присвячується дослідженню живого методами традиційної біо- 
фізики. У рубриці «Медичні та біологічні застосування» будуть друкуватися статті 
прикладного характеру, які мають відношення до тематики журналу. 

Статті публікуються англійською, російською та українською мовами у відповід- 
ності з обраною в рукописі мовою. Обсяг статті має не перевищувати одного дру- 
кованого аркуша. 


ЗАПРОШУЄМО ДО СПІВРОБІТНИЦТВА! 


INFORMATION FOR AUTHORS 


The journal is specially designed to publish articles and reports that show the tendency of 
transition from traditional Biophysics to Physics of the Alive. This transition is based on 
the use of synergetic and quantum principles in studying life of all levels of its hierarchic . 
organization. 

The journal contains three main sections. Under the «Synergetics and Quantum Physics 
of the Alive» heading we will publish general articles on problems of synergetics, quan- 
tum mechanics and their possible interrelation as well as their application to the Alive. 
The «Biophysics» section will be devoted to the traditional issues of this science. Under 
the heading of «Medical and biological applications» we will publish articles of applied 
character, relating to the Journals subjects. 

According to the language of the original, articles will be published in English, Ukrainian 
and Russian. The size of an article cannot exceed I printer's sheet. All articles aré to be ac- 
companied by abstracts. 


YOU ARE WELCOME TO COOPERATION! 


ИНФОРМАЦИЯ ДЛЯ АВТОРОВ 


Главным назначением журнала является публикация статей и сообщений, выяв- 
ляющих тенденцию перехода от традиционной биофизики к «физике живого». 
Этот переход основан на применении синергетических и квантовых принципов в 
исследовании живого на всех уровнях его иерархической организации. 

Журнал содержит 3 основные рубрики. В рубрике «Синергетика и квантовая фи- 
зика живого» будут печататься статьи общего характера по фундаментальным 
проблемам синергетики, квантовой механики и их возможной взаимосвязи. Рубри- 
ка «Биофизика» будет посвящена методам традиционной биофизики в'исследова- 
нии живого. В рубрике «Медицинские и биологические приложения» будут печа- 
таться статьи прикладного характера, имеющие отношение к тематике журнала. 
Статьи будут публиковаться на английском, русском и украинском языках в соот- 
ветствии с языком рукописи. Объем статьи - до одного печатного листа. 


ПРИГЛАШАЕМ К СОТРУДНИЧЕСТВУ! 


Contents 


QUANTUM PHYSICS OF THE ALIVE AND ITS 
THEORETICAL FOUNDATIONS 


Парадигма повздовжніх електромагнітних хвиль в концептуальних 
основах фізики живого: критерії та умови формоутворення, 
інформаційно-діагностична значимість 

Ю.В. Човнюк, Т.М. Овсяннікова, Б.Ф. Рудько 


Оновлений підхід до теорії фізичних процесів в середовищах з пам'яттю 
Р.В. Луцик, Ю.Л. Ментковський 


Synaptic Transmission as the Cooperative Process 
A.V. Chalyi, A.N. Vasil’ev 


BIOPHYSICS 


Modulation of the Membrane Cunductance by Calcium Release Events 
in Vascular Myocytes D.V. Gordienko 


NMR Investigation of Sequence-Specificity of Daunomycin Binding 
with Deoxyhexanucleotides 
A.N. Veselkov, V.I. Pahomov, L.N. Djimant, R.J. Eaton, D.B. Davies 


Protein folding asa problem of theoretical biophysics 

A.P. Demchenko 

Влияние цианидных комплексов переходных металлов 

на электроповерхностные свойства и энергетические характеристики 
бактерий В.И. Подольская, Л.Н. Якубенко, 

З.Р. Ульберг, Н.И. Грищенко, А.И. Ермоленко 


Концентраційний потенціал Нернста з врахуванням 
термодифузійних ефектів С.П. Самко, О.В. Чалий 


Экспериментальное моделирование нелинейных эффектов 
взаимодействия переменных магнитных полей с живыми организмами 
М.В. Нецветов, П.К. Хиженков, В.И. Авсянкин 


MEDICAL AND BIOLOGICAL APPLICATIONS 


Микроволновая резонансная терапия в профилактике и лечении парезов 
желудочно-кишечного тракта после операций на толстой кишке 
АП. Кузьменко, И.Е. Соловьев, А.В. Topan 


Некоторье аспекть нарушения гомеостаза при критических состояниях 
и их коррекция с помощью КВЧ-терапии и ксеноимплантации 
эмбриональных мозгоспецифических тканей A.B. Малыхин. 


Изучение пероксидазной активности крови у больных церебральным 
атеросклерозом в процессе микроволновой резонансной терапии 
В.М. Кузьменко 


